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LIST    OF   SYMBOLS 


a  extrapolated    x    dimension  of  pile. 

b  extrapolated    y    dimension  of  pile. 

B    2  material  buckling  of  loaded  pile. 
m 

B  2  geometric  buckling  of  pile. 
G 

c  extrapolated    z    dimension  of  pile,   measured  from  source. 

d  extrapolated    z    dimension  of  multiplying  lattice. 

D  diffusion  coefficient. 

f.  fraction  of  source  neutrons  in  ith  energy  group. 

f  thermal  utilization. 

J  neutron  current,   neutrons  per  cm      -     second. 

k  multiplication  factor  in  an  infinite  reactor. 
oo 

k    „  effective  multiplication  factor  in  a  finite  reactor, 
eff 

L  diffusion  length  of  thermal  neutrons. 

L  lattice  diffusion  length. 

L  reflector  diffusion  length, 
r 

q  slowing  down  density,   neutrons  per  cm3  per  second. 

r.  Gaussian  range  of  neutrons  in  the  ith  energy  group, 
l 

Greek  Letters 

y  reciprocal  relaxation  length  of  the  mnth  mode  of  the  flux, 
mn 

6  reflector  savings. 

e  fast  fission  factor. 

rj  fast  neutrons  produced  per  thermal  fission. 


LIST    OF   SYMBOLS   (concl.  ) 

k  ratio  of  absorption  cross  section  to  diffusion  coefficient. 

A  transport  mean  free  path. 

S  macroscopic  absorption  cross  section. 

T  Fermi  age. 

$  thermal  neutron  flux,   neutrons  per  cm2  -  second. 

<(>  thermal  neutron  flux  corrected  for  harmonics  and 

end  effects. 


INTRODUCTION 

The  diffusion  length  of  thermal  neutrons  in  a  moderating  medium 
is  a  necessary  parameter  in  most  reactor  calculations.     The  deter- 
mination of  this  constant  is,  however,    subject  to  the  assumption  of  a 
particular  mathematical  model.     The  types  of  solutions  most  widely 
used  can  be  classified  in  two  broad  categories:     I)    One-group  dif- 
fusion, in  which  all  source  neutrons  are  considered  born     at  thermal 
energy.     In  this  case  the  source  is  treated  as  a  boundary  condition 
on  the  thermal  diffusion  equation.     Z)    Age -diffusion,   in  which  Fermi 
age  theory  is  used  to  describe  the  source  term  in  the  thermal  dif- 
fusion equation  as  the  slowing  down  density  from  a  source  of  high  energy 
neutrons. 

Discussions  of  one-group  diffusion  theory  using  specific  thermal 
source  boundary  conditions  can  be  found  in  several  references 
(1,    10,   13,    25).     Descriptions  of  age-diffusion   can  be  found  in  refer- 
ences (5,   17,   18).     A  paper  written  by  Wallace  and  LeCaine  (31)  pre- 
sents mathematical  solutions  to  the  neutron  diffusion  problem  in  several 
geometries  and  subject  to  several  different  source  considerations. 
Lee  (15),   of  Hanford,   discusses  techniques  designed  to  correct  for  the 
effect  of  various  inhomogenieties  in  the  moderating  medium. 

The  most  elementary  technique  used  in  the  determination  of  dif- 
fusion length  is  to  use  the  negative  reciprocal  of  the  slope  of  a  plot  of 
inc|>  vs.     z    for  data  taken  in  an  exponential,   or  sigma  pile.     This  sim- 
ple analysis  is  described  by  Hoag  (13)  and  assumes  an  infinite  moderat- 
ing medium  and  thereby  avoids  consideration  of  the  higher  harmonics 


of  the  flux  function.  Glasstone  and  Edlund  (II)  describe  a  harmonic 
analysis  procedure  whereby  the  higher  harmonics  may  be  taken  into 
account. 

Methods  used  in  this  work  include  two  different  assumed  thermal 
source  boundary  conditions;     a  method  whereby  consideration  of  the 
source  condition  is  avoided  completely,  and  a  method  which  calcu- 
lates a  series  of  constants  which  would  result  from  a  fictitious  thermal 
source  equivalent  to  the  existing  fast  source.     Two  approaches  based 
on  age-diffusion  were  used;     1)    the  assumption  of  a  monoenergetic 
fast  point  source,   and  2)  the  use  of  an  empirically  determined  Gaus- 
sian range  source  which  gives  a  closer  approximation  to  the  Pu-Be 
spectrum  than  a  monoenergetic  source.     Results  of  each  method  are 
compared,   and  the  extent  of  validity  of  the  approximations  discussed. 

Previous  work  at  Kansas  State  University  leading  to  this  work 
has  been  the  determination  of  the  Fermi  age  of  Fu-Be  neutrons  both 
from  theory  and  experiment  by  Steichen  (26)  and  the  empirical  fit 
of  a  sum  of  exponential  terms  to  the  slowing  down  density  in  the  KSU 
pile  by  Foulke  (9). 

THEORY 

One-Group  Diffusion  Model 

The  theory  of  thermal  neutron  diffusion  has  been  studied  by  many 
investigators  (I,    10,    II,    16).     The  following  paragraphs  will  serve  to 
outline  the  basic  theory  with  particular  attention  to  the  assumptions 
involved. 


The  change  in  neutron  density  in  a  volume  element  of  a  moderat- 
ing medium  is  the  result  of    (a)  the  net  flow  of  neutrons  through  the 
boundary  of  the  element,   (b)  the  number  of  neutrons  absorbed  by  the 
medium  per  unit  time,   and  (c)  the  production  of  neutrons  by  sources 
within  the  element.     Assuming  Fick's  law  to  be  valid,   the  neutron 
current,   defined  as  the  number  of  neutrons  per  second  flowing  through 
a  unit  area  normal  to  the  direction  of  flow,   can  be  expressed  as 

J  =    -D  grad  4>.  (l) 

The  net  flow  of  neutrons  into  a  volume  element  can  be  written 
in  terms  of  the  current  as 

-  div  J  =  div  (D  grad  <(>). 

Assuming  D,   the  diffusion  coefficient,   to  be  independent  of 
position,   it  can  be  factored  out  of  the  above  expression.      Thus 
the  leakage  term  can  be  written    DV^    and  the  equation  of  con- 
tinuity becomes 

dvz<!>   -  s  ♦  +  s  =-|t-  (Z) 

For  the  purpose  of  this  work,    only  steady  state  dif- 
fusion will  be  considered.     The  neutron  source  will  be 
treated  as  a  boundary  condition  in  thermal  source  theory 
so  that  Eq.   (2)  becomes 

DV!<f.    -    Sa$  =  0.  (3) 


s 

Dividing  both  terms  of  (3)  by    D    and  defining    -^ 
results  in  the  equation 


V2  <(>    -     k2  $  =  0  (4) 

23  , 

where  k2  =   _1_  =    -i_  . 

D  L2 

Conditions  not  previously  mentioned  which  are  necessary 

to  the  validity  of  elementary  diffusion  theory  include  the 

following: 

1.  There  are  no  collisions  between  neutrons,  which  means 
that  each  neutron  diffuses  independently  of  all  others.  Weinberg 
and  Wigner  (32),   and  Hughes  (14)  point  out  the  importance  of 

this  property  due  to  the  resulting  simplification  of  the  mathematics 
involved. 

2.  The  neutrons  diffuse  with  a  constant  (average)  energy 
and  no  energy  is  gained  or  lost  (on  the  average)  in  a  collision 
with  a  nucleus.  This  assumption  is  supported  by  Hughes  (14) 
who  states  that  experimental  evidence  indicates  that  the  Max- 
wellian  distribution  is  maintained  reasonably  well  throughout 
the  diffusion  process. 

3.  The  flux  is  a  slowly  varying  function  throughout  the 
pile,  with  no  sharp  dips  or  spikes.     This  assumption  is  essential 
to  the  validity  of  Fick's  law,  and  means  that  diffusion  theory  is 

not  applicable  in  close  proximity  to  concentrated  sources,  absorbers, 
or  the  boundaries  of  the  diffusing  medium.     It  is  also  necessary  that 
the  medium  be  only  slightly  absorbing,    since  a  high  absorption  cross 


section  would  cause  the  neutron  density  to  vary  substantially 
within  one  mean-free  path. 

4.     The  neutron  scattering  is  spherically  symmetric 
or  isotropic,   allowing  the  velocity  vector  to  be  treated  as 
a  scalar  in  the  equation  of  continuity.     Since  scattering  is  not 
isotropic  near  sources,   boundaries,    and  absorbers,   this  as- 
sumption necessitates  restrictions  similar  to  those  stated  in 
condition  (3). 

The  solution  of  Eq.  (4)  in  rectangular  geometry,  with 
the  origin  at  the  source  location,  is  derived  in  Appendix  A, 
and  given  here  as 

00         CO 

<Hx,y,z)=     )     /        A         cos  — SS-    cos — J-^    sinh  y       (c  -  z) 
'       Lj  L-  mn  a  b  rmr  ' 

m=l    n=l 

where    m    and    n    are  odd,    and  where  a,  b,   and    c    are  the 
extrapolated  x,   y,   and     z  dimensions  of  the  pile.     The  auxi- 
liary separation  constants  equation  is 


where 


Equation  (5)  is  derived  with  the  origin  at  the  center  of  the 

source  plane.     The  constants    A  may  be  determined  analy- 

mn 

tically,  through  the  use  of  a  particular  source  boundary  con- 
dition,  or  experimentally. 


-a    2  -  0  *+  7    2 

m        'n          mn 

=    K2 

m  ir 

•_  = ; 

m           a 

(5) 


(6) 


Point  Thermal  Source.      Numerous  investigators  (5,   6,    22,  23) 
have  used  the  assumption  of  a  point  source  of  thermal  neutrons  at  the 
location  of  the  physical  source  in  diffusion  length  experiments.     Other 
investigators  (3,   7)  have  used  several  point  sources  superimposed 
upon  one  another  to  describe  the  thermal  flux  in  a  diffusing  medium. 

The  application  of  the  point  thermal  source  condition  involves 
the  use  of  the  Dirac  delta  function  (11).     By  expanding  the  source  in 
a  series  of  orthogonal  functions  which  satisfy  the  boundary  condi- 
tions,    S  6  (x,   y)    at    z  =  0    can  be  written 


mux 

3  COS 

mn 
m=l      n=l 


_,,         .           \         \       _                  mix  niry  i-\ 

S6(x,y)    =    /        /       S        cos cos  '  (7) 

v    '  "         Li      Li        ran  a  b 


where    m    and    n    are  odd.     The    S        's    may  be  regarded  as 

sources  for  each  mode  of  the  flux.     To  obtain  the  values  of  the 

S       's,     each  side  of  Eq.  (7)  should  be  multiplied  by 
mn 

cos  X— cos    V  '         and  integrated  over  the  interval  of 

a  b  e 

orthogonality,   from    -  —    to    —       and to    —  .     By 

virtue  of  the  Dirac  delta  function,  the  left-hand  side  of  Eq.   (7) 
will  reduce  to    S.      Due  to  the  principles  of  orthogonal  functions 
the  right  side  will  become    S  times  one-half  the  orthogonality 

interval  in  each  of  the    x    and    y    directions,   giving 

S         =4*  (8) 

mn  ab  '    ' 

Since  only  positive  values  of    z    will  be  used  in  this 
work,  the  current  density  in  any  one  mode  will  be  equal  to  one -half 
the  total  number  of  neutrons  produced  in  that  mode.     The  current 


density  in  the    mn    mode  can  be  found  by  using  Fick 's  law,   Eq.  (l), 
and  the  expression  for  the  thermal  flux,   Eq.   (5),  to  give 

j         =.D(3Jmil_)  =    Dv       A        cos  BI£-  cos  EJLY-  cosh  y      (c-z) 

mn  v    9  z '  ta      mn  a  b  mnv         ' 

7=0 

(9) 
m ,   n    odd 

Setting  this  expression  equal  to  half  the  number  of  neutrons 
emitted  by  the  source  in  the    mn    mode, 

2S  mux  niry 

— | —     cos  cos  — t— '—    , 

at>  a  b 


results  in  an  expression  for    A 


mn 


(10) 


A 2S      -■ 

mn        ab  D  7       cosh  y       c 
mn  mn 

Substituting  Eq.   (10)  into  Eq.   (5)  will  give  an  expression 

for  the  flux  due  to  a  point  source  of  thermal  neutrons  located 

at  the  origin  of  the  co-ordinate  system 

^      V*  ,c                                                        sinh  r      (c-z) 

,,             »\\  2  S                 mux               niry              mn1          '         / , ,  x 

♦  (*.y.")=2,      Z  abbT       cos^—    co.-^I     cosfa                           (II) 

,           ,  mn                                                           mn 
m  =  I      n=I 

where    m    and    n    are  odd. 

If  a  single  point  source  is  located  at  some  arbitrary 

position  (x. ,   y. ,   0)  in  the    xy  plane,   then  the  integration  of 

a/Z     b/2 

SJ  J  «(x.,  y.Jcos™^-    cos^dxdy 

-a/2       -b/2 

mux,  niry. 

would  result  in    Scos — - —  cos — .— i-  .    Since    x.    and    y.    are 


constants,  these  cosine  terms  would  carry  through  the  analysis  so  that 
It,  niry 


2  S  cos 


ab  D  7        cosh  7        c 


For    N    sources  located  at  various  positions  throughout  the  xy  plane, 
superposition  of  fluxes  yields 


j\j  m  irx.  n  it  y. 

y2  S  cos cos  — r — 
5 ISZ_ 
mn 


i  =  l 


ab  D  7        cosh  y        c 
mn  mn 


An  equation  describing  the  flux  due  to    N  point  thermal  sources 
located  in  the    xy  plane  is  then 

N 


♦  ( 


*>y,z)  *2j   i_j   2. 


2  S 


n  Try. 


,    ab  D  7 
i  =  I     m  =  l     n  =  l  mn 


(12) 


sinh  7       (c  -  z) 
mix  niry  rmv  ' 

cos  cos  -^ '-    ■     ,, 

a  b  cosh  7         c 

mn 

where    m    and    n    are  odd. 

To  determine    k    from  Eq.   (5)  it  is  necessary  to 

determine    7I(    (or  any  other  of  the    7        s).      In  order  to 

accomplish  this  purpose  the  hyperbolic  functions  in  Eq.   (11) 

can  be  rewritten  as 


(13) 


00         00 

+(x,y,  z)  =2y    Z 


,     ab  D  7 
m=l    n=l  mn 


where    m    and    n    are  odd. 


2S  mirx  niry         'mn 

cos  cos     .    '     € 

a  b 


1-e 


-27       (c-z] 

ran1        ' 


.,    -27       c 
I+e        mn 


(14) 


It  is  desired  to  express  the  flux  simply  in  the  form 


i  (x,y,z)  =  C„  C 


-?„  z 


H     E 


(15) 


so  that  a  linear  least  squares  fit  of    i  n  <|>    vs.     z    will 
have  a  slope  equal  to  minus    ?,,  .      In  order  to  write  Eq.   (14) 
in  the  form  of  Eq.   (15)  it  is  necessary  to  factor  the  funda- 
mental mode,     m  o  n  =  1,  out  of  the  summation.     Defining 


1  -  e 


-27„(c  -  z) 


71  \x  c 


(16) 


Eq.   (14)  can  be  written 


♦**•*•»)  "aETTn,   cos-i-  cosV" 


-Tii  z 

a.         cE{i  + 


m=l    n=l      mn 


-y      z 


m  ix  niry 

cos cos        ' 

a  b 


1  -  e 


•2y       (c-z) 


T~.    :r2r      c~ 
1  +  e        mn 


_  irx  try 

C_.  cos cos     .  ' 

E  a  b 


■}         (I?) 


where    m    and    n    are  odd,  but  not  simultaneously  equal  to  one. 


Thus,     0       for  a  point  thermal  source  is  defined  by 


00    *£$ 


mix 
Tin  cos cos 


niry 


-zy    (c-z) 

mn1        ' 


I  +  e        mn 


(18) 


m=l    n=l 


_  irx  iry 

Tmn  CE  COS  T"    COB  T" 


where    m    and    n    are  odd,  but  not  simultaneously  equal  to  one. 


The  fundamental  mode  end  correction  factor,     C_,     and  the  com- 

E 

posite  end  and  harmonic  correction  factor,     C„,     for  the  point 


10 


thermal  source  boundary  condition  are  listed  in  Table  1,  along  with 
the  corresponding  correction  factors  for  all  other  boundary  con- 
ditions discussed  in  this  paper. 

The  correction  factors  derived  above  are  used  in  an  iteration 
procedure  for  determination  of  an  accurate  value  of    Yn    starting 
from  an  initial  estimate  of    Yn  .     The  initial  value  can  be  determined 
from  a  semi-log  plot  of  neutron  count  rates  vs.    z  .     The  steps  in  the 
iteration  procedure  ares 

1.  Using  the  initial  estimate,     (7,,)    ,   C_,    and    C„    are 
calculated  for  each  data  point  in  a  series  of  measurements  along  any 
vertical  axis  in  the  pile. 

2.  A  series  of  corrected  count  rates  is  calculated  accord- 
ing to           j,  <fr 

9corr         CE  CH 

3.  The  corrected  count  rates  are  used  to  obtain  a  new 

value,     (Tm)     ,     as  the  least  squares  slope  of  I  n  d>  vs.     z. 

2.  r  Tcorr 

4.  The  new  value,     (Y11)  •     ls  compared  with    (Yn)    , 

z  1 

the  value  used  to  calculate  the  correction  factors    C„    and    CTT. 

E  H 

If  the  difference  is  not  within  a  previously  determined  precision, 
(Y11)      replaces    (Yii)     in  step  one  and  new  correction  factors  are 
calculated  leading  to  a  third  value,(Yn)  •     Tne  precision  is  again 
checked  until  two  successive  values  of    Yn    are  within  the  desired 
precision. 

5.  When  a  final  value  of    Yn    has  been  obtained,  it  is  used 
in  Eq.  (6)  to  calculate  the  diffusion  length. 


11 


r 

g" 

rJ^ 

r^ 

^ 

H 

o 

o 

j 

i 

o 

-U 

u 

>s 

tf 

A 

g 

K 

^~^J 

^ 

* 

H 

V 

I 

i 

¥• 

% 

r^ 

PC" 

^c 

« 

i 
< 

I 

-S 

+ 

+ 

J 

i 

1 

J 

•ft 

i 

* 

t 

1 

+ 

s 

t 

» 

% 

(= 

N 

i 

f 

IL 

a 

3  + 

K 

i„ 

>; 

*» 

5 

— 

$ 

f 

~15 

+ 

a 

u 

i 

* 

i 

^ 

i 

*■ 

+ 
ft 

+ 

+ 

l0 

$ 

i 

i 

* 

0 

i 

i 

1° 

o 

Sh 

+ 

I 

s 

o 

< 

z 

He 

He 

Mb 

s 

g 

g 

a° 

Q 

+ 

ME 

■ 

." 

I. 

o 

« 

i* 

T{ 

4! 

o 

Jl 

Ul 

W  = 

4! 

WC 

W- 

< 

WE 

WE 

a 

c| 

+ 

% 

+ 

< 

El 

* 

_l 

Z 

>- 

< 

o 

K 

X 

Ul  u 

< 

_l 

< 

u 

C9  K 

zz 

< 

a 

u 

X 

0£ 

gs 

Z  3 
<  O 

30 

K 

H 

O  DC 

<e  co 

o 

WZ 

1- 

Z 

U 

0C3 

Is 

oo 
ecu 

1- 

P 

03 

o. 

«  E 

7) 

Z 

z> 

"  it 

o 

5 

o 

O 

o3 

<Z 

v> 

a. 

o 

a 

X 

e 

(D 

12 


An  IBM  650  program  for  the  determination  of  the  diffusion 
length  by  the  above  method  is  described  in  Appendix  C. 

Constant  Thermal  Source.       The  assumption  Of  a  constant 
source  of  thermal  neutrons  in  the    xy  plane  at    z  =  0    is  not  often 
used  since  it  is  an  obviously  poor  geometric  assumption  in  most 
cases.     It  is  included  here  primarily  for  purposes  of  comparison. 

To  determine  an  expression  for    A  based  on  a  constant 

r  mn 

source  of  thermal  neutrons  let    z  =  0    in  Eq.   (5),  and  let    (|>(x,y,  0)  =  <f>0. 

Thus-  00        00 

"<f>0=yy      A         cosH^cos^HIsinhr  c  (19) 

T0     Li    L         mn  a  b  'mn  v      ' 

m=I    n=l 
where    m    and    n    are  odd.      As  for  the  point  thermal  source, 
multiplying  both  sides  of  the  equation  by    cos    "  -       cos    jjl 
and  integrating  over  the  interval  of  orthogonality  yields 

^j.~i_                m+n+2 
4<t,0ab  ^—jj 


mti 


(-1)  =    A         sinh  Y         c  (-^2-). 

i  mn  mn      v   4    ' 


The  resulting  expression  for    A  is 

m+n+2 

A        =  16  oV0(-l)~-2— 

mn  »     ,   , 

m  n  ir    sinh  y        c 
mn 


(20) 


Substituting  Eq.   (20)  into  Eq.   (5)  gives  an  expression  for  the 
flux  due  to  a  constant  thermal  source  as 

oo       oo  m+n+2 

<Kx,y,z)=   J    J      LLigHOZ-,    cos^JL2Lcos^lsinhy      (c-z) 

V    Ll,     mn*Jsinh7       c  a  b  mn 

m  =  I    n=l  'mn 

(21) 

where    m    and    n    are  odd. 


13 


Using  the  same  procedure  outlined  in  the  section  dealing  with 
the  point  thermal  source  boundary  condition,    correction  factors 
C„    and    C        can  be  derived  from  Eq.    (21),    and  the  flux    expressed 

e  h 

in  the  form  of  Eq.   (15).     These  correction  factors  may  then  be  used 
in  an  iteration  procedure  identical  to  that  described  previously  to 
determine  the  diffusion  length  based  upon  the  assumption  of  a  con- 
stant thermal  source.     The  correction  factors    C_    and    Cu    for 

ht  rl 

this  case  are  listed  in  Table  1. 

An  IBM  650  program  for  the  determination  of  diffusion  length 
by  the  above  method  is  described  in  Appendix  D. 

Ratio     Method.      In  the  use  of  the  ratio  method,   no  source 
boundary  condition  need  be  assumed.     This  method  was  originally 
presented  by  Uhrig  (29)  in  1959. 

The  procedure  involves  the  application  of  orthogonality  con- 
ditions to  Eq.   (5).     Multiplying  both  sides  of  Eq.   (5)  by 
cos  cos    ,  <■     and  integrating  over  the  interval  of  orthogonality 

yields 

a/2       b/2 

r\  <Kx»y.zi)  cos  "TT"  cos  TT  dx  dy  =  A"  ~K  sinh  YnC0-2!  ) 

-a/2     J-b/2  a  b  4 

(22) 
where    Zi    is  held  constant.     If  the  left  side  of  Eq.   (22)  is 
defined  as    Fn  (z,)  ,   then 

a/2       b/2 

Fn(z,)=    \  \  (j^Xjy.Z]  )  cos cos  J-   dx  dy.  (23) 

J-a/2    J-b/2  a  b 


14 


F(z|)  can  be  evaluated  by  using  numerical  integration  techniques 

such  as  Simpson's  rule.     If  this  integration  is  performed  at  two 

different  levels    z,    and    z2    the  ratio  of    F,,  (z,)    to    Fu  (zz)  can 

be  written 

F„  (z,)  sinh  yu  (c  -  z,) 

Fn  (z2)     "    sinh  ylx  (c  -  zz)    ' 

By  expanding  the  hyperbolic  functions  in  terms  of  exponentials  the 
above  equation  becomes 

-2y„  (c-z,) 


F„  (z,) 
F„  (z2) 


-Tn(zi  -  zt) 


I  .  e"«Tii  lc-z2j 


Assuming  the  bracketed  term  in  Eq.   (24)  to  be  equal  to  one, 
a  first  approximation  of    ytI    is 

F,(z.  ) 


Yu 


in 


"tf-1  ) 


Using  an  initial  estimate  provided  by  Eq.   (25)  an  iterative 
procedure  based  upon  Eq.   (24)  can  be  set  up  to  solve  for    yu  . 
IBM  650  programs  for  the  evaluation  of    Fn  (z)    and  for  the 
iterative  solution  for    yu    are  described  in  Appendix  E. 

Double  Iteration  Method.      Previous  methods  described 
have  used  assumed  source  boundary  conditions  to  evaluate  the 
constants    ^mn    in  Eq.   (5).     In  the  ratio  method  described  here, 
the  evaluation  of  these  constants  was  avoided  completely.     The 

double  iteration  method  evaluates  the  constants    A  exoeri- 

mn        r 

mentally. 


(24) 


(25) 
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To  express  each  of  the  constants    A  in  a  form  in  which  they 

r  mn  ' 

may  be  readily  calculated  both  sides  of  Eq.  (5)  should  be  multiplied 
by  cost cos  .  i  a.nd  integrated  over  the  interval  of  ortho- 
gonality.    Thus, 

a/2       b/Z 

5a/2  ib/z  *  (x'y' z)  c°s ^  c°s  ^dx  dy  =  APq^  si"h  Ypq(c"z)- 

(26) 
Defining    F         as 

ra/z  rb/2 

F        =    J  J  <f>(x,y,z)  cos£^i    cos  q^Y  dx  dy, 

-a/2       -b/2 

the  constant    A        can  be  written 

pq 

4F 

A         =   ES. .  (27) 

pq  ab  sinh  y      (c-z) 

pq  ' 

F       can  be  evaluated,  using  numerical  integration  techniques, 
from  data  taken  at  several  points  at  a  given  elevation  in  the  pile. 
Hence,  it  is  possible  to  evaluate  any  constant    A         subject  to 
an  initial  assumption  of    7n  . 

The  double  iteration  process  may  be  described  as  follows: 

1.  The  values  of  the  constants    A        are  calculated 

pq 

using  an  initial  estimate  of    ytt. 

2.  Using  these  values  of    A      ,     correction  factors 
C_,    and    C„,  derived  in  the  same  manner  as  those  for  the  point 

£j  ri 

thermal  source  and  listed  in  Table  1 ,   can  be  calculated.     Thus 
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a  correction  factor  iteration  identical  to  that  described  previously 
can  be  used  to  determine  a  more  accurate  value  of    Tn. 

3.  The  value  of    yu    obtained  from  the  correction  factor 
iteration  process  is  then  compared  with  the  value  used  to  calculate 

the  constants    A      .      If  the  two  values  are  not  within  a  specified 

pq  r 

precision,  the  new  value  of    Tn    resulting  from  the  correction 
factor  iteration  process  is  used  to  calculate  another  set  of    A      's. 

pq 

4.  The  new  set  of   A      's  is  used  in  another  correction 

pq 

factor  iteration  process  to  obtain  a  third  value  of    Yn. 

5.  The  comparison  described  in  step  (3)  is  again  made, 
and  the  process  is  repeated  until  the  precision  check  is  satisfied. 

An  IBM  650  program  written  to  perform  the  double  iteration 
analysis  is  described  in  Appendix  F. 

Experimental  Pile  Size.      An  experimental  method  of  deter- 
mining the  effective    x    and    y    dimensions  of  a  sub -critical  assembly 
is  outlined  by  Babb,   et.   al.   (3).     This  method  involves  a  trial  and 
error  solution  for  the  effective  pile  size. 

The  variation  of  the  thermal  flux  along  the    x  axis  where  y  =  0 
and  z  =  constant,  is  described  by 


*  w  ■  2  An cos  ILF-  •  (2f}) 

n=l 

where    n    is  odd,  and  where  a  represents  the  effective 
x  dimension  of  the  pile.    Assuming  a  value  of    a,     least 
squares  values  of  the  constants    A      can  be  calculated, 
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along  with  the  error  squared  of  the  least  squares  fit.     The  best  value 
of    a    is  found  by  changing  the  initial  estimate  by  an  increment     Aa 
until  a  minimum  is  found  in  the  curve  of  error  squared  versus    a. 
An  IBM  650  program  written  to  perform  this  analysis  using  a  vari- 
able number  of  harmonics  is  described  in  Appendix  G. 

Age-Diffusion  Model 

In  all  of  the  previously  described  methods,  the  source  neutrons 
were  assumed  to  be  of  thermal  energy.     In  actual  practice,  however, 
sources  are  used  which  produce  neutrons  of  much  higher  energies. 
To  account  for  this  discrepancy,  a  combination  of  Fermi  age  theory 
and  thermal  diffusion  theory  is  used.     Fermi  (18)  used  age  - 
diffusion  theory  to  describe  the  thermal  flux  due  to  a  plane  source 
of  fast  neutrons  in  a  semi-infinite  medium.     Discussions  of  the  ap- 
plication of  age -diffusion  theory  to  the  determination  of  diffusion 
length  can  be  found  in  several  references  (5,7,    17). 

The  continuity  equation  for  steady  state  thermal  neutron  dif- 
fusion is 

VJ<f,    -    k*+    +    -§-    =    0.  (29) 

Assuming  a  point  source  of  fast  neutrons,   the  source  term 
in  Eq.  (29)  can  be  described  in  terms  of  the  thermal  slowing 
down  density,     qt,(r,r).     The  slowing  down  density  is  defined 
as  the  number  of  neutrons  being  slowed  to  thermal  energy  per 
unit  volume  per  second,   and  is  a  function  of  the  distance  from 
the  source  and  the  age  to  thermal  of  the  source  neutrons. 


18 


Substituting  this  quantity  into  Eq.   (29)  gives 

V2.).    -    k*<|>    +    -^ =0  (30) 

The  slowing  down  density,     <L,  »     *s  described  by  the  Fermi 
age  equation 

VJq(r,    r)    =iS^ji  (31) 

The  solution  of  Eqs.  (30)  and  (31)  for  an  infinitely  tall 
rectangular  parallelepiped  is  derived  in  Appendix  B 
and  stated  here  as 

00         00  2 

±i            \\\           **  rmrx  niry  mn     f,  , 

<t>(x,y,z)=  /      )     e        cos-——   cos-r-i-     e  <1  + 

m=l    n=l       mn 

erf  (-5—  -y        kTT  )}    +  e^""1*  fl  -  erf (-5 vy      ,/T")} 

2^—  J  l  2^         mn         JJ 


(32) 


where    m    and    n    are  odd.     Qois  the  source  strength  in 
neutrons  per  second  and    t  is  the  age  to  thermal  of  the 
source  neutrons. 

Monoenergetic  Fast  Source.    Equation  (32)  describes 
the  thermal  flux  due  to  a  point  source  of  monoenergetic  neutrons. 
This  relationship  can  be  used  in  two  different  ways  to  determine 
the  diffusion  length  (7,   17).     First  corrections  to  the  original 
C_    and    Cu    for  the  point  thermal  source  may  be  derived. 
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Second,  Eq.   (32)  may  itself  be  written  in  the  form  of  Eq.   (15) 

-Yn  z  ,      , 

$    =    CE  CHe  (15) 

The  second  and  more  direct  method  will  be  used  in  this  work. 

On  the  basis  of  Eq.   (32)  the  harmonic  correction  factor, 

C„,     can  be  derived  by  a  procedure  analogous  to  that  out- 

lined  for  the  point  thermal  source.     Since  Eq.  (32)  is  based 

upon  an  infinitely  tall  pile,     C„    cannot  be  derived  in  the  usual 

manner.     Davenport,   et.  al.   (7)   state  that  an  appropriate  end 

correction  would  be  the  same  as  that  used  in  thermal  source 

theory.     It  was  therefore  decided  to  include  the  terms 

-  2  y      (c-z) 
(I  -  e  ) 

in  Eq.   (32).     This  provides  an  end  correction  for  each  mode  of 
the  flux  such  as  would  result  from  replacing  an  exponential 
function  (which  approaches  zero  as    z    approaches  infinity)  by 
an  hyperbolic  function  (which  can  be  made  equal  to  zero  at  a  finite 
value  of    z).     Thus  the  boundary  condition  that  the  flux  go  to  zero 
at  a  finite  value  of    z    can  be  met  intuitively.     The  correction 

factors    C_    and    Cu    for  the  monoenergetic  fast  source  con- 

E  J.  1 

dition  are  listed  in  Table  I.  An  IBM  650  program,  written  to 
determine  the  diffusion  length  using  these  correction  factors, 
is  included  in  Appendix  H. 
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GauSBian  Range  Empirical  Source.  The  energy  spectrum  of 
Pu-Be  sources,  shown  in  Fig.  I,  indicates  that  these  sources  are 
definitely  not  monoenergetic.  This  fact  introduces  a  difference  in 
the  slowing  down  density  from  that  expected  from  a  monoenergetic 
fast  source. 

The  one  dimensional  solution  to  the  Fermi  age  equation  for 
a  monoenergetic  point  source  of  fast  neutrons  at  the  origin  of  an 
infinite  column  of  rectangular  sides    a    and    b    is  derived  in  Appen- 
dix B  and  given  here  as 

a  r\  .22.     22,  J     /  ni      .    a     »  z2 


q(z,T)= >       >       e  a  b 

abNRV^T    ^    ^, 
m=I    n=l 


where    m    and    n    are  odd.     Defining 

*~i  i_ 

,2  1,2 


n 


■<■=!    *=r) 


ab\M  it  t 

m=I    n=l 


the  expression  for  slowing  down  density  becomes 

z2 
q(z,   t)    =    Q0c  e  (33) 

for  a  monoenergetic  fast  source. 

The  vt    is  often  refered  to  as  the  slowing  down  length. 
The  quantity    2  nTt"    will  be  referred  to  as  the    Gaussian  range 
denoted  by 

r  =    2-sTF  (34) 


z\ 


UJ 


O    STEWART    EXPERIMENTAL 
DATA  (27) 

A     CURVE    FOR    EXPERIMENTAL 
DATA    AS    GIVEN   BY    STEWART 

B     CALCULATED      SPECTRUM 
OF     HESS  (12) 


10         II 


E,  MEV 

Fig.  I.   Calculated  and  measured  neutron  energy  spectrum  for 
a    Pu-Be  neutron  source. 
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Thus  Eq.   (33)  can  be  written  in  terms  of  the  Gaussian  range  as 

-zz/r* 
q(z,    t)  =Q0c  e  .  (35) 

A  plot  of  the  natural  logarithm  of  the  slowing  down  density 
vs.    z      for  a  monoenergetic  source  would  he  expected  to 
give  a  straight  line  of  slope    -  l/rz.     However,   in  an  experi- 
ment performed  by  Foulke  (9)  a  plot  of  the  natural  logarithm 
of  the  activity  of  cadmium  covered  indium  foils  (proportional 
to  the  slowing  down  density  to  indium  resonance)  vs.    zz    was 
shown  to  be  curved,   due  to  the  broad  energy  spectrum  of  the 
Pu-Be  sources.     This  curve  could  be  considered  to  be  due  to 
the  superposition  of  an  infinite  number  of  monoenergetic  sources. 
Hence  the  slowing  down  density  should  be  expressed  as 

f  -  z2/r.2 

q(z.    t)  =Q0   2     f;  c.     e  >  (36) 

i=  I 

where    f     represents  the  fraction  of  the  source  neutrons  in 
the    ith    energy  group,  and    r.     represents  the  Gaussian  range 
of  the    ith    energy  group.     Foulke  made  a  series  of  empirical 
fits  of  Eq.   (36)  to  a  plot  of  the  natural  logarithm  of  the  activity 
of  the  cadmium  covered  indium  foils  vs.    z2    for  the  KSU  Pu-Be 
sources. 

If  the  Gaussian  ranges  used  are  defined  from  source  energy 
to  thermal  energy,  then  Eq.   (36)  may  be  used  as  the  source  term 
in  the  thermal  diffusion  equation.     The  resulting  expression  for 
the  thermal  flux  would  be 
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f"    f<     V             Qo  k2  r.J/4  f  -  7      zr 

/             \\NN  i  mux  niry  mn  f, 

(x,y>z)=^    2,      Z         7 SFD     e  cos— j-cosAp    e  j  1  + 


i=l    m=l    n=I 


7        z  T 

7         r.  ,  'mn      r  Y        r.    ->1 

,  ,    z         mn     i  \)     ,  I,  c  i  z      i      mn    i  xl 

erf  (  t:  "  —2—  )J  +  e        I1  " erf  ("7T +  — —  m 

where    m    and    n    are  odd.     In  order  to  express  the  flux  in  the 

form  of  Eq.   (15),  the  above  equation  must  be  treated  similarly 

to  Eq.   (32)  to  derive  the  correction  factors    C_    and    C„. 

kj  H 

These  correction  factors  for  the  Gaussian  range  empirical 
source  are  given  in  Table  I.     An  IBM  650  program  for  the 
determination  of  diffusion  length  by  this  method  is  described 
in  Appendix  I. 


Determination  of  Material  Buckling  and 
Effective  Multiplication  Factor 


The  continuity  equation  for  neutron  diffusion  in  a  multiply- 
ing medium  differs  from  the  non-multiplying  case  by  the  addition 
of  a  neutron  production  term.  Thus  Eq.  (3)  for  steady  state  dif- 
fusion in  a  multiplying  medium  with  an  external  source  becomes 

DV2*    -  Sa<(>  +  k^  Sa  <(>  =  0  (38) 

where    k       is  the  multiplication  factor  for  an  infinite  reactor 
(no  leakage).     Equation  (38)  can  be  written  in  the  form 

V*  <f>    +  B    2  <p  =  0  (39) 


(37) 
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where  the  material  buckling,     B    2,     is  defined  bv 
°  m  ' 

The  solution  of  Eq.   (39)  in  rectangular  geometry 
and  subject  to  the  same  boundary  conditions  as  the  non- 
multiplying  case  is 
oo       oo 

♦<*.*■>«!    I    ^".ajSi-eo-afLrinhT^-. )  (41) 

m=l    n=I 

where    m    and    n    are  odd.     The  auxiliary  separation  constants 
equation  is  now 

-«r2-|32    -f    Y    2    =    -Bz.  (42) 

m       'n  mn  m  *■      ' 

If  a  point  source  of  thermal  neutrons  is  assumed  at 
the  origin  of  the  pile  co-ordinates,  Fig.  3  as  the  physical 
source  location,   the  constants    A  are  again  given  by  Eq.    (10) 

The  material  buckling  can  thus  be  calculated  by  the  same  pro- 
cedure as  was  K2      in  a  non-multiplying  medium. 

If  the  reactor  were  critical,  with  no  external  source,  the 
boundary  conditions  impressed  on  the  solution  for    <(i(x,y,z) 
would  be  that  the  flux  be  zero  at  all  extrapolated  boundaries. 
Thus,  for  a  critical  reactor  at  steady  state  (II), 

<t>(x,y,z)  =  A  cos^cos^cos^-  (43) 

where    a,  b,   and    d    are  the  extrapolated    x,  y,   and    z  dimen- 
sions of  the  critical  assembly.     The  reasons  that  only  one  mode 
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is  considered  here  are  (a)  that  criticality  can  only  be  achieved  in  the 
first  mode,   and  (b)  that  once  criticality  is  achieved  in  this  mode, 
the  higher  modes  die  out  exponentially.     The  auxiliary  separations 
constants  equation  for  the  critical  solution  is 

Equation  (43)  also  defines  the  geometric  buckling.     Hence 
it  can  be  seen  that  a  condition  for  criticality  is 

Bm   -    BG-  <44> 

Glasstone  and  Edlund  (II)  define  the  effective  multi- 
plication factor    k    ..,    as 

eff  -  B  2  T 

k       e         G 

k  „  =  — = (45) 

eff         I  +  L2  B  « 

where  the  condition  for  criticality  is    k   ,,  =  1.     Assuming 
criticality,   and  simplifying  the  exponential  term  gives 

k 


I   = 


(1+  L2B2)(1  +  tB2  ) 
Since    B2    is  generally  very  small,   further  simplification  yields 


I  +  B2  (L2+  t) 


(47) 


where  L  represents  the  diffusion  length  in  the  multiplying 
medium.  L  may  be  calculated  from  the  diffusion  length  in 
the  pure  moderator  by 

L?    =    L2  (I   -  f)  (48) 
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where    f    is  the  thermal  utilization  of  the  lattice. 

The  infinite  multiplication  factor,     k   ,     can  be  determined 

00 

from  Eq.   (47)  once    B   2  ,     t,     and    L.    are  known.      B   2    can  be 
'  m  t  m 

used  in  this  case  due  to  Eq.  (44)  and  the  assumption  of  criticality 

in  deriving  Eq.   (47).     Once    k      has  been  determined  it  can  be  used 

oo 

in  Eq.   (45)  to  determine    k   ,,    for  any  size  reactor  with  a  constant 

material  buckling  by  varying    B*  .      An  IBM  650  program  for  the 

G 

determination  of  the  material  buckling,     k     ,     and    k    ,,    by  the 

°         oo  eff       ' 

methods  described  above  is  included  as  Appendix  J. 

EXPERIMENTAL   FACILITIES 

General  Pile  Description 

The  Kansas  State  University  graphite  pile,   shown  in  Fig.   2, 
consisted  of  a  rectangular  parallelepiped,   68  in.   square,   and 
100  in.  high,   resting  on  a  concrete  foundation.     The  pile  was  con- 
structed of  machined  reactor  grade  graphite  blocks  approximately 
four  inches  in  cross  section  and  of  various  lengths.     In  stacking,  the 
long  dimension  of  the  blocks  was  alternated  90°  from  layer  to  layer. 
As  can  be  seen  in  Fig.   2,   certain  of  the  graphite  blocks  were  drilled 
through  their  entire  length  to  a  diameter  of  1.75  in. .      The  purpose 
of  the  holes  was  to  accommodate  fuel  elements  and  they  were  located 
to  provide  an  8-in.  lattice.     For  this  study,  in  the  effort  to  make  a 
solid,  homogeneous  moderator,   each  of  the  fuel  ports  contained  three 
graphite  cylinders  measuring  1.625  in.   in  diameter,   and  22.68  in. 
in  length.     The  cylinders  were  of  the  same  material  as  the  graphite  blocks 
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Fig.    2.     KSU  graphite  pile. 
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Fig.   3.    Schematic  diagram  of  KSU  pile. 
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The  graphite  used  in  the  construction  of  this  pile  was  purchased  from 
the  Great  Lakes  Carbon  Corporation,     The  blocks  were  machined 
from    R  -  IHLM    Nuclear  Grade  Graphite  with  a  thermal  neutron 
absorption  cross  section  between  3.  7  and  4.  5  millibarns. 

The  graphite  cylinders,  being  of  smaller  diameter  than  the  fuel 
ports,   left  crescent-shaped  air  gaps  of  0. 125  in.  between  the  tops  of 
the  cylinders  and  the  tops  of  the  fuel  ports.     The  blocks  along  the 
central  vertical  axis  of  the  pile  contained  horizontal  foil  slots, 
shown  in  Fig.   2,  with  a  cross  section  of  1 .  281  in.  by  0.  343  in.  . 
These  slots  could  accommodate  foil  stringers  as  well  as    BF,    or 
scintillation  probes.     The  density  of  the  solid  graphite  blocks  was 
1.683  g/cc.        Steichen  (26)  calculated  the  percentage  of  air  voids  to 
be  approximately  0.3  per  cent.     The  resulting  effective  density  of  the 
pile  was  1.678  g/cc. 

A  summary  of  the  physical  specifications  of  the  one  curie  Pu-Be 
neutron  sources  is  given  in  Table  2. 

Non-Multiplying  Configuration 

For  the  determination  of  diffusion  length,   a  solid,  homogeneous 
moderator  was  needed.     To  accomplish  this,  the  fuel  ports  were  filled 
with  graphite  rods  as  mentioned  above.     The  source  configuration  used 
was  a  five-source  cluster  at  the  center  of  the  source  plane  (positions 
5,    6,   7,    8,   and  9  in  Fig.    4).     In  this  configuration,   one  source  was 
located  at  the  exact  center  of  the  source  plane  and  the  other  four  were 
j.    four  inches  from  the  center  in  either    x    or    y    direction. 
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Table  2.     Summary  of  physical  specifications  of 
Plutonium-Beryllium  neutron  sources   (4) . 


Source         i  :  •   Neutron  emission* 

number         :     Grams  Pu       :       Grams  Be     s       rate   (n/sec) 


365 

7.87 

15.99 

1.54  x  106 

366 

7.87 

16.01 

1.73  x  106 

367 

7.86 

15.89 

1.82  x  10 

368 

7.86 

15.88 

1.69  x  106 

369 

7.86 

16.09 

1.71   x  106 

Recalibrated  September  1961   at  Mound  Laboratories 
during  Integrity  tests  and  recannlng. 


The  inner  container  of  each  source  is  0.85  in.  in  diameter, 
0.  90  in.   in  height,  and  is  made  of  tantalum.     The  outside  container  of 
each  source  is  1.02  in.   in  diameter,   1.30  in.  in  height,   and  is  made 
of  stainless  steel.     All  sources  are  sealed  by  welding.     The  melting 
point  of  tantalum  is  5,425°  F.      The  melting  point  of  18-8  stainless 
steel  is  2,600°  F. 

The  sources  were  calibrated  by  comparison  to  within    ±2  per  cent 
of  standards  calibrated  at  Los  Alamos  Scientific  Laboratory.     The 
absolute  accuracy  of  the  Los  Alamos  standard  is  reported  as    ±5  per 
cent,  thus  giving  an  absolute  accuracy  of    ±7  per  cent  to  the  sources 
listed  above. 
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ONE  OF  13  POSSIBLE  SOURCE  LOCATIONS 
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Fig.   4.     Source  locations  in  KSU  pile. 
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A  schematic  diagram  of  the  pile  is  presented  as  Fig.    3.     The  source 
plane  is  located  six  inches  above  the  concrete  foundation. 

Multiplying  Configuration 

In  this  configuration,   the  fuel  holes  contain  aluminum  tubes 
68  in.   long,   and  I.  312  in.   in  outside  diameter.     Each  tube  is  sup- 
ported in  the  fuel  port  by  three  aluminum  rings,   one  fixed  at  the 
center,  and  two  removable  rings  at  the  ends.     Eight  depleted  uranium, 
Savannah  River  type  MK  VII  fuel  elements,   as  pictured  in  Fig.    5,  are 
placed  in  each  tube.     Figure  6  illustrates  a  portion  of  the  loaded 
pile.     The  source  configuration  used  with  the  loaded  pile  was  the 
five-source  cluster  described  above. 

Neutron  Detection  System 

The  neutron  detection  system  used  throughout  this  work  consisted 
of  a  Nuclear  Chicago  Model  NC-202  BFj  probe,   Nuclear  Engineering 
Inventory  No.   97;     a  Nuclear  Chicago  Model  1062  preamplifier,  NEI 
No.   210;     a  B.   J.   Electronics  Model  DM1 -D  count  rate  meter,        NEI 
No.    363;     a  John  Fluke  Model  400  BDA  high  voltage  power  supply, 
NEI  No.    188;     a    Baird  Atomic  Model  132  scaler,  NEI  No.    144;    a 
Baird  Atomic  Model  960  timer,  NEI  No.    147;  an  Esterline  Angus  Model 
A.    W.    strip  chart  recorder,   NEI  No.    360;  and  a  traversing  mechanism, 
NEI  No.   568. 

The  BF3  probe  active  volume,  with  a  length  of  1/2  inch  and 
diameter  of  3/l6  inch,   contained    B10  Fj    gas  at  a  pressure  of 


Fig.    5.     Fuel  elements  for  KSU  pile 


6.     K.SU  exponential  pile. 
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70  cm  Hg  and  at  an  enrichment  of  96  per  cent.     The  Model  1062 
transistorized  preamplifier  had  a  length  of  6.  25  in.   and  a  diameter 
of  1.  5  in. .     The  BF3  probe  with  the  associated  counting  system  is 
shown  in  Fig.   7. 

The  traversing  mechanism,  illustrated  in  Fig.   8,  was  cap- 
able of  supporting  the  BFj  probe  in  a  variety  of  positions  in  the  pile. 
This  mechanism  was  designed  and  built  at  Kansas  State  University, 
and  a  detailed  description  of  it  is  provided  by  L.   R.   Foulke  (9). 

In  order  to  accommodate  the  BFj  probe  in  the  maximum  num- 
ber of  positions  and  with  the  minimum  amount  of  created  void  space, 
three  1.625  in.   diameter  graphite  cylinders,    shown  in  Fig.    9,   were 
drilled  along  their  central  axis  to  a  diameter  of  3/8  in. .     Their  use 
to  accommodate  the  probe  may  be  seen  in  Fig.    4. 

EXPERIMENTAL   PROCEDURE 

Cadmium  Shutter  Technique 

The  cadmium  shutter  technique  involves  the  use  of  a  cadmium 
shutter  imposed  between  source  and  detector  for  the  purpose  of  elimi- 
nating counts  caused  by  epithermal  neutrons.     Data  is  taken  without 
the  shutter,   and  also  with  the  shutter  in  place.     The  difference  be- 
tween the  two  sets  of  measurements  is  taken  (28)  as  proportional  to 
the  thermal  flux. 

The  cadmium  shutter  technique  is  generally  used  when  the 
detector  used  is  indium  foil.     The  same  purpose  can  be  accomplished 
in  this  case  as  when  cadmium  covers  are  placed  on  the  individual  foils  (28). 
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Fig.  9.   Special  graphite  cylinders  for  accomodation  of  neutron 
probes 
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The  cadmium  shutter  technique  is  often  preferred,   however,   since  it 
avoids  error  induced  by  overlapping  of  the  indium  and  cadmium 
resonances. 

When  individually  covered  foils  are  used,   the  cadmium  takes 
out  a  certain  portion  of  the  indium  resonance  nuetrons  before  they 
reach  the  foil.     Thus  the  covered  foil  count  rates  are  not  as  high  as 
they  should  be.     This  results  in  insufficient  correction  for  epithermal 
neutrons  in  calculating  the  cadmium  difference. 

With  the  cadmium  shutter  technique,   any  neutron  which  would 
be  at  resonance  energy  at  the  foil  elevation  would  be  at  a  consider- 
ably higher  energy  at  the  lower  elevation  of  the  shutter.     Thus  it 
would  be  passed  by  the  shutter  and  absorbed  by  the  foil.     In  this  way 
the  absorption  of  resonance  neutrons  by  the  cadmium  does  not  affect 
the  cadmium  difference  count  rates. 

There  is  another  effect,  however,  which  leads  one  to  believe 
that  the  use  of  the  cadmium  shutter  does  not  produce  accurate  results. 
This  is  due  to  the  fact  that  some  neutrons  are  thermalized  above  the 
shutter  but  below  the  elevation  of  the  foils.     These  neutrons,    since 
they  are  epithermal  at  the  shutter  location,  would  not  be  absorbed 
by  the  cadmium.     They  would  be  highly  absorbed  by  the  indium, 
however.     Since  the  activity  of  the  foil  with  the  cadmium  shutter  in 
place  is  supposed  to  be  proportional  to  the  indium  resonance  activity 
induced  in  the  foils  without  the  shutter,   one  can  see  that  any  activity 
induced  by  neutrons  thermalized  above  the  shutter  is  in  excess  of  that 
desired.     Thus  the  cadmium  difference  flux  calculated  in  this  manner 
will  be  too  low,   and  not  exactly  proportional  to  the  actual  thermal  flux. 
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Since  the  BF3  probe  used  in  this  work  had  been  shown  (9)  to  give 
count  rates  consistently  proportional  to  the  thermal  flux,   even  in  close 
proximity  to  the  source,   it  was  considered  adequate  for  determination 
of  thermal  flux.     Further  discussion  of  the  cadmium  shutter  technique 
can  be  found  in  references  (5)  and  (6). 

Vertical  Traverse  Measurements 

All  data  were  obtained  using  the  previously  described  BF3  probe 
and  associated  counting  system.     Measurements  were  made  in  the  foil 
stringer  slots  and  in  convenient  fuel  ports  using  the  specially  adapted 
cylinders  shown  in  Fig.    9.     All  measurements  were  made  with  the 
active  volume  of  the  BF3  probe  centered  on  the  vertical,    or  z  axis 
of  the  pile  by  means  of  the  traversing  mechanism. 

The  counter  was  operated  at  1400  volts  with  a  pulse  height  sensi- 
tivity of  0.8.  Some  tendency  toward  day  to  day  variation  in  counter 
sensitivity  was  encountered  so  that  normalization  counts  were  taken 
each  day  at  a  position  10  in.  above  the  source  plane  at  the  center  of 
the  pile.  Several  short  counts  were  taken  in  this  position  and  their 
average  used  in  the  normalization. 

Two  separate  vertical  traverses  were  made,   one  with  the  non- 
multiplying  configuration,   and  one  with  the  multiplying  configuration. 
In  each  case,  data  were  taken  such  that  approximately  10,  000  counts 
were  obtained  in  each  position.     This  was  done  with  10  series  of 
measurements  getting  1,000  counts  in  each  position.     Figure  3  is  a 
schematic  diagram  of  the  pile  showing  data  locations. 


Double  Integral  Mesh  Data 
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The  data  used  ito  evaluate  the  double  integral 


a/2       b/Z 
J  J  *(x,y,z) 

-a/2        -b/2 


niry 
b 


for  the  ratio  method  and  double  iteration  analyses  were  obtained 
by  running  horizontal  traverses  in  several  fuel  ports  at  a  given 
elevation  in  the  pile.     Equally  spaced  measurements  were  made 
along  the  center  fuel  port  from  the  center  of  the  pile  outward.     The 
same  procedure  was  then  repeated  for  each  of  the  remaining  four 
fuel  ports  on  one  side  of  the  pile.     The  space  between  successive 
data  points  was  determined  by  dividing  half  the  extrapolated  pile 
dimension  by  the  number  of  measurements  to  be  made.     The  value 
of  zero  at  the  extrapolated  boundary  was  included  as  an  additional 
data  point  in  evaluating  the  integral.     The  data  were  taken  in  three 
sets  of  1,  000  counts  per  position,    giving  a  total  of  3,000  counts  per 
position. 

Due  to  the  symmetry  of  the  pile,   it  was  considered  sufficient 
to  evaluate  the  integral  over  one-quarter  of  the  pile  and  multiply 
the  result  by  four.     Figure  10  illustrates  the  manner  in  which  these 
data  were  taken.     These  data  were  normalized  in  the  same  manner  as 
the  vertical  traverses. 
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Fig.    10.     Schematic  diagram  of  KSU  pile  illustrating 
method  of  obtaining  double  integral  mesh  data. 
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Horizontal  Traverses 

Horizontal  traverses  for  the  experimental  determination  of  the 
effective  pile  size  were  taken  such  that  approximately  5,  000  counts 
were  obtained  in  each  of  nine  positions  from  the  center  of  the  pile 
outward.     This  was  done  with  five  series  of  measurements  getting 
1,  000  counts  in  each  position.     Certain  of  the  graphite  blocks  were 
removed  and  inverted  so  that  it  was  possible  to  make  both  north-south 
and  east-west  traverses  at  elevations  of  38  in. ,   46  in. ,   and  54  in. 
above  the  source.     These  data  were  normalized  in  the  same  manner 
as  the  vertical  traverses. 

The  strip  chart  recorder  was  used  in  conjunction  with  the  count 
rate  meter  to  monitor  the  counting  during  periods  of  operator  absence, 
thus  insuring  that  no  erratic  data  were  recorded. 

DATA  PRESENTATION  AND  ANALYSIS 

Treatment  of  Raw  Data 

The  data  obtained  for  each  of  the  vertical  traverses  were  cor- 
rected for  background  and  are  listed  in  Table  3  in  the  form  of  counts 
per  minute.     The  deviation  reported  is  the  standard  deviation  of  the 
mean  in  each  case. 

The  double  integral  mesh  data  presented  some  difficulty  due  to 
low  count  rates  and  poor  statistics  near  the  outside  of  the  pile.     An 
attempt  was  made  to  correct  this  difficulty  through  a  statistical 
analysis  of  the  data  based  upon  the  equation 


Table  3.  Vertical  traverse  data  seta. 
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Distance  from 
source  plane 
(Fig.  3) 


BFj  Count  Rates 


Inches 


Non-Multlplylng 
configuration 


Multiplying 
configuration 


counts  per  minute   .  counts  per  minute 


20.000 
26.156 
28.000 

3009.8 
1906.2 
1621.4 

+ 

± 
± 

20.8 
15.8 
14.2 

30.156 
34.156 
36.000 

1400.0 
998.5 
841.1 

+ 

± 
+ 

12.7 
9.0 
8.1 

946.0 

+ 

7.7 

38.843 
42.156 
44.000 

737.3 
513,7 

434.5 

+ 
I 

± 

7.7 
2.6 
2.9 

756.0 
597.0 

+ 
+ 

4.1 
6.3 

46.156 
50.156 
52.000 

367.1 
262.4 
226.5 

+ 
+ 
+ 

3.3 
3.0 
2.3 

351.3 

± 

5.0 

54.156 
58.156 
60.000 

196.1 
138.8 
121.1 

+ 
+ 
± 

2.1 

1.7 
1.0 

218.3 

+ 

2.1 

66.156 
68.000 
74.156 

71.4 
62.1 
36.3 

+ 
+ 

+ 

0.8 
0.9 
0.4 

136.8 
79.4 

+ 
+ 

1.4 
1.1 

76.000 
82.156 
84.000 

31.2 
16.0 
12.6 

+ 
+ 
+ 

0.4 
0.2 
0.1 

44.4 

+ 

2.9 
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(x)  =  f 


m  a 

m=l 

where    m    is  odd.     The  constants    A        were  determined  by 
a  least  squares  analysis  of  the  data  in  any  one  of  the  five 
horizontal  traverses.     Each  experimental  value  was  then 
compared  with  its  least  squares  value,   and  if  the  difference 
was  greater  than  the  calculated  90  per  cent  confidence  limits 
for  that  point,  the  point  was  rejected.     All  rejected  points 
were  replaced  with  values  calculated  from  a  least  squares 
analysis  of  the  accepted  points.     An  IBM  650  program  written 
to  perform  this  analysis  is  described  in  Appendix  K.     The  pro- 
cessed data,   corrected  for  background,  appears  in  Table  4. 
The  deviation  reported  is  the  standard  deviation  of  the  mean 
in  each  case. 

Horizontal  traverse  data,  used  in  the  determination  of 
effective  pile  size,  was  corrected  for  background  and  listed 
in  Table  5.     The  deviations  reported  are  the  standard  devia- 
tions of  the  mean. 

Analysis  of  Model  Behavior 

For  reasons  stated  previously,   it  is  known  that  elemen- 
tary diffusion  theory  is  not  applicable  near  sources  or  boundaries. 
The  amount  of  error  induced  in  close  proximity  to  the  boundary  of 
the  diffusing  medium  is  discussed  quantitatively  by  Davison  (  8  ), 
in  terms  of  transport  theory.     Very  little  information  is  available 
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Table  5a.  Horizontal  traverse  data  In 
east-west  direction  at  three  elevations. 
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Distance  from 
center  of 
pile 


0 
4 
8 
12 
16 
20 
24 
28 
32 


BF,  Count  Rates 


z  =  38  Inches  J  z  =  46  Inches  !  z  =  54  Inches 


Inches    rcounts  per  mln.rcounts  per  mln.:counts  per  mln. 


859.9  i  22.4  394.3  ±  4.0 

808.7  i  22.2  387.5  ±  6.9 

779.7  i  21.5  363.5  ±  8.6 
690.6  *   14.8  333.7  J  3.2 

594.8  t     5.2  288.6  ±  6.8 
481.6  i     7.9  234.9 


351.1 

223.4 
104.7 


*     3.5        174.7 


5.9 

4.2 


i     4.0        113.2     t  3.9 
i     4.6         51.8     ±  4.0 


201.2  ±   1.6 

203.0  i  7.4 
187.7  t  2.4 
167.9  i  3.0 
150.7  *   1.3 

119.1  -  2.0 
92.1  *  1.4 
58.5  i  0.5 
26.3  *  0.3 
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Table  5b.  Horizontal  traverse  data  In 
north- south  direction  at  three  elevations. 


Distance 

center 

pile 

froi 
of 

•! 

BP,  Count  Rate 

IB 

J  z  =  3£ 

i  inches 

• 

|    z  =  46  Inches 

J   z  =  54 

inches 

In chef 

i 

J counts 

per  mln, 

i  > counts 

per  min. 

8  counts 

per  min. 

0 

723.6 

±    4.2 

373.5 

±  5.2 

202.1 

*  1.6 

4 

710.3 

±    9.4 

365.9 

±  6.4 

191.3 

±  3.5 

8 

666.6 

±  10.2 

339.0 

±  8.2 

185.3 

±  2.9 

12 

605. A 

±  11.5 

308.9 

±  8.0 

168.5 

±  2.0 

16 

526.4 

±    5.7 

269.2 

*  5.9 

14*. 6 

±  2.4 

20 

408.0 

±   4.9 

218.0 

*  5.3 

118.7 

±   1.0 

24 

308.6 

*    3.5 

166.7 

±  2.6 

90.7 

±  0.8 

28 

195.1 

±    2.9 

103.6 

±  2.2 

58.1 

±  1.0 

32 

79.9 

±   0.4 

47.6 

4    1.0 

26.4 

±  0.2 
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regarding  the  amount  of  error  induced  by  using  data  points  in  close 
proximity  to  the  source  in  a  diffusion  length  analysis. 

In  attempting  to  determine  the  effect  of  using  measurements 
made  near  the  source  to  determine  diffusion  length  two  sources  of 
error  must  be  examined!     (I)    The  possibility  that  a  given  assumed 
source  boundary  condition  might  lead  to  errors  near  the  source, 
and  (2)    the  possibility  that  the  thermal  diffusion  model  itself  may 
be  invalid. 

A  preliminary  problem  which  was  expected  was  inconsistency 
in  results  obtained  using  points  too  far  from  the  source.     The  primary 
reason  for  such  inconsistency  being  poor  counting    statistics  near  the 
top  of  the  pile.     For  this  reason,   a  series  of  determinations  of    L 
was  made  in  which  the  point  at  the  highest  elevation  was  dropped  from 
the  data  set  for  each  determination.      Consistency  in  the  values  ob- 
tained was  used  as  a  criterion  for  judging  the  highest  level  in  the  pile 
at  which  useable  data  could  be  obtained.      The  lowest  data  point  used 
in  this  analysis  was  at  least  two  diffusion  lengths  above  the  source. 
At  this  elevation  the  harmonic  content  is  negligible  (11),    so  that  the 
one-group  model  and  a  given  boundary  condition  could  be    assumed  to 
be  valid.     When  an  effective  upper  limit  had  been  established,  the 
behavior  of  the  model  close  to  the  source  could  be  studied. 

To  examine  a  given  assumed  boundary  condition  based  on  either 
the  one-group  diffusion  model  or  the  age-diffusion  model,  a  procedure 
was  used  whereby  the  effect  of  adding  additional  points  near  the  source 
could  be  observed.     Successive  determinations  of  diffusion  length  were 
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made,   each  using  one  additional  data  point  closer  to  the  source.     It 
is  logical  to  assume  that  as  long  as  all  points  used  were  far  enough 
from  the  source,   neither  the  assumed  boundary  condition  nor  the 
model  would  cause  any  deviation  in  the  values  obtained.     Hence  the 
apparent  deviations  of  these  values  of    L    should  be  random  about  the 
average  of  the  series  of  determinations. 

The  appearance  of  a  decisive  trend  in  several  determinations 
of    L,     either  increasing  or  decreasing,   was  taken  as  evidence  of 
error  induced  by  one  of  the  two  causes  listed  above. 

RESULTS   AND   CONCLUSIONS 

Experimental  Effective  Pile  Size 

The  experimental  determination  of  effective  pile  size  described 
previously  was  attempted  in  the  KSU  pile.  For  this  purpose,  several 
horizontal  traverses  were  made  and  these  data  were  processed  using 
the  IBM  650  code  described  in  Appendix  G. 

The  effective  size  in  each  of  the  horizontal  directions  was 
determined  at  three  elevations;     38,    46,   and  54  in.   above  the  source. 
Since  it  was  expected  that  the  harmonic  content  would  not  be  the  same 
at  all  elevations,   provision  was  made  in  the  IBM  650  code  so  that 
successive  analyses  could  be  performed  with  different  assumed  harmonic 
content. 

Table  6a  lists  the  results  obtained  in  the  east-west  direction 
(parallel  to  pile  void  channels)     with  various  assumed  harmonic  con- 
tents.     Table  6b  lists  the  corresponding  results  in  the  north-south 
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Table  6a.     Results  of  determination  of  experimental 
pile  size  in  the  x  direction,   actual  size 
68  in.   (theoretical  value    a  =  69.  375  in. ). 


Effective  Pile  Size     "a"    in  inches 

Eleva- 
tion 

Assumed  Harmonic  Content 

inches 

1. 

3 

I     1,3,5                !        1,3,5,7          :         1,3,5,7,9 

38  68.9±3.71      69.5±3.71         68.9±3.71  68.8  ±.3.71 

46  69.1  ±0.73      69.5±0.73         69.  9  ±  0. 73  69.7  ±0.73 

54  71.0  ±2.96      70.  1  ±  2.  96         71.  3  ±2.  96  70. 1  ±2.96 


Table  6b.     Results  of  determination  of  experimental 
pile  size  in  the  y  direction,   actual  size 
68  in.  (theoretical  value    b  =  69.  375). 


'.                                        Effective  Pile  Size     "b"    in  inches 

Eleva- . 
tion      '. 

Assumed  Harmonic  Content 

inches  • 

1. 

3 

:       1,    3,   5          :       1,3,5,7           :          1,3,5,7,9 

38  69.8  ±0.72      69.5±0.72         68.7  ±0.72  68.3±0.72 

46  69.9*1.34      70.0±1.34         70.5  ±1.34  68.9±1.34 

54  70.5  ±3.64      70.7  ±3.64         72.  0  ±  3.  64  70.3  ±3.64 


direction.     Deviations  reported  in  these  tables  were  calculated  from 
the  experimental  data.     Eight  values  of  pile  size  were  calculated 
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from  the  equation 


-I       4 

cos         -J- 


which  was  derived  from  Eq.   (28)  by  assuming  that  all  higher  modes 
were  insignificant.     The  value  of    <(>    at  the  center  of  the  pile  was  used 
as    A].     The  standard  deviation  of  the  mean  of  the  eight  values  of    a 
is  reported  as  the  expected  deviation  of  the  calculated  value  of  pile 
size.     Since  these  deviations  were  determined  directly  from  the  raw 
data  using  only  the  fundamental  mode,   no  difference  in  confidence 
limits  appears  with  variation  of  assumed  harmonic  content. 

The  data  in  Tables  6a  and  6b  are  much  too  random  in  nature  to 
accurately  define  the  effective  pile  size.      The  only  noticeable  trend 
is  the  increase  of  experimental  pile  size  with  increase  in  the  elevation 
at  which  it  was  determined.     This  would  indicate  that  the  pile  was 
shaped  similar  to  an  inverted  pyramid,   expanding  toward  the  top. 
However,   the  confidence  limits  placed  on  these  values  indicate  that 
a  more  extensive  analysis  should  be  performed  before  any  valid  con- 
clusions may  be  drawn. 

Since  definite  values  of  a  and  b  could  not  be  determined  ex- 
perimentally due  to  the  difficulties  mentioned  above,  the  theoretical 
extrapolation  distance  was  used  in  the  remainder  of  this  work.     This 

quantity  is  described  from  transport  theory  as     0.71   X,      ,   where     X 

tr  tr 

represents  the  transport  mean  free  path  in  the  medium.     The  value 
of    X        used  was  2.58  ±0.09  centimeters,    the   same  value  used  by 
Foulke  (9)  in  the  standardization  of  the  KSU  pile.     The  value  of  the 
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effective  pile  size  obtained  in  this  manner  was  69.  375  in.  in  both 
horizontal  directions. 

Diffusion  Length  Results 

Point  Thermal  Source.     The  results  obtained  by  assuming  a 
point  thermal  source  are  listed  in  Table  7.     The  first  four  deter- 
minations do  not  show  any  definite  tendency  to  increase  or  decrease, 
but  as  points  were  considered  closer  to  the  source  than  38  in.  ,   a 
definite  increasing  trend  became  apparent. 


Table  7.     Results  obtained  for  the  diffusion  length  based 
on  assuming  a  point  thermal  source. 


Elevation  of    '.  Number  of 

Lowest  Data   :     Data  Points  Used 
Point  ; 


Diffusion  Length 
p  b  1.678  g/cc 


inches 


50 

44 
42 
38 
36 
34 
28 
26 


10 
11 
12 
13 
14 
15 


53. 

32 

± 

0.87 

53. 

84 

± 

0.72 

53. 

71 

± 

0.  59 

53 

69 

± 

0.48 

54 

14 

± 

0.51 

54 

23 

± 

0.45 

55 

01 

± 

0.61 

55 

98 

rfc 

0.84 

The  above  values  were  calculated  using  the  1,   3,   5,  and  7 
harmonics  with  the  highest  data  point  76  inches  above  the  source. 
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The  type  of  behavior  displayed  by  these  calculations  was  taken  as 
evidence  that  the  point  thermal  source  boundary  condition  was  valid 
only  beyond  a  distance  of  approximately  38  in.   from  the  source. 
Probable  causes  of  this  phenomenon  are  the  fact  that  the  assumption 
disregards  the  high  source  energy  and  that  the  point  source  boundary 
condition  exaggerated  the  harmonic  content  and  would  therefore  be  in 
error  for  all  elevations  low  enough  for  the  harmonics  to  be  significant. 

Constant  Thermal  Source.       Diffusion  length  determinations 
based  on  an  assumed  constant  thermal  source  boundary  condition 
are  listed  in  Table  8.     Contrary  to  the  case  of  the  assumed  point 
thermal  source,  no  consistent  results  were  obtained  using  this  bound- 
ary condition.     Instead,   a  consistent  decrease  in    L    was  noticed  as 
additional  points  were  included  close  to  the  source,   and  the  average 
value  was  considerably  lower  than  that  obtained  using  any  other  method. 


Table  8.     Results  obtained  for  the  diffusions  length  based 
on  an  assumed  constant  thermal  source. 


Elevation  of 
Lowest  Data 
Point 


Number  of  Data 
Points    Used 


Diffusion  Length 
p  =  1.678  g/cc 


inches 


50 
44 
42 
38 
36 
34 
28 
26 


8 

9 

10 

11 

12 
13 
14 

15 


50.38  ±    0.58 

50.17  ±   0.45 

49.62  ±   0.47 

48.92  ±   0.52 

48.78  ±   0.45 

48.37  ±   0.46 

47.86  ±   0.47 

47.58  ±   0.44 


The  above  values  were  calculated  using  the  1,   3,   5,   and  7 
harmonics  with  the  highest  data  point  76  in.   above  the  source.     The 
conclusion  was  drawn  from  these  results  that  since  the  constant 
thermal  source  boundary  condition  was  not  at  all  compatible  with 
actual  conditions,  and  since  the  results  obtained  by  its  use  showed 
none  of  the  desired  consistency,  that  this  boundary  condition  should 
be  considered  completely  invalid.     This  was  not  unexpected  in  that 
this  particular  approximation  disregards  the  high  source  energy 
and  also  uses  an  incorrect  geometric  approximation  of  the  source. 
This  boundary  condition  was  used  for  comparative  purposes,   and 
to  determine  the  effect  of  an  erroneous  boundary  condition  on  dif- 
fusion length  calculations. 

Ratio  Method.      Evaluation  of  the  diffusion  length  by  the  ratio 
method  was  accomplished  using  integrals  evaluated  at  four  different 
elevations.     These  integrals  were  used  in  six  separate  pairs,  pro- 
viding six  independent  determinations  of  diffusion  length.     The 
results  of  these  determinations,   at  a  pile  density  of  1.678  g/cc, 
are  given  in  Table  9. 
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Table  9.     Results  of  the  ratio  method  for 

determination  of  diffusion  length. 


z,    : 

zj         : 

Ratio      F 

/f2 

\     Diffusion  Length  ' 

inches  \ 

inches  : 

[                  cm 

28 

20 

0.584 

126.52 

36 

28 

0.544 

61.69 

44 

36 

0.532 

55.36 

36 

20 

0.  318 

79.19 

44 

28 

0.289 

58.29 

44 

20 

0.169 

68.48 

1    Confidence  limits  not  calculated  due  to  erratic  nature  of  the  results. 


It  is  immediately  obvious  that  the  range  of  these  results  is    ex- 
ceptionally wide.     In  fact,  the  difference  between  the  lowest  and 
highest  values  is  71.12  cm. 

The  basic  difference  between  this  method  and  all  others  is  the 
absence  of  a  vertical  traverse.     Since  the  diffusion  length  is  known 
fundamentally  to  be  strongly  related  to  the  attenuation  in  the 
z  direction,   it  is  logical  to  assume  that  a  larger  number  of  measure- 
ments at  different  elevations  would  serve  to  establish  the  value  of    L 
more  firmly.     The  ratio  method  uses  only  two  z's,   and  the  measure- 
ments made  at  these  elevations,  the  integrals,   cannot  be  calculated 
to  the  high  degree  of  accuracy  necessary  because  of  the  poor  count- 
ing statistics  near  the  edge  of  the  pile. 

Whereas  the  ratio  method  appeared  desirable,   since  it  avoided 
the  necessity  of  imposing  an  erroneous  boundary  condition,  the 
consideration  described  above  and  the  actual  results  obtained  indicate 
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that  it  is  not  a  practical  method  for  the  determination  of  the  diffusion 
length. 

Double  Iteration.       Diffusion  length  results  from  the  double 
iteration  process,  at  a  pile  density  of  1.678  g/cc,  are  listed  in 
Table  10.     All  of  these  results  were  obtained  using  the  1,    3,   and  5 
harmonics.     Horizontal  traverses  for  the  evaluation  of  the  con- 
stants   A  were  taken  at  four  different  levels  in  the  pile.     Hence 
mn  r 

there  are  four  sets  of  values,    each  based  on  constants  determined 

at  a  different  elevation. 

It  can  be  seen  from  Table  10  that  the  results  obtained  using 

constants  evaluated  at  the  higher  elevations  do  not  display  the  same 

consistency  as  those  obtained  using  constants  evaluated  at  the  lower 

elevations.      This  can  be  explained  as  follows:      The  value  of    Yn  , 

and  hence  the  value  of  L,   obtained  from  this  process  is  the  direct 

result  of  a  least  squares  analysis  of  a  plot  of  i  n  c|>  vs.     z. 

^       ,     is  determined  as  explained  previously  by  essentially  removing 

the  higher  harmonics  from  consideration  by  means  of  the  correction 

factor,     C„,     derived  previously.     In  order  to  do  this  accurately,   it 
ri 

is  generally  advisable  to  include  as  many  harmonics  in  the  analysis 

as  is  feasible.     This  is  particularly  true  when  deriving    C       for  a 

n 

point  close  to  the  source  where  the  harmonics  are  known  to  be  of 
significance. 

However,  in  determining    A  for  the  higher  harmonics  at 

z  =  44  in. ,  for  example,   one  is  attempting  to  evaluate  a  quantity  which 
is  numerically  insignificant.     This  is  known  to  be  the  case  because  at 
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Table  10.     Diffusion  length  results  from  the  double  iteration 
technique  using  the  1,    3,    5  harmonics. 


Elevation  of 
Lowest  Point 


Diffusion  Length,     cm 
Elevation  at  which  Constants  were  Determined,   inches 


inches 


20 


28 


36 


44 


50 

53. 

34  ±1. 

64 

53. 

35  ±1. 

65 

53. 

58  ±  1. 

68 

46 

53. 

54  ±1. 

26 

53. 

55  ±1. 

28 

53. 

83  ±1. 

29 

44 

53. 

96  ±1. 

06 

53. 

98  ±1. 

05 

54. 

31  ±1. 

09 

42 

53. 

67  ±0 

88 

53. 

68  ±0. 

88 

54. 

05  ±0. 

90 

38 

53. 

37  ±0 

72 

53 

.  38  ±0 

.72 

53 

81  ±0 

.74 

36 

53 

75  ±0 

67 

53. 

76  ±0. 

67 

54. 

27  ±0. 

70 

34 

53 

65  ±0 

58 

53 

65  ±0 

58 

54. 

23  ±0. 

60 

30 

53 

63  ±0 

49 

53 

63  ±0. 

50 

54. 

33  ±0 

52 

28 

54 

23  ±0 

59 

54 

28  ±0 

59 

55 

14  ±0 

67 

26 

54 

71  ±0 

60 

54 

70  ±0 

58 

55 

73  ±0 

71 

54.47  ±1.80 
54.88  ±1.40 
55.52  ±1.21 
55.37  ±0.99 
55.34±0.82 
56.05  ±0.83 
56.20  ±0.73 
56.62  ±0.68 


such  distances  from  the  source  the  harmonic  content  is  negligible  (11). 
The  numerical  integration  process  used  in  this  analysis  was  found  to 
be  overestimating  the  magnitude  of  these  harmonics  at  higher  ele- 
vations.    This  was  believed  to  have  been  caused  by  the  efforts  of  the 
numerical  integration  process  to  give  credence  to  the  statistical  noise 
inherent  in  the  data.     At  lower  elevations,   however,  where  the  harmonic 
content  was  appreciable,  the  estimate  was  considerably  better.     This 
statement  is  supported  by  Fig.   11,  which  illustrates  the  predicted 
relative  contribution  of  the  13  and  31  harmonics  as  a  function  of  dis- 
tance from  the  source.     Analytical  curves  based  on  the  point  thermal 
source  and  constant  thermal  source  boundary  conditions  are  included. 
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The  fast  source  used  can  presumably  be  represented  by  an 
equivalent  thermal  source  which  will  be  neither  a  point  nor  constant, 
but  some  intermediate  configuration.     It  is  therefore  logical  to 
assume  that  the  actual  harmonic  content  of  theflux  function  would 
be  something  in  between  these  two  specific  thermal  boundary  con- 
ditions.    It  can  be  seen  from  Fig.   11  that  as  lower  elevations  are 

used  to  determine  the  constants    A        ,     the  predicted: harmonic  con- 

mn 

tent  is  bounded  by  the  point  thermal  and  constant  thermal  source 
predictions. 

Further  evidence  of  the  effect  of  harmonic  cpntent  on  the  be- 
havior of  the  double  iteration  technique  is  presented  in  Figs.    12, 
13  and  14.     Figure  12  is  a  plot  of  the  data  given.in  Table  10,  and 
shawS'the  behavior  of  this  technique,fQr  all  four  elevations  when  the 
1,    3,   and  5  harmonics  are  used.     Figure  13  shows  the  effect  of 
neglecting  the  5th  harmonic  with  constants  determined  at  elevation 
z  =  44.     Figure  14  shows  this  same  effect  plus  the  effect  of  taking 
the  7th  harmonic  into  consideration.     The  data  used  in  these  plots 
are  presented  in  Table  11.     It  is  interesting  to  note  that  the  effect 
of  using  additional  harmonics  in  the  analysis  is  much  more  pro- 
nounced at  the  higher  elevation. 

The  conclusions  to  be  reached  from  the  observations  mentioned 

above  are  that  the  double  iteration  technique  gives  valid  results  only 

when  the  constants    A  are  determined  from  data  taken  at  a  level 

mn 

in  the  pile  where  the  harmonics  contribute  significantly  to  the  total 
thermal  flux. 


64 


H 

T 

•l 

* 

a 

Li 

^ 

?• 

O 

S 

E 

LI 

H 
U 
Q 

> 

Ul 

-J 

UJ 

</> 

o 

t- 

^ 

<r 

p 

OT 

IT 

O 

O 

</> 

UJ 

x 
o 

■z 


< 

u 

t   a: 
ui 

> 

O     u. 

*   o 


10  o 
o. 

«,  £ 

ro    5 

o 


to    ,_ 
tu 

N    -, 

o 

CO 

CO 

evi    2 

o 
c 

U>    "- 

~  u 
o 

z 


UJ 

a: 
< 

CO 
H 

z 

O 
0. 

CO 
< 

CO 
CO 

< 
< 


o 


§ 


Q 

q 

q 

m 

* 

IO 

in 

in 

in 

IU 


UJ 
UJ   -r 

m 


o: 
o 

> 
< 

x 

UJ 
CO 


CM  O 


O 


no  'i 


65 


It 

o 

1            1 

■ 

> 

<\J 

< 
Ill  "* 

K 

m 

z 

ui      m 

CD  ♦ 

z  io  io 

1 

!  n 

l    _ 

1. 

O 

to 

Ul 

_  '1 

i 

/ 

I    "^ 

I 

Z  m 

o  _'  _" 

X 

O 

o 

u 

h- 

-o  °   • 
z 
o 

CO 

z 

•> 

Z« 

.2 

I 

<D 

Hi 

2  Z 

Q* 

1 

/              ' 

sc 

(0 

2 

-1 

°  5 

o  a 

!  a       I 

/ 

< 

o: 

1 

°l 

o£ 

Cf 

Z  W 

- 

i-i/1 

J 

CM 

Ui 

> 

o  o 
SE  en 

i 

1 

u. 

K  H 

" 

/      _         i 

1 

o 
eg 

o 

h- 

z 

o 

<  z 

x<t 

1 52 
z  z 

UJ  o 

1 

1      / 

[Q 

Q. 

< 

cc  o 

/ 

■ 

t 

10 
IO 

i 

Kg 

1"  ' 

/ 

1 

a 

u.  ? 

1 

1 

10 

10 

CC 

o£ 

1 

1 

/ 

U-  Q 

u. 

oo 

~                         J 

CM 

ro 

o 

X 
Z  h- 
O  Ld 

,  / 

| 

- 

O 

IO 

UJ 

o 

F2 

1           / 

I 

CL 

3 

S  z 

u 

O 

3  O 

■   / 

1 

a/          ' 

*j 

CO 

cm 

to 

CO  y- 

B 

s 
o 

V>  < 

/ 

c» 

<  CC 

1 

/       1 

CD 

CO 
CM 

u. 

.,    W 

fc  1- 

(0 

o  — 

II 

'J- 

o 

ttH 

Q. 

CM 
CM 

z 

I 

to 

o  -J 

LU  CO 
U.  3 
U.  O 

1          1 

1 

I             I 

1 

I               1 

I 

CM 
O 

o 

UJ  Q 

IO     Q 

o 

s 

Q 

Q 

Q 

C 

5N 

N 

m 

* 

K) 

(VI 

O 

m 

in 

in 

in 

in 

1 

■» 

WO'l 


o 

Q 

o 

o 

f>-' 

m 

n 

* 

in 

10 

m 

a 

NO  'T 


66 

CJ 

10 

z 

I 

o  a 

o 

UJ 

in 

1-  H 

Z  UJ 

iii 

a 

■9 

H 
Z 

(0 

U 

O 

h- 

to    u 

o 

z 

*   z 

< 

*  s 

o 

CO 

z  z 

*  <J) 

o 

o 

so 

S3 

< 

o 

H 

T 

z 

<T 

?2 

(O 

H  = 

U- 

z 

Q 

O 

LJ 

o 

X 

9  . 

10  1- 

ae 

z 

UJ 

1- 

o 

u. 

Ul 

<D     «■ 

u. 

2 

i- 
< 

a 

z 

ro    |_ 

u. 
o 

o 

cc 

z 

< 

CM    U. 

o 

DE 

I*>    — 

UJ 

O 

1- 

H 

u_ 

— 

O    £ 

1 

(o   o 

UJ 

K 

3 
O 

to 

ffl 

8  M 

2 

< 

3 
O 

O 

u. 

o 

Q    10 

ro 

111 

H 

°N 

o 

*  z 

<M    S 

o 
III 

K   1- 

»- 

u. 

O  < 

(0 

u. 

•> 

CM    o 
N 

.  i1-  Z 

O 

* 

UJ   - 

m  5 

« 

u. 

67 


a 

i 

o 

<d 

■ 

01 

01  -p 

3 

a 

id 

-t-> 

Ft 

c 

8 

o 
o 

<M 

a 

■aw 

a) 

a 

a 

o 

•H 
01 

i 

■p 

Hi 

."£1 

O 

+3 

co 

a 

■P 

ai 

rH 

t, 

3 

<ii 

(0  <t-i 

O  <H 

h 

■H 

n 

3 
* 

3 

•H 

a) 

s 

r4 

Tt 

a 

O 

o 

5 

en 

(1) 

1 

S 

V-i 

a 

■H 

o 

Ot4 

*J 

nl 

• 

^ 

•— 

01 

«— 

■p 

•H 

0) 

■H 

^1 

01 

E-> 

o 

■H    (0 
4->    0) 

0)43 
3  2 

H  C 
CO  «H 

> 

0  - 

c 

01  O 

+>«4 

0  +> 
o)  «J 

43   > 

01  01 
BH 

O  01 

o 


c 
o 

ag> 

<«   01 


o 

5  to 
a 

<H  01 
■HH 
Q 


a 

O 

■**  43 
KV    co  4-> 

C 

%4    01 

a 


Cm    01 


a 
o 

2  c 

<H    01 


CO 
01 

IS 
o 

o  o 
p. 
n 

o  o) 

<->    0) 

S» 

CD 

H 


oo 
vo 

CM 

Ov 

o 

o\ 

t- 

CM 
VO 

8 

o\ 

*■ 

•" 

•" 

o 

o 

o 

o 

o 

o 

41 

41 

41 

+i 

41 

41 

41 

41 

41 

i 

i 

in 

CM 
CO 

• 

in 

CM 

in 

o 

• 

in 

CO 

m 
in 

in 

• 
in 

■sf 

in 

VO 

• 

in 

rn 
o 

VO 

in 

i 
i 

CO 

ON 
CM 

ov 
o 

o 

o 

s 

CM 

in 

t- 

VO 

c^ 

•"I 

~ 

— 

o 

o 

o 

o 

o 

o 

o 

41 

41 

m 

41 

41 

41 

+i 

41 

41 

•H 

co 

in 

m 
in 

CO 

m 
in 

tn 

-* 
in 

in 
o 

in 

oo 

in 

C- 

in 

CM 

in 

tn 

■4 

in 

m 
in 

KV 

in 
in 

CM 

in 

ov 

CM 

ov 
o 

ov 

t- 

Ov 
vo 

8 

CM 

in 

m 

VO 

CO 

— 

'- 

- 

o 

o 

o 

O 

o 

o 

o 

* 

■H 

41 

41 

41 

41 

41 

41 

41 

•H 

CO 

in 

CM 

co 

o\ 

CM 

o 

in 
t- 

OV 

m 

t- 

Ov 

5 

in 
in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

o 

co 

§ 

CM 

Ov 
OV 

CM 
00 

CO 

to 

t- 

CO 

VO 

- 

- 

•- 

o 

o 

o 

o 

o 

41 

41 

41 

41 

41 

41 

41 

41 

i 

I 

in 

00 
00 

in 

CM 

in 

in 
in 

t- 

m 

in 
in 

1*1 

in 
in 

in 
o 

VO 

in 

3 

VO 

m 

CM 
VO 

• 
VO 

in 

i 

i 

i 

I 

o 

m 

- 

CM 
Ov 

VO 

t- 

CM 

vo 

in 

CO 
VO 

CVJ 

- 

- 

- 

o 

o 

o 

o 

o 

o 

o 

41 

41 

•H 

41 

+i 

41 

-H 

41 

+1 

■H 

VO 

in 

in 

O 

^* 
in 

KV. 

in 

■* 
in 

t- 
CM 

in 

o 

in 

m 
in 

• 

in 

o\ 

in 

in 

in 
in 

CO 
Ov 

in 
in 

o     <o 

in     •* 


CVJ 


» 


VO       -=*• 

tn     m 


R 


VO 
CM 


68 


Using  the  values  of    A  determined  at  elevation  z  =  20  in.  ,   the 

magnitude  of  the  fundamental  and  the  higher  harmonics  at  the  source 
plane  could  be  calculated.     Thus  an  equivalent  thermal  source  condi- 
tion can  be  described  based  upon  the  actual  flux  distribution  in  the 
pile.     The  constants    A  used  to  calculate  this  equivalent  source 

are  listed  in  Table  12.     The  equivalent  source,  with  calculated  flux 
plots  at  distances  of  10,   20,    30,   and  40  in.   above  the  source,   is 
plotted  in  Fig.   15.     A  similar  plot  is  made  for  the  point  thermal  source 
boundary  condition  in  Fig.    16.     The  double  iteration  equivalent  bound- 
ary condition  was  plotted  using  the  11,    13,    31,   and  33  harmonics.     As 
indicated  on  the  plot,   these  curves  fit  the  normalized  data  accurately. 
It  is  believed  that  this  equivalent  thermal  source  condition  could  be 
duplicated  within  limits  of  experimental  error  using  constants  deter- 
mined at  elevation    z  =  28  in.    and  possibly  also  using  constants  evalu- 
ated at    z  =  36  in. 


Table  12.     Values  of  the  constants    A  evaluated 

by  the  double  iteration  technique  at 
elevation    z  =  20  inches. 


m  : 

n   : 

mil 

I 

1 

15.70 

0.  040 

1 

3 

7.59 

0.002  x  10"2 

3 

1 

,7.59 

0.002  x  I0~* 

3 

3 

3.28 

0.  004  x  10"6 

(0*'0£=Z)    03S-7WD/N     'XfHJ    1VWM3H1 
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Monoenergetic  Fast  Source.      Diffusion  length  results  based  on 
age -diffusion  theory  assuming  a  monoenergetic  fast  source  at  a  pile 
density  of  1.678  g/cc  are  listed  in  Table  13.      These  results  show  a 
fair  degree  of  consistency,  but  have  a  slight  increasing  trend  not 
noticed   in  the  double  iteration  results  based  on  constants  evaluated 
at  elevation    z  =  20  inches. 


Table  13.     Diffusion  length  results  based  on  a  mono- 
energetic point  source  of  fast  neutrons. 


Elevation  of 
Lowest  Data 
Point 


Number  of  Data 
Points     Used 


Diffusion  Length 


inches 


46 

42 
38 
36 

30 
28 
26 


10 
II 
12 
13 
14 
15 
16 
17 
18 


53.09  ±    0.81 

53.52  i   0.74 

53.50  ±    0.64 

53.45  ±    0.55 

53.79  ±    0.53 

53.76  ±    0.47 

53.64  ±    0.42 

53.90  ±    0.40 

53.92  ±   0.37 


Examination  of  the  results  listed  in  Table  13  indicate  that  the 
assumption  of  a  monoenergetic  fast  source  for  the  determination  of 
diffusion  length  is  reasonably  good  over  the  region  studied.     In  order 
to  more  firmly  establish  the  validity  of  this  assumption,   the  effect  of 
variation  of  the  age  of  Pu-Be  neutrons  should  be   studied.     The  value 
used  in  this  work  was  426  cm2  as  reported  by  Steichen  (26).     Table  14 
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lists  diffusion  length  results  for  several  values  of  age. 


Table  14.     Effect  of  variation  of  age  on  determinations 
of   Li  assuming  a  monoenergetic  fast  source. 


Number  of       : 
Data  Points     : 

Fermi  Age  of        : 
Pu-Be  Neutrons  : 

Diffusion 

^ength 

cm2 

cm 

10 

446 

53.08 

± 

0.81 

10 

436 

53.09 

± 

0.81 

10 

426 

53.  09 

± 

0.81 

10 

416 

53.09 

± 

0.81 

10 

406 

53.10 

± 

0.81 

10 

396 

53.10 

± 

0.81 

10 

386 

53.11 

± 

0.81 

10 

376 

53.11 

± 

0.81 

The  obvious  conclusion  to  be  drawn  from  Table  14  is  that  the 
results  obtained  by  assuming  a  monoenergetic  fast  point  source  are 
almost  completely  insensitive  to  normal  uncertainty  in  the  value  of 
the  Fermi  age. 

Gaussian  Range  Empirical  Source.       The  results  of  the  deter- 
mination of  diffusion  length  based  on  an  empirically  determined  fast 
source  described  by  Gaussian  ranges  and  range  fractions,   at  a  pile 
density  of  1.678  g/cc,   are  given  in  Table  15.     The  uniformity  of  the 
results  is  quite  good,    however,    the  average  is  slightly  lower  than  that 
of  the  point  thermal  source,     double  iteration,   and  monoenergetic  fast 
source  methods. 


Table  15.   Diffusion  length  results  from 
Gaussian  range  empirical  source. 
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Elevation  of   !    Number  of 
lowest  data  point; data  points  used  j 


Diffusion 
length 


Inches 


: 


50 
46 
44 
42 
38 
36 

30 


7 

8 

9 

10 

II 

12 

13 

14 


52.82  ±  1.56 

52.96  ±1.19 

53-33  *  0.92 

53.01  ±0.83 

52.66  ±  0.70 

52.99  ±0.63 

52.85  ±  0.55 

52.76  ±0.47 


Table   16.     Effect  on  diffusion  length  of  using 
varying  Gaussian  range  sets. 


Set     :No.   of:Error  squared  of  Equivalent  Age   :     Diffusion 
number:terms   :   fit  to  slowing     ss  f^v^/h,  cm2  :  length,  cn- 
:  i     down  density       :  : 


cm 


1  3 

2  3 

3  3 

4  3 

5  3 


7.94  x  103 
8.17  x  103 
9.08  x  103 
9.17  x  103 
1.51  x  104 


378.1  53.53  ±1.58 
379.6  53.56  ±1.58 

397.2  53-56  ±1.58 
396.5  53.12  ±1.57 
411.9  52.82  ±  1.56 
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Foulke  (9)  determined  several  sets  of  ranges  from  one  set  of 
indium  foil  data.     Each  set  was  determined  by  an  exponential  stripping 
method.     The  differences  between  the  several  sets  of  ranges  was  in 
the  treatment  of  the  raw  indium  foil  count  rates.     Table  16  lists  dif- 
fusion length  results  obtained  using  several  of  these  sets  with  nine 
data  points.     The  elevation  of  the  lowest  point  was  50  in. 

Set  number  five  was  used  in  the  calculations  for  this  work  be- 
cause it  appeared  to  be  the  best  fit  of  the  experimental  slowing  down 
data  (9).     The  empirical  curve  based  on  this  set  passed  within  the  ex- 
perimental deviations  of  the  slowing  down  data  for  a  greater  portion 
of  the  points  than  did  any  of  the  other  curves.     Details  on  the  data  treat- 
ment and  analysis  may  be  found  in  reference  (9).      This  evidence  is 
presented  here  to  point  out  the  fact  that  the  determination  of  diffusion 
length  in  this  manner  requires  a  very  precise  knowledge  of  the 
Gaussian  ranges  and  range  fractions.     The  small  variations  in  the  age 
and  the  equivalent  age  for  the  two  methods  indicate  that  neither  can 
be  said  to  be  highly  sensitive  to  determinations  of  these  quantities. 

General  Comparison.      The  behavior  of  each  of  the  methods 
studied  is  illustrated  in  Fig.    17.     For  clarity,   confidence  limits  have 
been  eliminated  from  this  plot.     However,  the  confidence  limits  for 
each  value  are  listed  in  the  respective  tables  of  results  for  each  method. 
It  is  apparent  from  this  figure  that  all  of  the  methods,  with  the  excep- 
tion of  the  constant  thermal  source,   give  results  of  essentially  the  same 
magnitude  when  considered  only  in  the  region  where  consistent  results 
were  obtained.       Table  17  lists  the  value  of  diffusion  length  obtained 
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from  each  method.     These  values  were  obtained  by  taking  the  average 
of  all  determinations  which  exhibited  only  an  apparent  random  distri- 
bution about  the  average.     The  confidence  limits  quoted  with  these  values 
are  the  confidence  limits  obtained  using  all  points  in  the  region  of  con- 
sistent results.     Values  are  reported  at  the  KSU  pile  density  and  at  a 
density  of  1.6  g/cc  as  is  commonly  quoted  in  the  literature. 

Using  as  criteria  for  judgment;    (I)  the  degree  of  consistency  of 
results  as  data  are  included  closer  to  the  source,     (2)  the  number  of 
determinations  for  which  this  consistency  is  observed,   and    (3)    the 
extent  to  which  the  average  value,     reported  in  Table  17,   of  a  given 
method  conformed  to  the  averages  of  the  other  methods,   after  reject- 
ing the  ratio  and  constant  thermal  source  methods,   the  best  method 
appeared  to  be  the  double  iteration  procedure,   based  on  one-group 
diffusion. 

It  can  be  seen  from  Fig.    17  that  age-diffusion  results  did  not 
exhibit  the  erratic  behavior,   as  points  were  included  near  the  source, 
that  was  displayed  in  the  results  obtained  using  the  thermal  diffusion 
model.     This  leads  to  the  conclusion  that  age-diffusion  may  be  the  more 
reliable  method  for  determination  of  diffusion  length.     However,  the 
slight  increasing  trend  noticed  in  the  monoenergetic  fast  source  results 
as  data  were  included  closer  to  the  source  cast  some  doubt  on  the 
reliability  of  this  method.     The  Gaussian  range  method  has  been  shown 
to  be  sensitive  to  determination  of  ranges  and  range  fractions,     and 
gives  an  average  value  somewhat  lower  than  would  be  expected  after 
examining  the  results  of  the  other  methods. 

The  point  thermal  source  gave  consistent  results  only  at  a  con- 
siderably greater  distance  from  the  source  than  the  other  methods  and 
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was  ruled  out  for  this  reason.     The  constant  thermal  source  method 
and  the  ratio  method  failed  to  produce  satisfactory  results.     The 
double  iteration  method  gave  highly  consistent  results  to  within  30  in. 
of  the  source,   and  its  average  value  was  within  0.  08  cm    of  the 
average  values  of  the  point  thermal  and  monoenergetic  fast  source 
methods. 

Based  on  the  double  iteration  method,   the  value  of    L    reported 
here  will  be  53.  61  ±  0.  60  cm,   at  a  graphite  density  of  1 .  678  g/cc. 
At  a  density  of  1.  60  g/cc,   this  value  becomes  56.  22  ±  0.  60  cm. 
Table  18  lists  this  value  of    L    along  with  several  other  values  ob- 
tained from  the  literature. 


Table  18.     Comparison  of  measurements  of  the  dif- 
fusion length  of  thermal  neutrons  in  graphite. 


Diffusion  Length 

Graphite  Den: 

iity 

References 

L,      cm 

P, 

gr/cc 

50.0 

1.62 

(10) 

50.2 

1.62 

(II),  (32) 

52.39 

1.60 

(3) 

53.39    -    55.21 

1.60 

(23) 

54.40  ±0.5 

1.60 

(?) 

56.22  ±0.6 

1.60 

this  work 

The  range  of  diffusion  lengths  quoted  from  reference  (23)  include 
several  values  obtained  in  several  different  piles.     The  fact  that  the 
diffusion  length  measured  in  the  KSU  pile  is  higher  than  the  literature 
values  is  attributed  to  the  purity  of  the  graphite  (described  in 
"EXPERIMENTAL  FACILITIES"). 
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Confidence  Limits.      Unless  otherwise  specified,    all  confidence 
limits  placed  on  determinations  of  diffusion  length  in  this  work  were 
calculated  in  the  following  manner:    First,   the  standard  deviation  of 
7n    was  calculated  from  a  linear  regression  analysis  based  on  the 
corrected  count  rates  from  which  that  value  of    "Yn    had  been  deter- 
mined.    The  standard  deviation  of    L    was  then  calculated  by  means 
of  propagation  of  errors.     An  IBM  650  program  written  to  perform 
this  analysis  is  described  in  Appendix  L. 

Material  Buckling  and  Effective 
Multiplication  Factor 

The  material  buckling  of  the  fully  loaded  KSU  pile  was  determined 
experimentally  using  Eqs.   (41)  and  (42-)  of  the  Theory.     Equations  (45) 
and  (47)  were  used  to  evaluate  the  effective  multiplication  factor  and 
the  infinite  multiplication  factor  respectively.      The  calculation  of  the 
material  buckling  was  based  on  a  point  thermal  source,    so  the  data 
used  was  restricted  to  that  which  was  within  the  previously  determined 

region  of  validity  of  this  boundary  condition.     The  value  of    B   *     used 

-4  -2 

in  the  remainder  of  the  calculations  was  7.  755  ±  0.  41  x  10        in. 

This  value  is  the  average  of  five  separate  determinations  and  the  con- 
fidence limits  are  the  standard  deviation  of  the  mean.     The  infinite  multi- 
plication factor,     k     ,     was  calculated  to  be  1.104  ±0.05.     This  result 
was  then  used  to  extrapolate  the  KSU  pile  to  critical  size  by  varying 
B  *  .     This  was  done  by  writing    B*    as  a  function  of  the  extrapolated 
dimension    a    only.     Thus  by  increasing    a,     the  pile  was  expanded  pro- 
portionally.    Table  19    lists  the  extrapolated  values  of  geometric 
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buckling  and    k   ,,.     These  values  of    k   ,f    are  plotted  vs.   pile  size 
in  Fig.    18. 

Comparison  of  the  values  of  material  buckling  calculated  in  the 
KSU  pile  compare  favorably  with  values  calculated  by  Richey  (21) 
of  Hanford  in  an  eight-inch  lattice  with  natural  uranium  as  fuel. 

Effective  Size.      Certain  problems  were  encountered  in  deter- 
mining the  effective  size  of  the  multiplying  lattice  in  the  KSU  pile. 
Whereas  the  top  and  side  boundaries  could  be  defined  by  adding  the 
extrapolation  distance  of  0.71   \      ,     the  bottom  boundary  required 
more  detailed  considerations. 

It  was  decided  to  treat  the  unfueled  pedestal  of  the  pile  as  a 
reflector,     so  that  the  reflector  savings    6    could  be  approximately 
calculated  from 

fS!L      tanh    -i- 
t  L, 

r 

where     T    is  the  reflector  thickness,   and    L      the  reflector  diffusion 

r 

length.     The  above  equation  is  based  on  a  slab  reactor,   and  will 
thus  overestimate  the  reflector  savings.     The  error  thus  induced 
will  not  seriously  affect  the  results  of  this  analysis,  however.     The 
reflector  was  considered  to  be  present  throughout  the  entire  analysis, 
although  in  a  critical  experiment  it  may  be  desirable  not  to  have  the 
reflector  entering  into  the  analysis. 

The  result  of  the  reflector  savings  calculations  was  used  to  define 
an  equivalent  bare  reactor,    so  that  the  extrapolation  distance  0.71  X 
was  also  included  in  defining  the  lower  boundary  of  the  equivalent 
multiplying  medium. 
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Table  .19.     Extrapolation  values  of    k   ,,  for  the  KSU 

pile  (8 -in.    lattice,  p  .      =     1.689  g/cc. 

v  graph  "' 


eff 


Extrapolated  :        Geometric  Buckling        :  k 

Pile  Dimension  a    :  : 

inches : in  " 

69.375  5216  x  10"'  0.58 

86.375  3365  xlO"6  0.72 

101.375  2443  x  10"'  0.81 

116.375  1853  xlO"6  0.87 

131.  375  1454  x  10"'6  0.91 

146.375  1172  xlO-6  0.95 

161-  375  964  x  10"6  0.97 

171.375  855  xlO"6  0.98 

180.895  767  xlO"6  0.99 

Lattice  Diffusion  Length.       In  the  calculation  of    k      from  Eq.   (47) 
the  value  of    L    used  should  be  that  for  the  multiplying  lattice.     This 
value  may  be  found  from  the  diffusion  length  in  the  pure  moderator  by 
the  equation 

L*    =    L2(l   -f), 

where     L      is  the  lattice  diffusion  length  and    f    is  the  thermal  utilization. 

Thermal  utilization  was  calculated  as  shown  in  Appendix  M,   from 
equations  presented  in  ANL  5800  (20)  which  were  based  on  integral 
transport  theory  in  cylindrical  geometry.     These  equations  considered 
the  case  of  a  rod  of  fissionable  material  surrounded  by  an  annular 
cooling  gap  and  an  annular  moderator  region.     The  value  of    f    obtained 
in  this  manner  was  0.  820.     Glasstone  (1  0)  uses  a  value  of  I.  308  for  t\, 
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the  number  of  neutrons  produced  per  fission,    in  natural  uranium.     If 
the  product  of  the  fast  fission  factor    e,     and  the  resonance  escape 
probability    p    is  assumed  to  be  one  which  will  be  high  (10),   the  four 
factor  formula  for    k      yields    k     =1.073.     This  value  is  slightly 
lower  than  the  value  obtained  in  this  paper.     The  reason  for  this  dis- 
crepancy is  believed  to  be  in  the  calculation  of  thermal  utilization. 
Other  methods  of  calculation  of    f    were  tried,  but  proved  less  satis- 
factory than  that  referred  to  above. 

Since  the  thermal  neutron  flux  in  the  uranium  core  of  a  hetero- 
geneous cell  is  greatly  depressed,   the  absorption  of  neutrons  in  the 
fuel  is  less  for  a  heterogeneous  cell  than  for  a  homogeneous  cell. 
Hence  the  value  of    f    for  a  homogeneous  cell,    calculated  as  0.910, 
could  be  regarded  as  an  upper  limit  of  the  actual  value  of  the  thermal 
utilization. 

Fermi  Age.     The  value  of  the  Fermi  age  to  thermal  used  in  the 
calculations  of  material  buckling  and  effective  multiplication    factor 
was  350  cm2.       This  value  is  reported  by  Weinberg  and  Wigner  (32) 
and  Glasstone  and  Edlund  (II). 

Suggestions  for  Further  Work 

The  basic  problem  encountered  in  determination  of  diffusion 
length  is  in  obtaining  an  accurate  mathematical  descroption  of  the  ther- 
mal flux.     Because  of  the  nature  of  a  sigma  pile,   the  flux  attenuation 
in  the    z  direction  will  always  be  very  nearly  exponential.     Determina- 
tion of  the  fundamental  mode  attenuation  coefficient,     7n  ,    should 
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always  lead  to  knowledge  of  the  diffusion  length,     L. 

One  possible  model  for  further  research  in  diffusion  analysis 
is  the  use  of  a  first  collision  density  function  as  a  boundary  con- 
dition for  age -diffusion  theory.     Instead  of  describing  all  neutrons 
as  starting  the  slowing  down  process  at  the  location  of  the  physical 
point  source,   they  should  be  considered  as  starting  from  the  location 
of  their  first  collisions.      This  would  constitute  a  further  refinement 
of  the  age -diffusion  model. 

Another  such  possibility  is  the  use  of  two  group  methods  to 
determine    L.      In  a  moderating  medium  (no  multiplication)  the  fast 
flux  would  be  independent  of  the  thermal  flux  and  could  be  solved  for 
separately.     The  boundary  condition  of  a  point  source,   applied  to  the 
fast  group,  would  be  necessary  to  completely  define  the  fast  flux. 
The  fast  flux  would  then  be  used  as  a  source  term  for  the  thermal 
flux  equation,   depending  upon  the  knowledge  of  an  appropriate  slow- 
ing down  cross  section.     This  method  would  be  subject  to  the  same 
complications  arising  in  most  two-group  calculations,     that  of 
determining  constants  for  the  fast  group. 
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APPENDIX  A 

Solution  to  the  Thermal 
Diffusion  Equation 

The  geometry  for  the  solution  to  the  thermal  diffusion  equation 
is  shown  in  Fig.    19. 


a 

- b    «-<^ 

> 


A 


Fig.    19.      Pile  geometry  for  solution  to  thermal  diffusion  equation. 


The  equation  to  be  solved  is 
V2<j>    -  k24>    =  0 


(A-l) 


The  solution  is  subject  to  the  boundary  conditions  that  the  flux, 
<t>  ,   be  zero  at  all  extrapolated  boundaries.     This  may  be  stated  mathe- 
matically as 
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♦  (±y,  y.  *)  =  o 

+  (*■* -$.,#).  I) 

<|>  (x,  y,  c)  =  0 

where    a,   b,   and  c    are  the  extrapolated    x,   y,   an#   #    dimensions  of 
the  pile. 

Equation  (4)  may  be  written  in  the  form 

*t_  +  i^;£t_.lt,+  .,  (A.2) 

3x2  3y2  3z2 

Assuming  a  solution  of  the  form 

4>  (x,  y,  z)  =  X(x)  Y(y)  Z(z)  (A- 3) 

equation  (A- 2)  becomes 

X"YZ  +  XY"Z  +  XYZ"   -  k  2XYZ  =  0 
Dividing  both  sides  of  this  equation  by    XYZ    gives 

-y!l  v"  *7" 

4+4   +    4    -    k2=0  J  (A-4) 

In  order  to  meet  the  boundary  conditions  let  each  term  of 
equation  (A-4)  be  equal  to  a  constant  as  follows: 

(A- 5a) 
(A- 5b) 
(A-5c) 


X" 
X    = 

-  a1 
m 

Y" 
Y     " 

-   0* 
*n 

Z" 

T  = 

+  c 
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Rewriting  equation  (A-5a)  gives 

X"  +   a  2      X  =  0 
m 

A  solution  of  this  equation  is 

X(x)  =  A      cos  a       x 
m  m 

To  meet  the  boundary  conditions  that    X(a/2)  =  0, 

a       =  ,m=   1,    3,    5,    ••• 

ma' 

so  that  the  solution  for    X(x)    is 

X(x)  =  Am  cos     m^x  ,  m  =   I,    3,    5,    ■•',  (A-6) 

Similarly  for    Y(y), 

Y(y)  =  An  cos  J1^,n  =1,    3,    5,    •  •  •   ,  (A-7) 

Rewriting  equation  (A- 5c)  gives 

Z"  -    7  2         Z  =  0, 
mn 

which  has  as  a  solution 

Z(z)  =  C         sinh  y        z 
mn  mn 

This  function  may  be  made  to  fit  the  boundary  condition  by  altering  the 
argument  to  give 

Z(z)  =  Cmnsinh7mn(c-z).  (A-8) 

Combining  the  solutions  for    X(x),    Y(y)  and  Z(z)  in  the  manner 
indicated  previously  and  combining  coefficients  yields 
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CO  oo 

o>(x,  y,  z)  =    /  /        A         cos     m1TX  cos    "J^     sinh  7  (c-z) 

Li         Li  mn  a  b  mn  '  ' 

m=  1    n=  1  (A-9) 

where    m    and    n    are  summed  over  odd  integers.     Substituting 
(A-5a),    (A-5b)  and  (A-5c)  into  equation  (A-4)  yields 

-a  2       -   02     +    t  2       =  K  a  (A- 10) 

m         K  n  mn  '  ' 

which  shall  be  referred  to  as  the  auxiliary  separation  constants 
equation  in  this  work.     Substituting  in  the  expressions  for    a  and 

fi       the  above  equation  becomes 

.     (  JUJLyi       .     (      »2L)  *    +    y    2  =    K    2 

a  b  mn 

where    a    and    b    are  the  effective  lateral  dimensions  of  the  pile  and 
k     is  the  reciprocal  of  the  diffusion  length. 

Discussion  of  this  solution  can  be  found  in  Glasstone  and 
Edlund  (ll)  and  many  other  Nuclear  Engineering  texts. 

When  considering  a  Multiplying  medium,    equation  (A-l)  becomes 

Vs  <)>  +  B^  a>  =  0  |  (A-ll) 

By   a  completely  analogous  procedure,    the  solution  to  this  equation  is 

oo  oo 

<|>  (  x,  y,  z)  =    /  /         A         cos    — - —    cos  — t-^-  sinh   7         (c-z) 

'  Li        Li  mn  a  b  mn 

m=  1    n  =  1  (A-12) 

where    m    and    n    are  odd  and  the  auxiliary  separation  constants  equation 
becomes 


+  a  2       +    02     -     7  2        =    B2  (A- 13) 

m  ~  n  mn  m  x  ' 
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APPENDIX   B 

Solution  to  the  Fermi  Age  Equation  and 
the  Thermal  Diffusion  Equation 

The  geometry  for  this  solution  is  shown  in  Fig.  20  . 


Fig.    20.     Pile  geometry  for  solution  to  the  age-diffusion  problem. 
The  equations  involved  were  the  Fermi  age  equation 

and  the  thermal  diffusion  equation 

DVz<t>   -  S    4>  +  S  =  0       . 


(B-l) 


(B-2) 


The  Fermi  age  equation  was  subjected  to  the  boundary  conditions 
that    q=0at    x  =  ±a/2,   y=±b/2,     and    q  -*0  as  |  z[  *•   oo.     The  source 
condition  was  that  all  source  neutrons  entered  the  system  at 
x  =  y  =  z  =  0 ,   atanage   T  =  0 . 
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This  condition  was  written  inrto  the  Fermi  age  equation  in  the 
form 

V2q=-^3-     .     Qq    6    (T)    6i(x)    6{y)    6(z).  (B-3> 

The  boundary  conditions  to  the  thermal  diffusion  equation  were 
that  <t>  =0  at  x  =±a  /2,  y  =±b/2,  and  o>  — 0  as  |  z|  —  oo  .  The  source 
condition  was  that  S  =  q(x,   y,    z,   r  J, 

Q    ,    q    ajid    <|>     were  expressed  as  double  Fourier  series  satis- 
fying the  boundary  conditions 
oo  oo 

*=2         I       ♦««««•-    ^  (B-4) 

m  =1     n  =  l 


00  00 

* =  Z    Z    <w  (z)  cos  i2r^  cos  -nF  (B-5) 

m  =  l   n  =  l 

and 

oo  oo 

Qo«T)6(x)6(y)6(z)  =    J,        Z       Qmn(z) 


a  b 

m  =  l      n  =  l  (B-6) 

where    m    and    n    are  summed  over  odd  integers.     Each  side  of  Eq.   B-6 
was  multiplied  by  orthogonal  functions  and  then  was  integrated  over  the 
range  of  orthogonality  which  resulted  in 

oo         oo        .   _ 
Qo6(T)6(x)6(y)6(z)    n  £         £    —^  6(  7  )  6(  z)  cos    -E-^cos-^ 
m  =  l     n  =  l  i 

(B-7) 
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Substitution  of  Eqs.   B-5  and  B-7  into  the  age  equation,    Eq.   B-3 


gave 


GO  00 

I  I 

m  =  l     n  =  lL 


4Q 


Tt 


2      Hmn'    ab 


B6(T)6(zj=0 
(B-8) 


Taking  the  Fourier  trnasform  of  the  equation  for  a  particular 
m    and    n    and  letting 


(B-9) 


and 


led  to 


b 


4Q 


6  A".    +  (a2      +   02     +  u,2)  q  (co)-— j-°   6(t), 
where    q    (uj     represents  the  Fourier  transform  of    q         (z). 


(B-10) 


(B-ll) 


The  solution  to  this  differential  equation  was 
4Qo        _(„^+w*)T 
M  =  ^b—  e 


(B-12) 


Taking  the  inverse  transform  of  Eq.   B-12  and  substituting  it 
into  Eq.   B-5  gave  the  slowing  down  density  along  the    z    axis  as 


» 


00  DO 

■1   I 

m  =  l    n  =  l 


■(a2    +/32  )t 
v     m      'n' 


4Q     e 
o 


2/4t 


ab  "V  4  tt  T 


(B-13) 
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Substitution  of  Eqs.   B-4  and  B-5  into  the  thermal  diffusion 
equation.     Equation  A-2,   and  designating    2    /D  as  k  2  gave 


oo  oo 

1     I 

m  =  l      n  =  l 


a2o> 


mn 
T — 


dz 


m        'n  mn  D 


=  0   . 


(B-14) 


Taking  the  Fourier  transform  of  Eq.   B-14  for  a  particular    m 

and    n    and  designating  the  Fourier  transforms  of   +  and    q  by 

mn  mn 

<t>(u)  and    q(u),    respectively,    gave 


♦■(«)  = 


q  (<") 


1    m       Mn 


(B-15) 


Taking  the  inverse  transform  of  Eq.   B-15  by  means  of  the 
convolution  theorem  (30),    after  Eq.    B-12  had  been  used  to  express 
q(u  )  gave 


m 

2Q   e 
o 

+ 

/32)t 
n ' 

in 

D  ,     y 
ab      n 

n/4ttT 

f 


-z/4T     -|z-z|y 

I  I        '„ 


dz^ 


(B-16) 


where 


■y        =  n/<*'    +fi'  +k' 
mn  m     'n 


(B-17) 


Glasstone  and  Edlund  (ll),   on  pages  185  to  187,   give  the  develop- 
ment showing  how  the  integral  in  Eq.   B-16  may  be  expressed  in  terms 
of  the  error  function 


rf(t) 


<4      IT 


pt  -U2 

L  e 


du 


(B-18) 
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Writing  Eq.    B-16  in  terms  of  the  error  function  and  substituting 
it  into  Eq.   B-4  gave  the  thermal  flux  along  the    z    axis  as 

'mn  " 

ll  +  erf 


00  00 


♦<->-z  2  tA 


k't 


m  =  1   n  =  1 


ab    'mn 


_? —  v      *rr)] 


Y        z 
mn 


+  e 


f  1  -  erf ( 


+  y         sTt 

ZsIT        mn 


(B-I9) 


To  describe  the  thermal  flux  due  to  an  empirical  »    source  which 

r. 

is  defined  by  a  set  of  Gaussian  ranges,    the  quantity     -4 —      should  be 
substituted  into  Eq.    B- 19  for  the  age    T  .     The  flux  must  then  be 
multiplied  by  the  appropriate    f.     and    summed  over  the  number  of  terms 
in  the  empirical  source,    which  gave 


N       00         00 


♦<•>-!  I  I    Tdk 


i  =  1   m  =  l  n  =  1 


[,,       ,<  z  mn    i\1 

l+erf(_:  .  _ __ )j 


+  e 


r  y        r-  1 

L1  - erf  (—  *  —2 — }J 


(B-20) 


The  above  derivation  is  identical  to  that  presented  by 
L.   R.   Foulke  (9). 
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APPENDIX    C 

IBM  650  Program  for  Diffusion  Length, 
Point  Thermal  Source 

This  code  was  written  to  determine  an  experimental  value    of 
diffusion  length  based  on  the  assumption  of  a  point  source  of  thermal 
neutrons.     The  program  was  written  in  SOAP  II  and  floating  point 
form.     The  object  program  and  a  logic  diagram  are  given  in  this  section. 

The  method  of  calculation  was  an  iteration  procedure  based  upon 
the  equation 

00  00 

,   /  \     \         \         2S  mirx  niry         sinh     mn(c-z) 

9  I  x,  y,  zl  =  /         /  .  r.     cos  cos  — s-* —    p— ■ — *— ^ — '     ■ 

T  '       "    '     Li       Li       abD  a  b  cosh  y       C 

m  =  1  n  =  1 

Values  of  the  correction  factors    C„    and    C,, ,   were  calculated  using 
an  initial  estimate  of    7n  .   the  data  was  corrected  according  to  these 
factors,   and  a  new  value  of    7u  was  obtained  by  a  least  squares 
analysis  of  the  corrected  data.     The  new  value  of    fu  was  then  used 
to  calculate  a  second  set  of  correction  factors.     This  procedure  was 
repeated  until  the  least  squares  value  of    7u   obtained  at  the  end  of 
one  trial  was  within  a  specified  precision  of  the  value  obtained  from 
the  previous  trial.     At  this  point  the  last  value  of    7n  obtained  was 
used  to  calculate  the  diffusion  length  according  to 


where 


IT        0         ir 


b 
a    and    b    being  the  extrapolated  pile  dimensions. 


Input  for  this  code  consisted  of  the  data  and  accompanying 
position  co-ordinates  as  well  as  certain  parameters  listed  in  Table  20. 
Table  21  gives  the  form  of  the  two  types  of  output  of  this  program. 


Table  20.     Input  data  for  IBM  650  program  for  diffusion 
length  based  on  a  point  thermal  source. 


Symbol     :  Explanation  :  Drum  Storage 

:  :       Location 

FIRST  Initial  value  of    Til  0100 

POINT  No.   of  data  points,   in  form  00  0000  OOxx  0101 

DATPT         No.   of  data  points,   in  floating 

point  form  0102 

A  Extrapolated  x-dimension  0103 

B  Extrapolated  y- dimension  0104 

C  Extrapolated  z-dimension  0105 

X  x  co-ordinate  of  data  0106 

Y  y  co-ordinate  of  data  0107 

Z  z  co-ordinate  of  data,    to  be  stored 

consecutively  starting  at  0301 

N  neutron  count  rates,    to  be  stored 

consecutively  starting  at  0501 
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The  input  parameters  and  data  listed  in  Table  19  should  be 
punched  on  one-word  load  cards  and  fed  into  the  machine  with  the 
object  program.  Point  Stripping  may  be  accomplished  by  listing 
the  data,  N  and  z,  with  the  points  to  be  dropped  last  and  varying 
the  parameters  POINT  and  DATPT. 

The  form  one  output  is  punched  after  the  calculation  of  each 
corrected  data  point.     One  card  of  this  form  will  be  obtained  for 
each  data  point  during  a  given  trial.     When  the  desired  precision  is 
reached,    the  final  calculations  are  made  and  one  card  of  form  two 
is  punched. 
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Table  21.     Output  forms  for  IBM  650  program  for 

diffusion  length  based  on  a  point  thermal 
source. 


Word  No.    1 


Form  One: 

z  co-ordinate       data 


Form  Two: 

•y,,,     Last 
value 


corrected  data     Harmonic 
Correction 


7n  Next 
to  Last 
value 


POINT 


Yn,     Initial 
value 


End 
Correction 


Diffusion 
Length 


The  operating  time  necessary  to  calculate  the  correction  factor 
for  one  data  point  using  four  harmonics  is  approximately  45  seconds. 
The  capacity  of  the  program  is  200  data  points. 
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APPENDIX   D 


IBM  650  Program  for  Diffusion  Length 
Constant  Thermal  Source 


The  purpose  of  this  program  was  to  determine  an  experimental 
value  of  the  diffusion  length  based  on  the  assumption  of  a  constant 
source  of  thermal  neutrons.     The  program  was  written  in  SOAP  II 
and  floating  point  form.      The  object  program  and  a  logic  diagram  are 
given  in  this  section. 

The  calculations  in  this  program  are  based  on  the  equation 

m+  n+  2 
2 
l6<t>o(-l) 


y  C-  16<M-1J  m1TX 

<)>(x,  y,  z)=    /  /         cos  m1TX   cob— r-^  sinh?       (c-z) 

Cj        Lf       __  _  *  8. j_ v   „,  „  a  b  mn'        ' 


.     mni    sinh  y  c 

1    n  *  1  mn 


The  same  iterative  procedure  is  used  in  this  program  as  was  used 
in  the  point  thermal  source  program  described  in  Appendix  C. 

Input  parameters  and  data  are  identical  to  those  used  in  the 
point  thermal  source  program  and  are  listed  in  Table  20.     Output  is 
similar,   and  is  described  in  Table  22. 


Table  22.     Output  forms  for  IBM  650  program  for 
diffusion  length  based  on  a  constant 
thermal  source. 


Word  No.    1   - 
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Form  One: 

z  co-ordinate      data 


Form  Two: 

Tn.   Last 
value 


7n  Next 
to  Last 
value 


corrected 
data 


■yu  Initial 
value 


Harmonic 
Correction 


POINT 


End 
Correction 


Diffusion 
Length 


The  output  sequence  is  identical  to  that  for  the  point  thermal 
source  program.     The  operating  time  necessary  to  calculate  the 
correction  factor  for  one  data  point  using  four  harmonics  is  about 
45  seconds.     The  capacity  of  the  program  is  200  data  points. 
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APPENDIX   E 

IBM  650  Programs  for  Diffusion  Length, 
Ratio  Method 

The  first  of  these  programs  was  written  to  avaluate  the  double 
integral  function 


pa/ 2      pb/2 

Fli(  z)  =   \  \  <|>  (x,  y,  z)  cos  cos    .  '      dx  dy 

J-a/2     J-b/2  a  b 


The  program  was  written  in  SOAP  II  and  floating  point  form.     The 
object  program  and  a  logic  diagram  are  listed  in  this  section. 

The  program  was  written  to  accommodate  as  many  data  points 
as  were  desired.     The  data  was  taken  in  the  manner  described  in 
the  section  entitled  "EXPERIMENTAL  PROCEDURE".     An  odd 
number  of  data  points  is  required  in  each  of  the  horizontal  traverses 
making  up  one  set  of  data.     This  is  to  accommodate  Simpson's  rule, 
which  is  used  for  the  major  portion  of  the  integration.     The 
trapezoidal  rule  is  used  to  evaluate  the  integral  over  the  area  between 
the  fifth  horizontal  traverse  and  the  extrapolated  boundary. 

Input  constants  and  data  input  forms  are  listed  in  Table  23. 
Flux  data  should  be  stored  in  descending  order  of  traverse,    the 
outermost  traverse,    the  extrapolated  boundary,   being  first.     Each 
traverse  should  be  liste  d  with  its  outermost  point  first.     One  card 
of  form  one  is  needed  for  each  set  of  data.     The  output  form  is 
listed  in  Table  24. 

The  approximate  time  for  evaluation  of  the  integral  using  54 
data  points  is  three  minutes. 

The  second  program  was  written  to  determine  the  diffusion 
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Table  23.     Input  forms  for  double  integral 
evaluation  program. 


118 


Symbol 


Explanation 


Drum  Storage 
Location 


A  Extrapolated  x-dimension 

B  Extrapolated  z- dimension 

BPRIM        Distance  from  center  of  pile  to  outer- 
most fuel  port 


0000 
0001 


0002 


Word 


Form  1        Number  of  Data  Points 
(IR  Form) 

Form  2        Count  Rate 


z  Elevation 
of  Data 

x  co-ordinate       y  co-ordinate 


Table  24.     Output  forms  for  double  integral 
evaluation  program. 


Word 


Integral  Elevation  z 


Number  of  Data  Points 


length  by  the  ratio  method.     It  was  written  in  SOAP  II  and  floating 
point  form.     The  object  program  and  a  logic  diagram  are  included  in 
this  section. 

An  initial  trial  value  of  Yn   is  determined  by 

1  Fl 


Input  constants  and  data  input  forms  are  listed  in  Table  25. 
In  the  notation  used,     Z|    should  always  be  greater  than    z%  . 

Table  25.     Input  forms  for  ratio  method  program. 
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Symbol 


Explanation 


Drum  Storage 
Location 


A 
B 
C 


Extrapolated  x-dimension 
Extrapolated  y-dimension 
Extrapolated    z-dimension 


0100 
0101 
0102 


Word         1 
Form  1     z. 


2 
Z2 


3 


4 
F2 


The  two  output  forms  for  the  ratio  method  program  are  listed  in 
Table  26.  One  card  of  form  one  is  punched  for  each  trial.  Operating 
time  is  between  20  and  30  seconds  for  one  calculation  of  diffusion  length. 

Table  26.     Output  forms  for  ratio  method  program. 
Word    1  2        :  3  :  4  \         5  6 


Form  1: 

Gamma 

Form  2: 
zl 


Error 


Gamma 


Diffusion 
Length 


Ratio 


Number  of 
Trials 
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016  1 

Ul  1  5 

006  9 

0122 

006  9 

017  3 

002  6 

007  6 

02  2  9 

0055 

00  6  9 

0119 

017  2 

0  09  9 

00  6  9 

U169 

0  2  2  2 

014  9 

006  9 

0219 

027  2 

0199 

00  6  9 

026  9 

032  S 

024  9 

006  9 

0319 

037  3 

0299 

030  4 

0057 

020  4 

•  254 

030  4 

0107 

0  2  0  4 

U  304 

020  4 

0157 

016  1 

0165 

021  ti 

0  3  6  9 

0  3  7  2 

0  32  7 

020  4 

0354 

020  4 

036  9 

020  4 

020  3 

00  0  0 

0051 

y  y  9  b 

6H  50 

2  5  7  5 

03  4  9 

913  7 

1  24  R 

156  2 

004  7 

30  3  0 

0045 

060  0 

004  4 

195  1 

0  301 

195  1 

0404 

020  7 

0  160 

195  3 

010  5 

010  6 

0211 

195  3 

0  30  6 

0109 

00  6  2 

195  4 

1)257 

0210 

•  06  3 

04  0  0 

02  53 

02  5  6 

0159 

000  0 

0215 

0108 

0113 

0210 

0260 

016  4 

0017 

002  0 

05S0 

022  4 

0277 

010  9 

016  3 

020  7 

0133 

01  B  6 

014  1 

022  4 

027  9 

018  8 

02  38 

034  2 

0045 

000  1 

0351 

020  7 

0261 

010  9 

•  4  85 

02  4  3 

•  292 

0050 

0650 

800  3 

•  307 

031  0 

•  600 

016  6 

•  191 

010  2 

•  357 

0207 

0163 

034  2 

•  34  2 

0150 

•  700 

800  3 

•  407 

036  0 

0600 

026  5 

•  366 

0000 

•  1  55 

036  5 

•  341 

026  5 

•  318 

010  3 

•  457 

0109 

•  535 

024  3 

•  392 

01  50 

•  750 

800  3 

•  SOT 

0410 

•  600 

031  5 

•  36S 

0000 

•  205 

0315 

•  391 

031  5 

•  4  18 

0188 

•  093 

026  5  0365 

0315  9415 

0070  »323 

0164  *419 
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FBB 

FUNCT 

CHE  CK 

STU 

ER  ROR 

ERR  OR 

LOO 

E  OOCL 

LDO 

GAMMA 

BTO 

197? 

LOO 

ERROR 

8TO 

197B 

PUNCH 

f  CM 

1977 

TR|  AL 

Rill 

ERROR 

BU  i 

CQNT1 

CON  TU 

C  0  NT  i 

fMU 
FAD 

GAMMA 
DELTA 

STU 

GAMMA 

LOOP1 

C  ONT  3 

RIU 

r  aa 

STU 

«  »  u 
FO  V 
STU 

GAMMA 
DELTA 
GAMMA 
OE  L  TA 
TEN 
0  E  L  T  A 

F8R 

PH  E  CN 

CHE  CK 

BM  1 

CUNT3 

LOOP  1 

PRE  C  1  B 

C  OH   I  3 

RAU 
FOX 
STU 

P  1 

A 

ALPH3 

FMP 

A  LPH3 

ALPHA 

STU 

ALPH3 

SQUARE 

RAU 

P  1 

FO  V 

B 

a  tu 

BETAS 

FMP 

bE  T  A3 

BETA 

8  TU 

BETAS 

8  0  U  A  H  E 

RAU 

GAMMA 

GAMMA 

FMP 

G  AUM  A 

SQUARE 

FSB 

A  L  PM3 

FSB 

U  E  T  A3 

LOO 

EOO  A  U 

8TU 

KAPPA 

KAPPA 

RAU 

ONE 

FO  V 

KAPPA 

0  1  F  FUB 

FMP 

CON  V  T 

STU 

OL 

LENGTH 

LDO 

20  N  E 

STO 

1977 

LOO 

ZT  HO 

BTO 

1978 

LOO 

GAMMA 

STO 

1979 

LOO 

OL 

B  TO 

19B0 

LDO 

RATIO 

BTO 

196  1 

LDO 

BOOS 

BTO 

1983 

PCH 

1977 

SOOO 

ay& 

0419 

33 

0070 

S047 

307 

0047 

31 

000  3 

•  355 

308 

0255 

69 

0158 

•  4S0 

309 

0150 

69 

024  3 

8095 

310 

0095 

24 

197  7 

•  1  30 

311 

0130 

69 

000  a 

•  305 

212 

0305 

34 

197  8 

•  381 

313 

0281 

71 

197  7 

■  337 

314 

0337 

60 

000  3 

•  557 

315 

0557 

46 

046  0 

•  311 

316 

0460 

60 

034  3 

•  097 

317 

0097 

33 

035  6 

•  333 

318 

0233 

31 

034  3 

•  045 

319 

0  311 

60 

024  3 

•  147 

230 

0147 

33 

0256 

0363 

331 

02S3 

31 

024  2 

Ul  45 

233 

0145 

60 

025  6 

036  1 

223 

0361 

34 

030  0 

080  0 

334 

0000 

31 

025  6 

0309 

335 

0209 

33 

035  0 

0377 

236 

0377 

46 

01B0 

0045 

327 

01B0 

60 

0300 

0355 

238 

0355 

34 

010  0 

9850 

339 

0650 

31 

04  54 

•  607 

230 

0607 

39 

04  5  4 

U504 

331 

0504 

31 

045  4 

0657 

333 

0657 

60 

0300 

0405 

233 

0405 

34 

010  1 

•  40  1 

234 

0401 

31 

030  6 

0359 

335 

0359 

39 

0306 

0  356 

336 

03S6 

31 

030  6 

•  309 

337 

0309 

60 

024  3 

0197 

230 

0197 

39 

024  3 

04  43 

239 

0443 

33 

045  4 

0331 

24  0 

0331 

33 

030  6 

0333 

34  1 

0333 

69 

0086 

0500 

343 

0086 

31 

0090 

•  143 

24  3 

0143 

60 

0000 

•  455 

344 

0455 

34 

009  0 

0140 

345 

0140 

39 

0250 

0900 

246 

0900 

31 

055  4 

0707 

247 

0707 

69 

030  7 

•  510 

24  6 

0510 

34 

1977 

•  230 

249 

0330 

69 

0109 

•  118 

350 

0113 

34 

197  8 

•  381 

351 

03B1 

69 

024  8 

0195 

352 

0195 

34 

197  9 

0333 

353 

0333 

69 

0554 

•  757 

254 

0757 

34 

1960 

•  383 

255 

0383 

69 

016  4 

•  06? 

256 

0067 

34 

198  1 

•  134 

357 

0134 

69 

0005 

•  341 

358 

0341 

34 

1983 

•  585 

359 

058S 

71 

197  7 

8000 
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APPENDIX   F 


IBM  650  Program  for  Diffusion  Length 
Double  Iteration  Method 


This  program  was  written  to  evaluate  the  diffusion  length  by  the 
double  iteration  technique.     It  was  written  in  SOAP  II  and  floating  point 
form.     The  object  program  and  a  logic  diagram  are  included  in  this 
section. 

A  detailed  description  of  the  double  iteration  process  can  be 

found  in  the  THEORY.     The  constants    A  are  determined  by  the 

mn  ' 

equation 


4  F 


A 


mn 


mn  ab  sinh  y         (c  -  z    ) 

mn  v  o 

where    z      is  the  elevation  at  which  the  functions    F       ,   determined  by 
o  mn 

a/2         b/2 

Fmn  =    la/2     lb/2  *  (X'Y'  ^  C°S  JI^L    C°S  ^  ^  ^ 

are  evaluated.     The  first  portion  of  this  program,  which  calculates 

the    F        's    is  an  adaptation  of  the  double  integral  evaluation  program 

described  in  Appendix  E.     The  remainder  of  the  program  performs 

the  double  iteration  analysis.     Several  determinations  may  be  desired 

with  one  set  of    F        "s.     If  this  is  the  case,   the  data  for  each  additional 
mn 

determination  should  be  followed  by  a  transfer  card  which  sends  the 
machine  to  the  CON11  location.  This  will  avoid  unnecessary  recal- 
culation of  the    F        's.     The  correction  factors  used  are  listed  in 

mn 

Table  1. 

Input  to  this  program  is  partially  in  the  form  of  one  word  load 
cards.     Certain  parameters,  plus  the  data  for  the  vertical  traverse 
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are  of  this  form.     This  input,   Hated  in  Table  27,   should  be  followed 
by  a  transfer  card  and  the  data  for  evaluation  of  the  double  integral 
functions.     The  data  for  the  evaluation  of  these  functions  is  exactly 
the  same  as  the  input  to  the  double  integral  program  described  in 
Appendix  E,  forms  one  and  two. 


Table  27.     One  word  load  input  for  double 
iteration  program  . 


Symbol 


Explanation 


Drum  Storage 
Location 


A 
B 
BPRIM 

HARM 

FIRST 
POINT 

DATPT 

C 

N 


Extrapolated    x-dimension  0000 

Extrapolated    y-dimension  0001 

Distance  from  center  of  pile 

to  outermost  fuel  port  0002 

Number  of  highest  harmonic, 

floating  point  form  0003 

Initial  Estimate  of  Gamma  0004 

Number  of  Data  Points  (vert,   traverse) 

in  form    00    0000    OOxx  0005 

Number  of  Data  Points  (vert,   traverse) 

in  floating  point  form  0006 

Extrapolated    z-dimension  0007 

Count  Rates  for  vertical  traverse, 

stored  consecutively  starting  at  0650 


Co-ordinates  of  vertical  traverse 
data,   stored  consecutively  starting  at 


0675 


The  capacity  of  the  program  is  25  data  points  in  the  vertical 
traverse  and  100  data  points  in  the  horizontal  mesh.  The  program 
capable  of  using  up  to  the  ninth  harmonic. 
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There  are  three  output  forms,   listed  in  Table  28.     The  first  lists 

the  constants  and  may  use  two  cards.     The    A       's    are  listed  in  the 

'  mn 

order    An,   2A13,   2AI5,    ...  A33,   2A35,   ...  A55,   etc..      Where    m    and 
n    are  not  equal,  the  constant  is  multiplied  by  two  to  account  for  both 
the    mn    and    nm    terms  at  once.     This  can  be  done  because  of  the 
symmetry  of  the  pile.     Form  two  is  the  standard  correction  factor 
output,  and  form  three  is  a  final  output  form.     Form  three  is  punched 
after  the  completion  of  each  correction  factor  iteration. 

Table  28.     Output  forms  for  double  iteration  program. 

Word      1:2:3:4:5  :         6  ':  7:8 

Form  I:  (3  harmonics) 

A„  2A13  2Al5        A33         2A35  A55 

Form  2: 

Z  N  N  C„         C_  -  -  - 

„  .  corr        H  E 

Form  3: 

FIRST       Next  to      Last       Data        Diffusion       -  - 

last  7  7  Points        Length 

Five  to  ten  minutes  are  required  for  complete  double  iteration 

analysis  using  10  data  points.     An  additional  three  minutes  is  required 

for  evaluation  of  each    F  used, 

mn 
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OBJECT  PROGRAM  -  APPENDIX  F 
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BLR 

1951 

1960 

BLR 

1977 

19B  4 

Bl  R 

1900 

1999 

BLR 

0  100 

0790 

SfN 

FLUX 

0100 

SYN 

X 

0300 

8YN 

0300 

8VN 

F  UN  CT 

0400 

SIN 

illiPF 

OS  00 

8  TN 

TR  *  Pf 

06  00 

8  VN 

Z 

0650 

8YH 

0675 

8V  N 

N  COR 

07  00 

8V  N 

G  AHUM 

0730 

8VN 

*  UN 

0750 

8  V  N 

F  MN 

0770 

SYN 

* 

00  00 

SV  N 

B 

0001 

SYN 

BPR  1  M 

0003 

S  V  N 

HARM 

0003 

BY  N 

FIRST 

0004 

3YN 

POINT 

0005 

SYN 

0  *  TP  T 

0006 

8  YN 

C 

00  07 

SYN 

START 

1000 

ZERO 

00 

0000 

00  00 

ONE 

00 

0000 

0001 

TWO 

00 

0000 

oooa 

THREE 

00 

0000 

0003 

FIVE 

00 

0000 

00  05 

8  1  X 

00 

0000 

0006 

EIGHT 

00 

0000 

00  08 

TEN 

00 

0000 

00  10 

SIXTY 

00 

0000 

0060 

FP1 

10 

0000 

0051 

tpz 

20 

0000 

00  51 

pi 

31 

4  159 

0051 

MM 

4  0 

0  0  00 

0051 

FP6 

6  0 

o  o  no 

00  51 

FP8 

8  0 

0000 

00  5  1 

FP9 

9  0 

00  00 

005  1 

FP  16 

16 

00  00 

00  5i! 

C  ON  VT 

25 

3  600 

00  5  1 

1O00O 

10 

00  00 

00  55 

1  NOX  * 

00 

2  400 

0000 

1  NO* 

00 

2600 

00  00 

1  NOR 

00 

4600 

0000 

EOOE  * 

3  TO 

A  A  A  1 

R  A  U   A  A  A  2 
A  A  3 


FSB 

BH  I 
S  TU 

R  A  U 
FMP 

STU 

RAU 
LOD 

FMP 
S  TU 
R  A  U 
BM  I 


A  A  A  28 

A  A  A  2 

A  A  A  4 

A  A  A  9 

A  A  A  4 

A  A  A  2 

A  A  A  10 

A  A  A  11 

A  A  A  2 

A  A  A  4 

A  A  A  12 

A  A  A  4 

A  A  A  2 


A  A  A  4 

A  A  A  13 

A  A  A  14 

A  A  A  15 

A  A  A  13 

A  A  A  3 

A  A  A  13 

B003 

A  A  A  1  B 

A  A  A  3 

A  A  A  a 

A  A  A  19 

A  A  A  37 

A  A  A  20 

STU  A  A  A  31 

R AU   A  A  A  a 

A  A  A  a 

A  A  A  33 

A  A  A  31 

A  A  A  24 

A  A  A  19 

A  A  A  19 

A  A  A  34 

A  A  A  19 

A  A  A  25 


R  AU 
FO  V 
RAL 
8  TO 
R  AU 
FAD 
STU 
R  AU 
STU 


FMP 
9  T  kt 
FO  V 
STU 
F  AO 
STU 


F  I)  V 
F  9  | 
BM  I 
R  AU 


A  A  A  19 


1 

0000 

00 

0000 

0000 

2 

0000 

00 

000  0 

0000 

3 

0000 

00 

0000 

0000 

0000 

00 

000  0 

0000 

5 

0000 

00 

000  0 

0000 

6 

0000 

00 

0000 

0000 

7 

0000 

00 

0000 

0000 

8 

0000 

00 

0000 

0000 

9 

0000 

00 

000  0 

0000 

10 

0000 

00 

0000 

0000 

1  1 

0000 

00 

000  0 

000  0 

13 

0000 

00 

00  0  0 

0000 

1  3 

0000 

00 

0000 

0000 

14 

0000 

00 

0000 

0000 

15 

0000 

00 

0000 

0000 

16 

0000 

00 

0000 

000  0 

17 

0000 

00 

0000 

0000 

18 

0000 

00 

0000 

0000 

19 

0000 

00 

0000 

0000 

20 

0000 

00 

OOO  0 

0000 

31 

0000 

00 

0000 

0000 

33 

0000 

00 

0000 

000  0 

2  3 

0000 

00 

000  0 

00  0  0 

24 

0000 

00 

00  0  0 

00  0  0 

25 

0000 

00 

000  0 

0000 

26 

0050 

00 

000  0 

0000 

27 

oeoo 

00 

000  0 

000  1 

28 

0850 

00 

000  0 

OOO  2 

39 

0900 

00 

000  0 

000  3 

30 

0950 

00 

0000 

0005 

31 

1050 

00 

oooo 

0006 

32 

1100 

00 

00  00 

0008 

33 

1150 

00 

00  00 

0010 

34 

1200 

00 

0000 

0060 

35 

1250 

10 

00  0  0 

00  51 

36 

1  300 

20 

0000 

005  1 

37 

1350 

31 

415  9 

0051 

38 

1400 

40 

0000 

0051 

39 

1450 

60 

000  0 

005  1 

40 

1500 

80 

000  0 

005  1 

4  1 

1550 

90 

00  0  0 

0051 

42 

1600 

16 

00  0  0 

0053 

4  3 

1650 

35 

360  0 

0051 

4  4 

1700 

10 

0000 

005  5 

4  5 

1750 

00 

24  00 

000  0 

4  6 

1800 

00 

2600 

0000 

47 

1850 

00 

460  0 

0000 

EXIT  INSTR 

48 

0051 

34 

00  5  4 

0057 

STORE   ARC 

4  9 

0057 

30 

001  1 

0014 

ZERO 

SO 

0014 

60 

0017 

0021 

51 

0021 

37 

00  11 

0037 

52 

0037 

21 

00  4  2 

004  5 

ONE 

5  3 

004  5 

60 

00  4  8 

0053 

COE  FF 

54 

00  53 

31 

00  0  8 

0061 

55 

0061 

60 

004  2 

004  7 

FIVE 

se 

0047 

33 

080  1 

00  2  7 

57 

00  27 

46 

0030 

0031 

5B 

00  31 

31 

00  4  2 

009  5 

59 

0095 

60 

000  8 

0013 

E  TP  i 

60 

0013 

39 

00  16 

0066 

61 

0066 

31 

OOO  8 

006  1 

63 

0030 

60 

004  2 

0097 

ONI 

63 

0097 

33 

004  8 

003  5 

64 

0025 

46 

002  8 

002  9 

65 

0029 

31 

004  2 

0795 

66 

0795 

60 

000  8 

006  3 

E   T  0   1 

67 

006  3 

39 

0816 

0866 

60 

0866 

31 

000  8 

00  30 

69 

0028 

60 

004  2 

0797 

POINT   TWO 

70 

0797 

33 

085  1 

0077 

71 

0077 

46 

0080 

OOB  1 

72 

008  1 

31 

00  4  2 

08  4  5 

7  3 

0845 

60 

000  8 

08  13 

E   TO  P  T   2 

74 

0813 

39 

0916 

0966 

75 

0966 

21 

000  8 

002  B 

76 

0080 

60 

004  2 

OB  47 

77 

0B47 

69 

090  1 

0104 

78 

0901 

39 

000  8 

0058 

79 

0058 

21 

0012 

00  15 

BD 

0015 

60 

0011 

0065 

81 

0065 

46 

0018 

0019 

62 

0019 

65 

0013 

00  5  4 

83 

0018 

60 

004  8 

0803 

84 

0803 

34 

0012 

0062 

85 

0062 

65 

800  3 

0054 

EXIT   1  N8TR 

86 

080  4 

S4 

080  7 

0010 

ONE 

87 

00  10 

60 

004  8 

08  53 

80 

085  3 

33 

00  4  2 

0069 

89 

0069 

21 

002  4 

0B27 

TWO 

90 

0837 

60 

083  0 

0035 

91 

00  35 

31 

004  0 

004  3 

OENDH 

92 

004  3 

31 

009  8 

0951 

93 

0951 

60 

004  2 

0B97 

94 

0897 

39 

004  2 

0093 

NUMERATOR 

95 

009  3 

31 

00  4  6 

0049 

96 

0049 

34 

00  9  8 

0798 

NTH   TERM 

97 

0  7  98 

21 

00  5  2 

005  5 

99 

0055 

33 

00  2  4 

1001 

TOTIL 

99 

1001 

Bl 

00  2  4 

0877 

100 

0677 

60 

00  5  2 

08  57 

101 

0857 

34 

002  4 

0074 

CRITERIA 

102 

0074 

33 

092  7 

090  3 

103 

090  3 

46 

00  S  6 

0907 

104 

00  56 

60 

002  4 

0807 

133 


*  *  AS  6 

R  AU 

A  A  A  ao 

IDS 

0907 

60 

004  0 

0895 
007  5 

FAD 

A  A  A  3 

105 

0895 

38 

00  4  8 

STU 

A  A  A  ao 

NEI   N 

107 

0075 

ai 

004  0 

009  3 

FM  P 

A  A  a  ai 

108 

0093 

39 

009  8 

08  4  8 

STU 

A  A  A  21 

NED  DENOM 

109 

0848 

81 

009  a 

10  5  1 

R  A  U 

A  A  A  33 

110 

1051 

60 

00  4  6 

11  0  1 

FUP 

A  A  A  3 

111 

1101 

39 

004  3 

0792 

STU 

A  A  A  83 

A  A  A  83 

118 

0792 

81 

004  6 

00  4  9 

A  A  A3 

10 

0000 

0051 

ONE 

113 

0046 

10 

0000 

00  5  1 

A  A  AS 

50 

0000 

0081 

F  1  VE 

114 

0801 

SO 

00  0  0 

00  5  1 

A  A  A? 

14 

B  413 

1653 

E  T  0  5 

115 

0016 

14 

84  13 

16  5  3 

A  A  A9 

87 

1838 

1851 

E 

116 

0B16 

87 

183  6 

18  5  1 
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8  X  A 

8001 

CHE  CK 

394 

1418 

8  X  A 

0001 

CELLS 

395 

1274 

eu  A 

LOD 

CON  T5 

RF.UMNING 

396 

1080 

C  1  N02      LOOP  3 

3  97 

108  3 

C  ONTS     RAL 
6LO 

STL 

CELLS 
TWO 

3  96 
399 

0  83  4 
1313 

CELLS 

400 

1655 

Mil 

CO  N  T  6 

401 

1262 

C  1  H  Dl 

402 

1116 

ALO 
8  TL 

TWO 

403 

1363 

C  1  N  D  1 

ADJUST  TO 

404 

1705 

C  1  N  02 

NEXT   CELL 

405 

1312 

ALO 
STL 

TWO 
C  1  N02 

SET 

406 
407 

1413 
1755 

LDO 

C  ONT«     RAW 

8CT 

STL 

RAL 
ALO 
SL  0 
RAU 

C  1  N  Dl      L00P3 
DATA 

0000 
TEMPI 

CONVERT 

408 
409 
4  10 
411 

136  2 
1067 
1411 
1133 

8  003 

ROWS  TO 

4  12 

1070 

SIXTY 

FLOATING 

4  13 

1  477 

T  E  UP1 
8002 

POINT 

4  14 

4  15 

1808 
1271 

600  3 

1009 

1012 

120  1 

800  2 

117  1 

091  7 

0967 

210  0 

1252 

240  0 

12  5  3 

000  1 

1059 

106  6 

1163 

13  5  7 

1061 

09  5  0 

0070 

10  5  0 

1160 

1165 

1068 

116  5 

1118 

BOO  1 

107  4 

24  0  0 

130  3 

250  0 

135  3 

00  0  1 

1109 

107  4 

1313 

116  5 

12  19 

110  0 

0920 

115  0 

1210 

000  4 

122  1 

175  0 

1305 

095  0 

1111 

095  8 

116  1 

135  7 

1211 

OSE  0 

0970 

10  5  0 

1360 

800  2 

1017 

80  0  1 

0956 

999  9 

1303 

260  0 

14  0  3 

000  1 

1159 

095  b 

1263 

116  5 

1269 

095  0 

13  10 

1215 

1168 

085  0 

13  55 

1209 

1062 

116  5 

1319 

090  0 

1020 

09  5  0 

1360 

08  5  0 

140  5 

13  5  9 

1112 

13  5  7 

126  1 

10  50 

1410 

080  0 

14  55 

130  9 

1162 

130  9 

1213 

600  1 

1216 

250  0 

1505 

160  0 

13  52 

090  6 

08  33 

090  6 

1359 

000  1 

126  5 

250  0 

1555 

000  2 

1311 

250  0 

1177 

000  1 

0883 

1315 

1268 

800  1 

113  4 

350  0 

1237 

1215 

1318 

900  1 

1174 

600  1 

0980 

2500 

1377 

140  0 

1402 

090  6 

09  33 

090  6 

1409 

000  1 

1315 

2500 

1605 

000  2 

1361 

3500 

1327 

121  5 

1368 

800  1 

1324 

800  1 

1030 

250  0 

1377 

000  2 

0963 

2500 

1427 

090  6 

10  3  3 

090  6 

1459 

1215 

1418 

600  1 

1374 

000  1 

1060 

108  3 

08  34 

135  9 

1212 

1309 

1313 

085  0 

1655 

130  9 

1262 

1116 

1067 

120  9 

1363 

085  0 

1705 

120  9 

13  12 

12  5  9 

14  13 

08  5  0 

1755 

1259 

1362 

120  9 

1212 

13  5  7 

1411 

0000 

1133 

091  7 

1070 

600  3 

1477 

1200 

180  5 

091  7 

127  1 

800  2 

0929 

136 


FOV 

f  PS 

416 

0939 

34 

1450 

1453 

F  36 

FP1 

417 

1452 

33 

12  50 

1527 

BTU 

ROWS 

418 

1527 

21 

006  2 

1135 

R  *  Lt 

419 

1135 

60 

0000 

1855 

rov 

f  P  2 

430 

1B55 

34 

1300 

1502 

F  O  V 

RONS 

4  31 

1502 

34 

008  3 

0832 

rg  p 

FPB 

433 

0833 

39 

1500 

1553 

FOV 

F  P9 

433 

1S52 

34 

1550 

1602 

FHP 

1  N  TEG 

4  34 

1602 

39 

090  6 

1006 

8  TU 

1  N  TEC 

CON  T7 

435 

1006 

21 

090  6 

1509 

C0HT7 

Rii 

0  A  T  A 

4  26 

1  509 

65 

135  7 

14  6  1 

D  1  V 

B  t  X 

437 

1461 

14 

1050 

14  60 

STL 

1  NOEX 

BET 

4  38 

1460 

20 

136  5 

1468 

8LT 

0O04 

DATA 

4  39 

1468 

35 

000  4 

097  9 

ALO 

1  N  0  A 

AODRE  8SE  8 

430 

0979 

15 

180  0 

1056 

LOO 

A  TlO 

FOR 

4  31 

1056 

69 

1559 

1412 

BOA 

*TiO 

TRAP  RULE 

433 

1412 

32 

1559 

1462 

LOD 

ft  TREE 

433 

1462 

69 

14  15 

1516 

SO  A 

ftTREE 

434 

1518 

22 

1415 

156B 

R  At 

IN  DB 

435 

156B 

65 

165  0 

1106 

LDO 

BONE 

436 

1106 

69 

160  9 

1512 

SO  A 

BONE 

437 

1512 

32 

160  9 

1563 

LOO 

BT  BO 

438 

1563 

69 

146  5 

1616 

BOA 

BT  WO 

439 

1618 

22 

146  5 

1668 

RAU 

B 

4  40 

1  668 

60 

000  1 

1156 

FO  V 

FP  2 

4  41 

1156 

34 

1300 

16  52 

FSB 

B  P  R  1  M 

4  42 

1652 

33 

000  3 

1029 

FMP 

ROWS 

443 

1039 

39 

006  2 

0862 

FOV 

FP  4 

444 

0882 

34 

14  0  0 

1702 

is  T  U 

F  A  C  TR 

44  5 

1702 

31 

0096 

06  4  9 

LOO 

INDEX 

446 

0849 

69 

136  5 

17  18 

R  A  A 

8  001 

447 

1718 

80 

800  1 

1324 

1  AB 

8  001 

A  T  WO 

44  6 

1324 

83 

800  1 

15  59 

ATWO 

RAU 

9999 

BONE 

449 

1559 

60 

999  9 

1609 

BONE 

FAO 

9999 

4  50 

1609 

32 

999  9 

0975 

81  A 

0001 

CALCULATE 

4  51 

0975 

51 

000  1 

118  1 

3XB 

0001 

TRAP  RULE 

4  52 

1161 

53 

000  1 

0637 

H  7   A 

A  TREE 

CON  T8 

INTEGRAL 

453 

0  8  3  7 

40 

1415 

079  1 

ATREE 

FAO 

9999 

BTWO 

4  54 

1415 

32 

999  9 

14  65 

BTWO 

FAO 

9999 

455 

1465 

32 

999  9 

1025 

FUP 

F  A  CT« 

4  56 

1035 

39 

009  6 

084  6 

FAO 

1  N  TEG 

4  57 

0846 

33 

090  6 

1183 

BTU 

t  N  TEG 

A  TlO 

4  58 

1183 

21 

090  6 

1559 

C  ONTB 

RAU 

1  N  TEG 

4  59 

0791 

60 

090  6 

1511 

FMP 

FP  4 

4  60 

1511 

39 

1400 

1753 

STU 

1  N  TEG 

4  61 

1752 

21 

090  6 

1659 

RAU 

M 

4  62 

1659 

60 

08  5  6 

1561 

FSB 

H 

463 

1561 

33 

067  5 

1802 

NZ  E 

CON  T9 

464 

1803 

45 

120  6 

1407 

RAU 

1  N  TEG 

465 

1206 

60 

090  6 

1611 

FMP 

FP  2 

466 

1611 

39 

130  0 

18  52 

BTU 

1  N  TEG 

CON  T9 

467 

1652 

31 

090  6 

1407 

C0NT9 

AXC 

0001 

468 

1407 

58 

000  1 

14  6  3 

LOO 

1  N  TEG 

469 

1463 

69 

090  6 

170  9 

STO 

FUN     G 

470 

1709 

24 

6770 

087  3 

LOO 

8007 

COUNT 

471 

0873 

69 

800  7 

1079 

8T0 

T  E  RMS 

TERMS 

473 

1079 

34 

093  2 

118  5 

RAU 

HARM 

473 

1  1B5 

60 

000  3 

1457 

FS6 

M 

474 

1457 

33 

085  6 

123  3 

NZE 

CO  N  10 

475 

12  33 

45 

093  6 

0887 

RAU 

476 

0936 

60 

085  6 

1661 

FAO 

fp  a 

AO JUST 

477 

1661 

33 

130  0 

1577 

STU 

M 

CON  Tl 

M   A  Nil    N 

478 

1577 

21 

08  5  6 

0628 

C  OHIO 

RAU 

HARM 

479 

0887 

60 

000  3 

1507 

FSB 

N 

4  DO 

1507 

33 

067  5 

1453 

NZ  E 

CON  1  1 

481 

14  53 

45 

1356 

1557 

RAU 

N 

482 

1356 

60 

067  5 

1129 

FAO 

fp  a 

4  83 

1139 

32 

130  0 

1637 

STU 

M 

4  84 

1637 

31 

085  6 

1759 

BTU 

N 

C0NT1 

4  B5 

1  759 

21 

067  5 

0828 

C  ONI  1 

LOO 

FIRST 

4  86 

1557 

69 

0004 

1607 

BTO 

GAM11 

4B7 

1607 

24 

1510 

1513 

STO 

GAMMA 

C0N33 

488 

1513 

34 

073  0 

1263 

C0N83 

LOO 

FP  1 

489 

1383 

69 

1250 

1503 

STO 

M 

490 

1503 

24 

085  6 

1609 

STO 

N 

4  91 

1809 

24 

067  5 

087B 

RAU 

P  1 

4  92 

0878 

60 

1350 

1306 

FOV 

A 

493 

1306 

34 

0000 

155  3 

FMP 

6003 

4  94 

1553 

39 

800  3 

1657 

BTU 

A  LPH3 

4  95 

1657 

21 

1612 

1515 

RAU 

P  1 

496 

1515 

60 

1350 

1356 

FOV 

B 

497 

1356 

34 

0001 

1603 

FMP 

9003 

4  98 

1603 

39 

BOO  3 

1707 

STU 

BETAS 

499 

1707 

31 

166  2 

1565 

RAU 

G  AU11 

500 

1565 

60 

1510 

1615 

FMP 

6003 

501 

1615 

39 

800  3 

1369 

FSB 

A  LPH3 

CALCULATE 

502 

1369 

33 

161  2 

0B89 

FSB 

BE  T  A3 

KAPPA 

503 

0889 

33 

166  3 

0939 

STU 

K  A  P60 

SOU  ARE  0 

504 

0939 

21 

004  4 

1397 

RAC 

0001 

505 

1397 

B8 

000  1 

1653 

LOD 

a  Awn 

506 

1653 

69 

1510 

1563 

STO 

GAMMA   C 

507 

1563 

24 

6730 

13  33 

AXC 

0001 

508 

1333 

58 

000  1 

0969 

RAU 

M 

509 

0989 

60 

085  6 

1711 

FAO 

FP2 

LOOP  A 

510 

1711 

33 

1300 

1677 

LOOPA 

FMP 

P  1 

511 

1677 

39 

1350 

1703 

FOV 

A 

512 

1703 

34 

0000 

1753 

FMP 

8003 

513 

1753 

39 

BOO  3 

1757 

BTU 

A  LPH3 

514 

1757 

31 

1612 

166S 

RAU 

N 

515 

1665 

60 

067  5 

1179 

FMP 

P  1 

516 

1179 

39 

1350 

1603 

FOV 

B 

517 

1803 

34 

0001 

1853 

FMP 

BOOS 

518 

1853 

39 

800  3 

1807 

137 


FAD 

K  A  P9Q 

F  *0 

AL  PH3 

LOO 

EOO  AU 

S  TO 

GAMMA   C 

A  X  C 

OOOl 

R1U 

HARM 

FSB 

U 

N  I  £ 

C0N13 

ft  *  U 

N 

FAD 

FP  3 

sru 

M 

LOOP  A 

C0N1  S 

RIU 

HARM 

FSB 

N 

NZE 

CON  1  3 

R  A  U 

N 

FAD 

F  P  2 

S  T  if 

STU 

i 

LOOP  A 

C0N13 

nzb 

CO  N  15 

LDO 

TERMS 

RA  C 

8001 

LOOPB 

L  0  OPB 

RSU 

C 

F  AO 

ELEV 

F  UP 

GAMMA   C 

STU 

ARC 

R  AL 

ARC 

LOO 

EOOE  A 

3TL 

TEMPI 

ft  AM 

A  R  G 

LOO 

EOOE  A 

R  A  U 

8002 

FSB 

TEMPI 

FMP 

FMP 

| 

STU 

TEMPI 

R  A  U 

F  P  8 

FMP 

FMN     C 

FD  V 

TEMPI 

STU 

A  M  N     C 

3  X  C 

000  1 

NZ  C 

L  0  0  P  B 

C  0  N  1  4 

CONK 

R  A  C 

00  0  B 

LO  □ 

L  0  OPC 

E  00  CL 

L  OOPC 

LDD 

A  UN     C 

S  TO 

1976   C 

8XC 

0001 

NZ  C 

L  0  OPC 

PCH 

1977 

R  AL 

TERMS 

SL  0 

EIGHT 

BU  1 

CO  N  15 

LOD 

TERMS 

R  A  C 

8001 

C0N25 

C  0N2  5 

LOO 

A  MN     C 

STD 

1968   C 

BXC 

0009 

N  Z  C 

C  0  N  3  4 

PCH 

197  7 

C0N1  5 

C0N24 

A  X  C 

0008 

C0N25 

C0N15 

LDO 

POINT 

RA  A 

B001 

R  IB 

0001 

LOOP  0 

L  OOPD 

LOO 

TERMS 

R  A  C 

8001 

LDO 

ZERO 

S  TO 

CO  F  AC 

LOOPE 

LOOPE 

R  A  U 
FMP 

STU 

GAM  11 

TEMPI 

R  AU 

C 

FSB 

z      ' 

FMP 

GAMMA   ( 

FAO 

TEMPI 

R  AL 

8003 

LOO 

root  * 

R  A  It 

eooa 

FMP 

A  MN     C 

FO  V 

FP  3 

STU 

TEMPI 

RSU 

C 

FAO 

Z       A 

FMP 

GAMMA   C 

FMP 

FP  3 

R  AL 

B003 

LDO 

EOOE  A 

RSU 

8003 

FAD 

f  P  1 

STU 

END  CO 

FUP 

TEMPI 

FAO 

COP  AC 

STU 

COF  AC 

8  X  C 

0001 

NIC 

LOOPE 

C0N16 

C  ON  16 

R  AU 

N 

FO  V 

CO  F  AC 

STU 

N  COR 

LDD 

E  00  CL 

LDO 

Z 

3  T  0 

1977 

LOO 

N 

STO 

1978 

10  0 

N  COR 

9  TO 

1979 

519 
520 
CALCULATE  52  1 
GAMMA  MN  522 
523 
534 
5  25 
526 
527 
52P 
539 
530 
531 
5  32 
533 
534 
535 
536 
5  3? 
5  J  a 
539 
540 
541 
542 
54  3 
54  4 
5  4  5 
546 
547 
548 
549 
550 
I  S  1 
553 
5  5  3 
5S4 
555 
CALCULATE  556 
AMN  557 

558 
559 
560 
561 
562 
563 
564 
PUNCH  565 
CONSTANTS  566 
567 
56  B 
569 
570 
571 
572 
573 
574 
575 
576 
577 
STB 
579 
580 
5  Bl 
5B2 

5  8  J 
58* 
585 
586 
587 
588 
589 
590 
591 
592 
593 
594 
595 
596 
597 
598 
599 
600 

6  01 
6  0  2 
603 
604 

CALCULATE  605 
END  AND  606 
HARMONIC  607 
CORRECTION  608 
609 
6  10 
6  1  1 

cia 

61  3 
6  14 
6  15 
6  16 
6  17 
618 
DATA  619 


CALCULATE 
COR  Rt  C  TE  D 
OAT  A 


1807 
1321 
1039 
0  94  2 
1383 
1089 
1857 
143  3 
09B6 
176  1 
1737 
0937 
1008 
150  4 
10  SB 
1329 
1777 
1560 
18  59 
16  13 
12  35 
084  1 
1811 
1385 
1130 
0987 
1139 
0992 
1120 
1189 
1042 
155  4 
084  3 
1604 
1654 
1170 
1406 
1220 
1117 
1704 

16  10 
0864 
127  0 

0  923 
1754 
1279 
1335 
1339 
1827 

10  37 
14  56 
1660 

1  3G5 
089  1 
1804 
1371 
1877 
1331 

11  BO 

17  13 
1108 
0914 
1320 
1435 
094  1 
1854 
1861 
1715 
1  506 
1  370 
1762 

0  938 
1230 
0893 
1556 
1710 
1419 
1606 
1656 
1430 
1812 
0978 
1280 
1706 
1663 
1166 
1075 
1028 

1  4B  5 
1167 
1535 
186  2 
1768 
0922 
1  339 
1208 
1756 
1760 
1806 
1330 
1078 
1381 
1856 


33  004  4 

33  1612 

69  094  2 

31  673  0 
58  000  1 
60  0003 
33  08  5  6 
45  098  6 
60  OB  5  6 

32  1300 
21  0856 
60  000  3 

33  067  5 
45  105  8 
60  067  5 
32  130  0 
31  OB  5  6 


21 


067  5 
171  2 
093  3 
800  1 
000  7 

090  8 
673  0 
088  4 
088  4 
09  9  2 

091  7 
088  4 
104  2 
80  0  2 
0917 
0000 
000  1 
0917 
1500 
6770 
0917 
675  0 
000  1 
084  1 
000  8 

092  3 
67  5  0 
797  6 
0001 
092  3 


1712 
093  3 
800  1 
67  5  0 
796  8 
00  0  9 
118  0 
197  7 
000  8 
000  5 
800  1 
00  0  1 
093  2 
BOO  1 

00  5  0 

1  1  S  8 
60  1  SI  0 
39  2650 
21  0917 
60  00  0  7 

33  265  0 
39  673  0 
32  0917 
65  BOO  3 
69   1710 

60  SO  0  2 
39   675  0 

34  130  0 
21   0917 

61  0007 
32  2650 
39  6730 
39  130  0 
65  80  0  3 
69  116  6 
61  800  2 
32  125  0 
21  098  2 
39  0917 
32  115  8 
21   1158 

59  000  1 
4B  186  1 

60  267  5 
34  11 5  6 
31  270  0 
69  176  0 
69  265  0 
24  197  7 
69  26  7  5 
24  197  8 
69  270  0 
34   197  9 


69 


133  1 
10  39 
150  1 

13  83 
10B9 
1857 

14  33 

09  37 
176  1 
172  7 
1677 
1006 
1504 

10  59 
1229 

17  77 
156  0 

16  77 
1613 
123  5 
08  4  1 

18  11 
128  5 
1130 
0987 
1139 
00  5  1 
1120 
1189 
0051 

15  5  4 
0B4  3 
1604 
1654 
1170 
1406 
1220 
1117 
1704 
1610 
08  6  4 
1270 
115  1 

17  54 
127  9 

13  35 
1239 
1R27 
10  37 

14  56 
166  0 
13  8  5 
069  1 
160  4 
137  1 
1877 
1231 
1712 
089  1 
1108 
0914 

13  20 

14  35 
094  1 

18  5  4 
186  1 
1715 
1506 
1370 
1763 
0928 
1230 
089  3 
1556 
0051 
1419 
1606 
1656 
1420 
1812 
097  6 
1280 
1706 
1663 
0051 
1075 
1028 
1485 
1167 

15  35 
1862 
1768 
0922 
1329 
1208 

17  56 
1151 
1806 
13  30 
1078 
1281 

18  56 
10  32 


138 


LD  D 

COF  AC 

FACTOR 

622 

10  3  3 

69 

1158 

1713 

3  TO 

19B0 

633 

1713 

34 

I960 

1483 

LOO 

END  CO 

624 

148  3 

69 

098  3 

1585 

8  TD 

1901 

635 

1585 

14 

198  1 

0934 

PC  h 

1977 

636 

0934 

71 

197  7 

USB 

9H 

0001 

6  37 
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APPENDIX   G 


IBM  650  Program  for  Experimental 
Determination  of  Effective  Pile  Size 


This  program  was  written  to  determine  the  effective  pile  size  by 
an  experimental  analysis  of  data  taken  in  a  horizontal  traverse.     The 
program  was  written  in  SOAP  II  and  floating  point  form.     The  object 
program  and  a  logic  diagram  are  included  in  this  section. 

The  program  is  based  upon  the  equation 


*  W  =   £ 
m=l 
and  proceeds  as  follows: 

1.  Using  an  assumed  initial  value  for  the  effective 

pile  size,     a,     the  constants    A        are  evaluated  using  a  least 

m  ° 

squares  analysis. 

2.  The  least  squares  error  is  calculated  using  the 

equation    ,. 

V*     !■  »x.  3irx.  n2 

L,        Al  cos      a*  +  A3  cos —    +  •••    -  N.      . 

i=  1 

3.  The  value  of  a    is  increased  by  an  increment    Aa 

and  a  new  set  of    A    '  s    is  calculated. 
m 

4.  A  new    EJ    is  calculated  and  compared  with  the  old 

E  .     If    E      is  found  to  be  decreasing,     a    is  increased  again  by  Aa. 

5.  If    E*   is  increasing,     a    is  decreased  by    Aa    and 
Aa    is  divided  by  two.     This  process  is  repeated  until    Aa    is  less 
than  a  specified  precision.     The  last  value  of    a    is  then  the  desired 
value  of  extrapolated  pile  size. 


(G-l) 
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The  criterion  for  determination  of    a    is  that    E*,     described 
above,   should  be  minimized.     Input  to  this  program  is  entirely  on  one- 
word  load  cards  and  is  listed  in  Table  29. 


Table  29.     Input  to  IBM  650  code  for  deter- 
mination of  effective  size. 


Symbol 


Explanation 


Drum  Storage 
Location 


AINIT 
DELA 
INDXA 

ENN 
N 


Initial  value  of    a 

Increment      Aa 

Number  of  Data  Points 
form    -    00    0000    OOxx 

Number  of  Harmonics 
form    -    00    0000    OOOx 
(2,   3,   4,  or  5) 

BF3  Count  Rates    -    to  be 
stored  consecutively  starting 
with  location 

Co-ordinates  of  count  rates, 
to  be  stored  consecutively  in 
location 


1000 
1001 
1002 

1003 


0301 


0201 


The  capacity  of  this  program  is  100  data  points.     It  will  perform 
the  analysis  using  from  two  to  five  terms  of  Eq.   (G-l).     Output  forms 
for  this  program  are  listed  in  Table  30.     Form  one  is  punched  after 
each  trial  and  form  two  is  the  final  output. 


Word  1 


Form  2: 


Table  30.     Output  forms  for  IBM  650  program 
for  effective  pile  size. 


Form  1: 

A,  A, 


Precision       Number  of 
Data  Points 
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The  initial  estimate  of  a    should  always  be  less  than  the  final 
answer  or  the  program  will  not  converge.     Therefore,  the  actual 
size  is  recommended  as  an  initial  estimate  since  this  will  always  be 
less  than  the  extrapolated  size.     Approximately  45  seconds  is  required 
for  one  trial  using  nine  data  points. 
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OBJECT  PROGRAM   -   APPENDIX  G 
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APPENDIX   H 

IBM  650  Program  for  Diffusion  Length; 
Monoenergetic  Fast  Source 

This  code  was  written  to  determine  an  experimental  value  of 
diffusion  length  based  upon  the  assumption  of  a  monoenergetic  source 
of  fast  neutrons.     The  program  was  written  in  SOAP  II  and  floating 
point  form.     The  object  program  and  a  logic  diagram  are  given  in 
this  section. 

The  solution  to  this  problem  is  based  upon  the  equation 
««  q  -k2t["-T       z 


ab  Dy 
,  ,  mn 

m=l    n=l 


2nTT 


+  e    mnZ    [i  -«*£(_*_  +   y       «/—)] 


-2v      (c-z) 

/,  mn1        '  , 

(1  -  e  ). 


An  iteration  procedure  similar  to  that  used  in  the  point  thermal  and 
constant  thermal  source  programs  was  used. 

Input  parameters  and  data  are  listed  in  Table  31.     All  input 
should  be  punched  on  one  word  load  cards  and  read  into  the  machine 
with  the  object  deck. 

Output  from  this  program  is  in  two  forms,   as  shown  in  Table  32. 
Form  one  is  punched  after  the  calculation  of  each  individual  cor- 
rection factor.     Form  two  is  the  final  result  and  gives  the  last  value 
of  7ii    as  well  as  the  diffusion  length. 

The  operating  time  for  the  calculation  of  one  correction  factor 
using  four  harmonics  is  approximately  40  seconds.  The  capacity  of 
the  program  is  200  data  points. 


Table  31.     Input  to  IBM  650  program  for  diffusion 
length,     monoenergetic  fast  source. 
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Symbol 


Explanation 


Drum  Storage 
Location 


FIRST  Initial  estimate  of    yn 

POINT  Number  of  Data  Points,   in 

form  00    0000    OOxx 

DATPT  Number  of  Data  Points,   in 

floating  point  form 

A  Extrapolated  x-dimension 

B  Extrapolated  y-dimension 

C  Extrapolated  z-dimension 

X  X    co-ordinate  of  data 

Y  y    co-ordinate  of  data 

TAU  Fermi  age  of  source  neutrons 

HARM  Magnitude  of  highest 

harmonic  used. 


0100 
0101 

0102 

0103 
0104 
0105 

0106 

0107 

0108 

0109 


Table  32.     Outpxit  forms  for  IBM  650  program  for  dif- 
fusion length,   monoenergetic  fast  source. 


Word  No.     1 


Form  One: 
z  co-ordinate    data 


Form  Two: 
Yn,  initial 
value 


corrected     Harmonic  End  Diffusion 

data  correction       correction        Length 


Yn  last        Precision 
value 


POINT 


Yu,  next  to      Diffusion 
last  value  Length 
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APPENDIX  I 


IBM  650  Program  for  Diffusion  Length, 
Gaussian  Range  Empirical  Source 


This  program  was  written  to  determine  the  diffusion  length  based 
upon  an  empirically  determined  source  of  high  energy  neutrons.     The 
program  was  written  in  SOAP  II  and  floating  point  form.     The  object 
program  and  a  logic  diagram  are  included  in  this  section. 

The  solution  to  this  problem  is  based  upon  the  equation 

K2r.2 
•» .    »,      f.  Q                        '       -7       zr                                 y         r.     - 
„\       \       \         l      o                  4              'mn         ,   ,        ,,  z              mn     i   \\ 
»,*.«) -£     I       ahD  e  e  [l  +  erf(— g )j 


m=l    n=I 


+  7 
.           mn        ,            ,  ,  z  mn     i 

+  e  1   -  erf  (  +  ■ 


[l   -  erf  ( 


-2y       (c-z) 
t,  mnv        '. 

(1   -  e  ) 


Correction  factors    C       and    C       derived  from  this  equation  are  listed 
E  ri 

in  Table  1.     A  correction  factor  iteration,    such  as  has  been  described 
previously,   is  used  to  obtain  an  accurate  value  of    Yn    an<i  hence  the 
diffusion  length. 

Input  parameters  and  data  are  listed  in  Table  33.     All  input  should 
be  punched  on  one  word  load  cards  and  read  into  the  machine  with  the 
object  deck. 

Output  forms  are  listed  in  Table  34.     Form  one  is  punched  after 
calculation  of  each  individual  correction  factor.     Form  two  gives  the  end 
result  of  the  analysis. 

The  operating  time  for  the  calculation  of  one  correction  factor  using 
four  harmonics  is  approximately  one  minute  and  30  seconds.     The  data 
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capacity  of  the  program  is  200  data  points.     As  many  as  five  sets  of 
Gaussian  ranges  and  fractions  may  be  used. 


Table  33.     Input  parameters  and  data  for  IBM  650 
program  for  diffusion  length,   empirical 
source  method. 


Symbol: 


Explanation 


Drum  Storage 
Location 


GAMIN        Initial  Estimate  of    Tu 

POINT         No.   of  Data  Points,   form 
00  0000  OOxx 

DATPT       No.   of  Data  Points, 
floating  point  form 

A  Extrapolated  x-dimension 

B  Extrapolated  y-dimension 

C  Extrapolated  z-dimension 

HARM  No.    of  highest  harmonic, 

floating  point  form 

GAUSS         No.    of  Gaussian  Ranges  used, 
form    00  0000  OOOx 

QO  Source  strength 

DIFCO         Diffusion  Coefficient 

RI  Gaussian  Range,   stored  con- 

secutively starting  at 

FI  Range  Fractions,    stored  con- 

secutively starting  at 

N  Count  Rates,   stored  con- 

secutively starting  at 

z  Data  co-ordinates,   stored 

consecutively  starting  at 


0100 
0101 

0102 

0103 
0104 
0105 

0106 

0107 

0108 
0109 

0901 

0906 

0501 

0301 


Table  34.     Output  forms  for  IBM  650  program  for  determination 
of  diffusion  lenth,   empirical  source  method. 


WORD  1 


Form  1: 

z 
Form  2: 
Gamin 


N 


N  (corrected) 


Tii,   last        Precision 
value 


Combined  cor-  7n 
recti  on 

No.   of  Data  Diffusion 

Points  Length 
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A  L  PH2 

STU 

A  L  PH2 

R  A  U 

P  1 

FO  V 

e 

8TU 

If  TA2 

FMP 

BE  T  A3 

STU 

BE  T  A3 

RAU 

GAM  30 

FSB 

ALPH2 

FSB 

BE  T  A3 

LOO 

EOO  AU 

STU 

KAPPA 

RAU 

ON  E 

FD  V 

KAPPA 

FMP 

CON  V  T 

STU 

DL 

LOO 

GAMIN 

S  TO 

1977 

LOO 

NKAU 

STO 

1978 

LOO 

PR  E  C 

STO 

1979 

LOO 

POINT 

STO 

I960 

LOO 

OL 

STO 

1981 

PCH 

1977 

eo  oo 

413 

1477 

60 

370  0 

1655 

414 

1  655 

69 

020  B 

1100 

41S 

0308 

39 

2300 

1950 

416 

1950 

32 

01B0 

1107 

417 

1107 

31 

0180 

0283 

418 

0283 

51 

000  1 

138  9 

419 

1389 

40 

018  3 

024  3 

4  10 

034  3 

60 

026  2 

0917 

481 

0917 

39 

096  8 

1018 

433 

1018 

21 

0012 

16  15 

423 

1615 

60 

010  2 

1157 

424 

1157 

39 

0180 

0230 

425 

0230 

21 

006  2 

1665 

426 

1665 

60 

010  2 

1207 

427 

1207 

39 

0274 

0934 

420 

0924 

21 

102  8 

1231 

429 

1231 

60 

096  8 

0273 

4   30 

0273 

39 

096  8 

1068 

431 

1068 

31 

122  2 

0125 

432 

0125 

60 

102  8 

0933 

435 

0933 

33 

132  2 

1049 

4  34 

1049 

21 

102  B 

128  1 

435 

1281 

60 

001  2 

0967 

4  36 

0967 

33 

006  2 

14  39 

4  37 

1439 

34 

102  8 

1078 

438 

1078 

21 

133  2 

1585 

439 

1585 

34 

1000 

1301 

440 

1301 

21 

1.106 

1259 

441 

1259 

60 

110  6 

1311 

442 

1311 

33 

133  2 

1309 

443 

1  309 

67 

800  3 

1017 

444 

1017 

60 

800  2 

0175 

445 

0175 

33 

130  6 

0983 

446 

0983 

46 

023  6 

1187 

447 

1187 

60 

133  2 

1337 

44  8 

1237 

31 

1106 

0309 

449 

0236 

60 

133  2 

1387 

GAMMA 

4  50 

1287 

39 

133  3 

13B2 

SOU  AR  E  0 

451 

1382 

21 

0386 

1489 

4  52 

1489 

60 

0250 

1705 

4  53 

1705 

34 

010  3 

1703 

4  54 

1703 

21 

0112 

1715 

ALPHA 

455 

1715 

39 

0112 

0912 

3  iU Aft  E  D 

4  56 

0912 

21 

0113 

1765 

4  57 

1765 

60 

0250 

1755 

458 

1755 

34 

0104 

1504 

459 

1504 

31 

016  2 

1815 

BETA 

4  60 

1815 

39 

016  2 

0962 

SQUARED 

461 

0962 

31 

016  2 

1S65 

462 

1865 

60 

0386 

094  1 

KAPPA 

463 

094  1 

33 

0113 

15  39 

SOU  ARE  0 

464 

1539 

33 

016  2 

158  9 

465 

1589 

69 

099  2 

1350 

466 

0993 

21 

014  6 

1099 

467 

1099 

60 

0050 

1805 

468 

1805 

34 

014  6 

0196 

0  1  F  FU8  1  ON 

469 

0196 

39 

0950 

1351 

LENGTH 

470 

1351 

31 

1356 

1359 

471 

1359 

69 

0100 

1753 

472 

1753 

24 

1977 

0280 

473 

0280 

69 

1332 

1635 

474 

1635 

24 

197  8 

1331 

475 

1331 

69 

1306 

1409 

476 

1409 

24 

1979 

14  33 

477 

1432 

69 

010  1 

1554 

478 

1554 

24 

1960 

1033 

479 

1033 

69 

13S6 

1459 

480 

1  459 

24 

1981 

0834 

481 

0234 

71 

197  7 

8000 
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APPENDIX   J 

IBM  650  Program  for  Determination  of  Material 
Buckling  and  Effective  Multiplication  Factor 

This  code  was  written  to  determine  the  material  buckling, 

B      .   and  the  effective  multiplication  factor,     k   ,,,     in  the  KSU 
m  '  r  eff 

exponential  pile.     The  program  was  written  in  SOAP  II  and  floating 
point  form.     The  object  program  and  a  logic  diagram  are  included 
in  this  section. 

The  program  is  based  upon  the  point  thermal  source  solution 
to  the  thermal  diffusion  equation.     An  accurate  value  of    Yu    is  deter- 
mined by  the  iteration  procedure  outlined  previously  using  the  cor- 
rection factors  listed  in  Table  I.      B    z    is  calculated  from  the  equation 

m 

B    2    =    (-f-)2    +    (-£-)*-  Y„2. 
m  a  b 

The  program  then  calculates  the  infinite  multiplication  factor,     k   , 

oo 

from 

k     =   1  +   B   2  (L2  +  r). 

oo  m  t 

The  effective  multiplication  factor,    k   ,,,     can  then  be  calculated  from 
-B   2  t  eif 

k      e      G 


keff 


1  +    B, 


where    B2    is  the  geometric  buckling.     By  increasing  the  effective 
dimensions  of  the  multiplying  medium,     k  „    can  be  extrapolated  to 
critical  pile  size.     This  is  done  by  reading  in  new  values  for    a,  b, 
and    d     and  starting  the  program  over  at  CON11    by  means  of  a 
transfer  card. 


Input  forms  for  this  program,  are  given  in  Table  35.     All  input 
should  be  in  the  form  of  one  word  load  cards. 


Table  35.     Input  forms  for  IBM  650  program  for  material 
buckling  and  effective  multiplication  factor. 
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Symbol 


Explanation 


Drum  Storage 
Location 


FIRST                Initial  estimate  of    Tn  0100 

POINT               Number  of  Data  Points,  0101 
form  00  0000  OOxx 

DATPT             Number  of  Data  Points,  0102 
floating  point  form 

A                         Extrapolated  x-dimension  0103 

B                           Extrapolated  y-dimension  0104 

C                         Extrapolated  z-dimension  0105 

X                         x  co-ord  of  data  0106 

Y                         y  co-ord  of  data  0107 

L                           Lattice  diffusion  length  0108 

TAU                     Fermi  age  0109 

D                          z  dimension  of  active  lattice  0110 

N  Count  Rates,  to  be  stored 

consecutively,    starting  at  0501 

z  Data  co-ordinates,  to  be  stored 

consecutively,   starting  at  0301 


Output  from  this  program  is  in  three  forms.     Form  one  is 
punched  after  calculation  of  each  correction  factor.     Form  two  is  the 
results  of  the  buckling  calculation.     Form  three  gives  the  buckling 
and  the  multiplication  factors.     Table  36  lists  the  output  forms. 


Table  36.     Output  forms  for  IBM  650  program  for  material 
buckling  and  effective  multiplication  factor. 


WORD    1 
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Form  1: 


Form  2: 


Form  3: 


N 


N  (corrected)        C„ 


7„,   last        Tn,  next  to  Tfti   initial        No.   of  Data 

value  last  value  value  Points 


eff 


The  operating  time  necessary  to  calculate  the  correction  factor 
for  one  data  point  using  four  harmonics  is  approximately  45  seconds. 
The  capacity  of  the  program  is  200  data  points. 
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BLR 

19  51 

BLR 

1977 

6LH 

0  300 

'■;  *  n 

F  (  ft  ST 

19SB 
19  64 
09  00 
01  00 
0101 
0103 
0103 
0104 
0105 
01Q6 
0101 
0108. 
0109 

offfl 

0300 
0500 


8YN 

H  C  OR 

07  00 

8  YN 

START 

1999 

ZERO 

0  0 

0000 

0000 

ONE 

10 

000  0 

0051 

TO  0 

no 

0000 

00  5  1 

MT  NO 

0000 

00  51 

THREE 

30 

0  0  00 

0051 

P  1 

31 

4  160 

00  5  1 

9  E  YEN 

70 

0000 

0051 

C  ON  VT 

25 

3600 

00  51 

10000 

10 

00  00 

0055 

E  OOCL 

S  TD 

Z  Z  Z  1 

LOO 

Z  Z  Z  10 

&  TO 

1977 

STO 

1  97B 

N  TD 

1S79 

STU 

I96  0 

STO 

198  1 

&  TD 

1983 

STO 

19S  3 

STO 

198  4 

I  11   1 

Z  Z  Z  10 

00 

00  0  0 

0000 

B  TL 

LNX05 

STL 

L  N  X  1  1 

S-L  T 

0  0  08 

SUP 

F  1  F  T  V 

Nil 

LNX  04 

|.  NX  03 

KM  1 

R  SU 

B  0  0  3 

LO  D 

F  P  ONE 

STO 

LNXOH 

i-N<OI 

L  NX03 

a  r  t 

000  9 

H  T 

0000 

A  UP 

9  t  XT  Y 

SUP 

8  0  0  3 

R  A  U 

8003 

FMP 

L  N  X02 

FMP 

INTEN 

STU 

L  N  X05 

INHM 

LNX04 

R  AL 

L  N  X09 

SRT 

0003 

R  A  U 

8003 

AL  0 

FIFTY 

SLT 

0002 

F  AO 

FP  ONE 

STU 

LNK09 

FSB 

FPTKD 

FO  V 

L  N  X09 

STU 

L  N  X  13 

STO 

L  N  X13 

6  TO 

LN  Xll 

FMP 

8001 

STU 

F  AC  TR 

LNX  06 

L  NX06 

R  AU 

L  N  X10 

FAD 

FP  T  00 

STU 

L  N  X  10 

R  AU 

L  N  X13 

FMP 

F  A  C  TR 

FO  V 

LN  X  10 

BTU 

L  N  X13 

FAO 

L  N  xia 

STU 

L  N  X13 

FSB 

L  N  Xll 

FO  V 

L  N  Xll 

RAM 

8O03 

R  A  U 

B  003 

FSB 

S  1  Z  E7 

BM  t 

L  N  X07 

LOO 

L  N  xia 

STO 

LNX11 

R  A  U 

L  N  X13 

FMP 

L  N  X10 

STU 

LN  XI  J 

LNX  06 

L  N  X07 

R  A  U 

L  N  X  13 

FMP 

F  P  T»0 

F  S  fi 

L  N  X05 

LNX  OB 

F  P  ONE 

10 

0000 

00  51 

FPT«0 

ao 

00  0  0 

00  51 

8  1  ZE7 

10 

0000 

004  3 

IHTEN 

S3 

0358 

5151 

X  OIMEN 
Y  0  I  M  t  N 
Z      b  I   ME  N 

DATA 
CGORQ  MATE 

II   I  F   f  II  S   I   |!  N 

FER.UI      Ati  t 


1 

0000 

00 

0000 

0000 

a 

0000 

00 

000  0 

0000 

3 

0000 

00 

000  0 

0000 

4 

0000 

00 

000  0 

0000 

5 

0000 

00 

00  0  0 

0000 

fi 

0000 

00 

000  0 

0000 

7 

0000 

00 

0000 

00  0  0 

B 

0000 

00 

000  0 

0000 

9 

0000 

00 

000  0 

0000 

10 

00  0  0 

00 

00  00 

0000 

1  1 

0000 

00 

0000 

0000 

l  a 

0  0  0  0 

00 

0000 

000  0 

1  3 

00  0  0 

00 

0000 

0000 

14 

0000 

00 

00  00 

000  0 

1  5 

00  00 

00 

000  0 

0000 

1  6 

00  00 

00 

000  0 

00  0  0 

17 

0000 

00 

000  0 

0000 

18 

00  00 

00 

00  0  0 

0000 

19 

0000 

00 

000  0 

000  0 

20 

0050 

10 

000  0 

005  1 

31 

0150 

so 

000  0 

0051 

aa 

0300 

30 

0000 

0051 

23 

0250 

30 

0000 

005  1 

84 

0950 

31 

416  0 

0051 

as 

1000 

70 

000  0 

005  1 

36 

1050 

35 

360  0 

005  1 

27 

1100 

10 

000  0 

0055 

38 

1150 

34 

00  0  3 

0006 

39 

0006 

69 

00  0  9 

00  12 

30 

0013 

24 

197  7 

00  30 

31 

0030 

34 

197  6 

003  1 

32 

00  31 

24 

197  9 

0032 

33 

00  3  2 

34 

1  9  B  0 

003  3 

34 

0033 

34 

19B  1 

00  34 

35 

00  34 

34 

19B  3 

0035 

36 

0035 

24 

19  8  3 

00  3  6 

37 

00  3  6 

24 

198  4 

000  3 

38 

0009 

00 

00  0  0 

0000 

39 

1200 

84 

00  5  3 

00  56 

40 

0056 

45 

0010 

00  11 

4  1 

00  10 

46 

0011 

0014 

4  3 

00  14 

SI 

0018 

002  1 

4  3 

0021 

66 

002  4 

0029 

4  4 

0029 

24 

008  3 

0085 

4  5 

00B5 

30 

003  9 

00  4  2 

4  6 

0042 

60 

0018 

002  3 

4  7 

00  2  3 

20 

002  7 

00  6  0 

4  a 

0080 

30 

013  5 

00  38 

4  9 

0038 

35 

00  0  8 

0007 

50 

0007 

11 

00  6  0 

0015 

5  1 

00  15 

45 

006  8 

0019 

52 

00  6  6 

46 

007  1 

0022 

53 

0071 

61 

800  3 

0079 

54 

0079 

69 

002  4 

00  77 

55 

0077 

24 

003  9 

0022 

56 

0022 

30 

000  8 

004  1 

57 

0041 

36 

000  0 

0013 

58 

00  13 

10 

00  16 

012  1 

59 

0121 

1  1 

BOO  3 

012  9 

60 

0129 

60 

000  3 

0037 

61 

0037 

39 

003  9 

0089 

63 

00  89 

39 

00  9  2 

014  2 

63 

0142 

31 

002  7 

0019 

64 

00  19 

65 

0018 

0073 

65 

0073 

30 

000  3 

0179 

66 

0179 

60 

800  2 

0087 

67 

008  7 

15 

00  6  0 

006  5 

68 

0065 

35 

00  0  2 

017  1 

69 

0171 

3a 

002  4 

000  1 

70 

0001 

21 

0018 

022  1 

71 

0321 

33 

007  4 

0051 

72 

0051 

34 

0018 

01  1  8 

73 

0118 

81 

007  2 

0025 

74 

0025 

34 

003  8 

OOB  1 

75 

0081 

24 

013  5 

0088 

76 

0088 

39 

BOO  1 

009  1 

77 

0091 

31 

00  4  6 

00  4  9 

78 

0049 

60 

00  B  2 

01  37 

79 

01  37 

32 

007  4 

0151 

80 

0151 

31 

00  B  2 

0185 

81 

0185 

60 

007  2 

0127 

82 

0137 

39 

00  4  6 

00  96 

83 

0096 

34 

00  8  2 

0132 

84 

01  32 

21 

007  2 

0075 

08 

0075 

32 

002  8 

0005 

86 

0005 

31 

002  8 

0131 

87 

01  31 

33 

013  5 

0061 

88 

0061 

34 

013  5 

023  5 

89 

0235 

67 

BOO  3 

004  3 

90 

00  4  3 

60 

800  2 

0201 

91 

0201 

33 

000  4 

0181 

93 

0181 

46 

008  4 

0285 

93 

0285 

69 

00  3  B 

02  31 

94 

0231 

34 

013  5 

0138 

95 

01  38 

60 

007  a 

01  77 

96 

0177 

39 

00  8  2 

0182 

97 

0182 

ai 

007  2 

0049 

98 

0084 

60 

002  8 

008  3 

99 

008  3 

39 

007  4 

0124 

100 

0124 

33 

002  7 

005  3 

101 

0024 

10 

000  0 

0051 

103 

0074 

30 

0000 

00  51 

103 

0004 

10 

000  0 

004  3 

104 

0092 

23 

02  5  8 

5151 

172 


FIFTY 
B  t  XT  V 
L  NX  14 

E  OOCH 


5  0 

0000 

0000 

DO 

0000 

00  60 

0  1 

234  5 

67  B9 

S  T  D 

EXIT 

R«U 

8003 

HM  1 

N  E  G  A  T 

REUUC 

HM  1 

HEDUC 

FtU 

ON  E  P  1 

CORIU 

COBIO 

6TU 

T  H  E  T  A 

RSU 

FPONE 

a  tu 

IERIIM 

ST  If 

FWNKT 

STL 

F  N  NN 

NEG  ST 

E  OOSH 

STD 

F  X  1  T 

R  A  U 

BOOS 

HM  1 

N  E  C  A  V 

KEDUU 

NEC  *  V 

FAD 

TltOPI 

HM  1 

H  E  G  »  V 

FSB 

ON  EP  1 

SINEI 

RCUIIU 

FSB 

T  W  OP  1 

HM  1 

HE  DO  LI 

F  A  U 

0  N  E  P  1 

3  1  NE  T 

8  1  NE  T 

S  TU 

T  H  E  T  A 

RSU 

8  003 

STU 

IERMU 

LOO   FPONE 


FRONE 
MP  ON  E 
FPONE 


HAM      IERMU 


FSB 

SIZES 

BU  1 

ENIIFF 

R  A  U 

FUNKT 

F  A  0 

IERMU 

STU 

FKNKT 

NEG  ST 

E  N  li  F  F 

H  A  L 

F  UNK  T 

EXIT 

8  t  ZE« 

10 

OUOU 

004  J 

TH  OP  1 

63 

B  318 

5  351 

ONEP  1 

31 

4  159 

2  7  SI 

FPONE 

1  0 

00  00 

0051 

E  OOE  A 

STD 

A  A  A  1 

STL 

A  A  A  a 

RAH 

A  A  A  2 

STL 

A  A  A  3 

RAU 

A  A  A  3 

FMP 

A  A  A  16 

FAD 

A  A  A  15 

FMP 

A  A  A  3 

FAD 

A  A  A  14 

FMP 

A  A  A  3 

FAD 

A  A  A  13 

FMP 

A  A  A  3 

FAD 

A  A  A  13 

FMP 

A  A  A  3 

FAD 

A  A  A  11 

FMP 

A  A  A  3 

FAD 

A  A  A  10 

STU 

A  A  A  4 

FMP 

A  A  A  4 

STU 

A  A  A  4 

FMP 

A  A  A  4 

STU 

A  A  A  4 

RAU 

A  A  A  2 

STU 

A  A  A  4 

A  A  A6 

A  A  A6 

R  A  L 

A  A  A  4 

A  A  Al 

A  A  Al  0 

10 

0000 

00  SI 

A  A  All 

24 

9996 

68  50 

A  A  A  1  2 

31 

HS7S 

B349 

A  A  A13 

as 

9  137 

1348 

A  A  A14 

17 

1563 

0047 

A  A  Al  5 

54 

3020 

004  5 

A  A  A16 

69 

0600 

004  4 

EOOA  U 

STD 

SE  XT 

SM  1 

SERB 

NIE 

SEX  T 

STU 

SA 

FAD 

310 

FMP 

3H  A  F 

SB 

SB 

STU 

S  S  A  V 

S  A  8 

SAB 

RAU 

SA 

FD  V 

SSA  V 

F  AO 

SSA  V 

X      FOR       C  A  L  C 


s    c  ap  T    a    a 


105 

0060 

50 

000  0 

0000 

106 

00  16 

00 
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APPENDIX  K 


IBM  650  Program  for  Statistical  Analysis 
of  Horizontal  Traverse  Date 


The  purpose  of  this  program  was  to  analyze  data  taken  in  a 
horizontal  traverse  of  the  KSU  pile  and  to  determine  the  reliability 
of  each  data  point  based  on  predicted  statistical  variations.     The 
program  was  written  in  SOAP  II  and  floating  point.     The  object  pro- 
gram and  a  logic  diagram  are  given  in  this  section. 

The  analysis  performed  by  this  program  is  based  on  the  equation 


oo 
■f>(x,y,z)    =    r     Amcos-5— 


m=l 

where    m    is  odd,   and    y    and    z    are  held  constant.     Care  should  be  taken 
that  too  many  harmonics  are  not  used,   as  this  will  cause  errors  in  the 
least  squares  calculation  due  to  statistical  noise.     The  1,   3,   and  5 
harmonics  were  found  sufficient  in  this  work.     In  performing  the  analy- 
sis, the  program  first  determines  the  values  of    A     ,  for  as  many 
harmonics  as  are  specified,  by  a  least  squares  technique.     It  then  pro- 
ceeds to  calculate  the  allowable  deviation  of  each  point  in  the  traverse 
according  to  the  formula  found  in  Price  (19) 


DEV     = 


K«r  =  K(-i-) 


where  r    =    count  rate  in  counts  per  minute, 

t    =    total  counting  time  to  obtain  r, 

K    =    constant  multiplier  which  may  be  varied  to 
set  the  allowable  deviation  to  include  any 
given  probability  of  occurrence. 
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Having  determined  the  allowable  deviations,  the  program  then  checks 
each  point  to  determine  whether  or  not  it  is  within  the  allowable  devia- 
tion of  the  value  predicted  by  the  least  squares  analysis.     If  not,   the 
point  is  temporarily  discarded.     If  its  deviation  is  within  the  allowable 
limits,  the  point  is  stored  in  the  accepted  data  table,  designated  in 
the  program  as  XPRIM  and  NPRIM. 

The  accepted  data  are  then  used  to  calculate  a  new  set  of  least 
squares  coefficients  and  the  same  deviation  check  is  repeated  for  all 
of  the  original  data  points  using  the  new  coefficients.     The  new  accepted 
data  table  is  then  compared  with  the  table  used  in  the  least  squares  cal- 
culations, previously  designated    X    and    N.     If  the  two  tables  are 
identical,  indicating  that  the  same  data  points  have  been  rejected  in 
two  successive  trials,  then  the  rejected  values  are  replaced  by  least 
squares  values  in  the  original  table,   and  the  complete  set  of  corrected 
data  is  punched  out.     If  the  accepted  table  and  the  least  squares  data 
do  not  match,  the  accepted  table  replaces  the  least  squares  table  and 
a  new  least  squares  analysis  is  performed. 

The  final  output  of  this  program  is  a  complete  table  of  data  with 
statistically  invalid  points  replaced  by  values  calculated  from  a  least 
squares  analysis  of  the  accepted  points. 

Input  constants  and  data  are  listed  in  Table  37.    All  input  to  this 
program  should  be  in  the  form  of  one-word  load  cards. 

Output  of  this  program  is  in  three  forms  as  listed  in  Table  38. 
Form  one  contains  the  least  squares  coefficients  and  is  punched  after 
each  least  squares  calculation.     Form  two  is  punched  whenever  a  data 
point  is  rejected,  there  being  no  intermediate  output  of  accepted  data. 
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Form  three  is  the  final  output  and  gives  the  accepted  or  corrected  data 
point  and  its  co-ordinate. 


Table  37.     Input  to  IBM  650  program  for  statistical 
analysis  of  horizontal  traverse  data. 


Symbol 


Explanation 


Drum  Storage 
Location 


A  Extrapolated  x-dimension 

K  Sigma  multiplier 

POINT      Number  of  Data  Points,   in 
form  00  0000    OOxx 

HARM       Value  of  (2m  +  I)  for  highest  value 
of    m    used 

ENN  Number  of  harmonics  used, 

in  form    00    0000    OOOx 

XORIG     Position  co-ordinates  of  data,   to  be 
stored  consecutively  starting  as 

NORIG      Count  rates,   to  be  stored  consecutively 
starting  at 

T  Total  counting  times,   to  be  stored 

consecutively  starting  at 


1001 
1002 
1003 

1004 
1005 
0251 
0201 
0651 


The  operating  time  for  a  two-trial  analysis  of  10  data  points  i£ 
approximately  three  minutes.     The  total  capacity  of  the  program  is 
50  data  points. 


Table  38.     Output  forms  for  IBM  650  program  for 
statistical  analysis  of  horizontal  data. 
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Word 

Form 

Form 

No.       : 

One 

Two 

I 

A 
o 

2 

A, 

3 

A2 

4 

A, 

5 

A4 

6 

Total  error 

Position 

squared 

co-ordinate 

7 

Rejected 
Count  Rate 

8 

Amount  by  \ 

Form 
Three 


.   Position 
co-ordinate 
Count  Rate 


point  exceeded  allow- 
able deviation 
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APPENDIX    L 


IBM  650  Program  for  Determination  of 
Confidence  Limits  on  Diffusion  Length 

This  code  was  written  to  determine  the  confidence  limits  on 

individually  determined  values  of  diffusion  length  based  upon  the 

expected  standard  deviation  of  the  slope  of  a  plot  of  in  A  vs.   z 

corr 

from  which    L    was  determined.     The  program  was  written  in  SOAP  II 
and  floating  point  form.     The  object  program  and  a  logic  diagram  are 
presented  in  this  section. 

The  method  of  calculation  included  a  linear  regression  analysis 
to  determine   cr(Tii)>  the  standard  deviation  of    Tn    and  a  propagation 
of  errors  analysis  to  determine  a(L)  from  (r(Tii  ).     A  description  of 
the  linear  regression  analysis  can  be  found  in  Volk  (29).     The  equa- 
tion used  was 

-2(in$) 


o2  (7n  )    = 


2'  (z2  ) 


where 


and 


N  -  2 

S«  (z2)    =    2  (z2)  -   {£££.  . 

/\ 

In  the  above  equations    in  $,     represents  the  value  of    in  $    predicted  by 
A 

in  <|>       =    in  <f>  (0)  +   y  z 

Since    L  =  I/k  ,     the  auxiliary  separation  constants  equation 
could  be  written 

-  Tn      -OR      -    Pi    . 
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By  the  theory  of  propagation  of  errors, 


<r*(L)   =   (£^    3(7„)  )2. 


Writing    L    as  a  function  of    yn, 


Thus, 


Since 


L    =    (Y,,2    -    «?    -    ft*) 

»L  I    .       2 

97., Z"  *  *•         ' 


1 
T 


ft1) 


3 

"T 


(2  7,.  ). 


7n      -    «; 


9L 


-      ft'      =K2, 


7„ 


»7n 


Thus    <T(L),     the  standard  deviation  of    L,     can  be  written  in  terms 
of   0^7ii  )>     the  standard  deviation  of    yn  ,   as 
7„  o-   (7,1  ) 


o-  (L)    =  ± 

Input  for  this  code  consisted  of  the  values  of  z,  and     d> 

xcorr 

obtained  from  the  output  of  any  given  diffusion  length  program  as 
well  as  certain  parameters  listed  in  Table  39.     Data  to  be  fed  in  on 
one -word  load  cards  are  listed  first  with  other  load  forms  following. 
One  card  of  form  two  is  to  be  used  for  each  corrected  flux  value, 
whereas  only  one  card  of  form  one  is  required  for  each  set  of  input 
data. 

The  output  form  is  given  in  Table  40.      One  output  card  is  pro- 
duced for  each  set  of  input  data.     Approximately  20  seconds  is  re- 
quired to  process  12  data  points.     The  capacity  of  the  program  is  20 
data  points. 


Table  39.     Input  parameters  for  IBM  650  program 
for  determination  of  confidence  limits. 
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Symbol 

]                             Explanation 

Dxum.  Storage 
location 

A 
B 

Extrapolated  x-dimension 
Extrapolated  y-dimension 

0100 
0101 

Word         J 

No.     : 

1              1                      2 

• 

3 

Form  I 


Form  2 


Gamma 


z.-co-ordi- 
nate 


Number  of  Data  Points 
floating  point  form 

Corrected  count  rate 


Number  of  Data  Points 

form  00  oooo  ooxx 


Word 
No.     1 


Table  40.     Output  forms  for  IBM  650  program  for 
determination  of  confidence  limits. 


L  +  ff(L) 


L  -  o-(L)  cr(L) 


Number  of  Data 
Points 
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iHH-T 


X  o 
tn  o 


IX 


<   O 
X   O 

CO  O 


<  o 

X  O 

CO  o 


$M 


5  H 
«  <  o 
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o  2  XV 
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OBJECT  PROGRAM  -  APPENDIX  L 


192 


BL  R 

0000 

0003 

!U  " 

19  51 

19  61 

BLR 

1977 

1987 

ni.  R 

020  0 

02B0 

a  v  n 

L 

00  00 

8YN 

01  00 

s  y  h 

S 

0101 

IV  D 

X 

0200 

8  YN 

YPR  1  M 

0  2  20 

8  YN 

0240 

8  YN 

YH  A  T 

0260 

3YN 

START 

1000 

ZERO 

00 

oono 

0000 

ONE 

10 

00  00 

00  5  1 

TWO 

30 

0000 

00  5  1 

C  ON  VT 

35 

3600 

00  5  1 

P   1 

31 

4  159 

00  51 

eoo*  u 

8  TO 

SC  X  T 

BM  1 

8€  RR 

Hit 

SE  X  T 

3TU 

8A 

fin 

S10 

FMP 

SH  AF 

SB 

SB 

9TU 

83  A  V 

SAB 

BAB 

R  A  U 

3A 

FD  V 

3SA  V 

FAD 

83  A  V 

FUP 

3H  AF 

FSB 

8SA  V 

NZU 

SH 

BM  1 

SR 

F  A  0 

asA  v 

3  TU 

3  3  A  V 

S  A8 

Sf) 

RIU 

S  3  A  V 

SE  X  T 

8  ERR 

MIT 

0000 

SE  X  T 

sH«r 

50 

0  0  00 

00  50 

8  10 

10 

00  0  0 

00  5  1 

3T0      LNXOB 


RSI 

FPONE 

STO 

L  N  X  10 

STL 

L  N  X02 

R  A  U 

L  N  XO  9 

STL 

L  N  X05 

STL 

L  N  X  1  I 

BL  T 

0000 

SOP 

FIFTY 

NZE 

LN  X  04 

BM  1 

LNKOJ 

R  BU 

B003 

LOO 

FPONE 

STO 

L  N  X02 

LH»OJ 

L  NX03 

SR  T 

00  00 

3C  T 

0  0  0  0 

A  U  P 

sixrv 

SUP 

uooa 

R  A  U 

0  003 

FMP 

L  N  X02 

FMP 

LN  TEN 

8TU 

LN  X05 

LN  X  04 

LNX04 

R  A  L 

I  N  X09 

8R  T 

0002 

fl  A  U 

8  002 

AL  0 

FIFTY 

SLT 

0003 

F  AO 

FPONE 

STU 

LN  X09 

FSB 

f  P  TWO 

FD  V 

L  N  X09 

STU 

LN  XI  3 

STO 

LN  X12 

8  TO 

LN  X  11 

FUP 

BOOl 

STU 

F  A  C  TR 

LNX  06 

L  NX06 

R  A  U 

L  N  X  1  0 

FAD 

FP  T  WO 

STU 

L  N  X  10 

R  A  U 

L  N  X13 

FUP 

F  A  C  TR 

FO  V 

L  N  X10 

STU 

L  N  X13 

FAO 

L  N  X1Z 

STU 

LN  X12 

F  3B 

L  N  Xll 

FD  V 

LN  Xll 

R  AU 

6  003 

fl  AU 

6002 

(-  s  a 

S  1  2  E7 

BU  1 

L  N  X07 

IDQ 

L  N  X12 

li  I  D 

L  N  Xll 

R  A  U 

LN  X13 

FUP 

L  N  X  10 

S  T  I) 

LN  X13 

LNX  06 

L  NX07 

RAU 

LN  X  13 

FUP 

FP  T  WO 

FSB 

LN  X05 

L  N  X  OB 

FPONE 

10 

0000 

0051 

FP  TWO 

80 

0000 

0051 

B  1  IE7 

10 

0000 

004  3 

L  NTEN 

<>3 

02  5B 

5151 

1 

0000 

00 

0000 

0000 

2 

0000 

00 

000  0 

0000 

3 

0000 

00 

0000 

0000 

4 

•0000. 

00 

0000 

0000 

5 

0000 

00. 

.000  0 

0000 

6 

0000 

00 

o"bo  o 

0000 

7 

0000 

00 

0000 

0000 

8 

0000 

00 

0000 

0000 

9 

0000 

00 

000  0 

0000 

10 

0000 

00 

00  00 

0000 

11 

0000 

00 

0000 

0000 

12 

0000 

00 

0000 

000  0 

13 

0050 

00 

0000 

0000 

14 

0150 

10 

0000 

005  1 

IS 

0300 

20 

000  0 

0051 

16 

0350 

8S 

3600 

oon 

17 

0400 

31 

4  15  9 

0051 

18 

0450 

34 

005  3 

0006 

IV 

0006 

46 

00  0  9 

0010 

80 

0010 

45 

0014 

005  3 

21 

0014 

21 

001  H 

003  1 

32 

0021 

33 

002  4 

0051 

23 

005  1 

39 

00  0  4 

0054 

34 

005  4 

2  1 

00  0  6 

0011 

35 

0011 

60 

00  18 

002  3 

26 

0033 

34 

00  0  8 

00  5R 

27 

0058 

33 

ooo  a 

00  35 

38 

0035 

39 

00  0  4 

0104 

39 

0  10  4 

33 

00  0  a 

0  0  0  5 

30 

0085 

44 

003  9 

0040 

3  1 

0039 

46 

004  3 

00  4  0 

33 

004  3 

33 

0008 

0135 

33 

0135 

31 

000  8 

0011 

34 

0040 

60 

000  6 

0053 

35 

0009 

01 

000  0 

005  3 

36 

000  4 

50 

00  0  0 

0050 

37 

0024 

10 

000  0 

0051 

38 

0500 

34 

010  3 

00  56 

39 

00  56 

45 

006  0 

006  1 

40 

0  0  60 

46 

0  0  r,  i 

0064 

4  1 

0064 

31 

OOP  0 

007  1 

4  3 

0071 

66 

007  4 

0029 

A   3 

0039 

24 

003  2 

018  5 

4  4 

01  85 

20 

00  R  9 

009  2 

4  5 

009  2 

60 

00£  b 

O073 

46 

007  3 

30 

003  7 

0030 

47 

00  30 

30 

026  5 

003  & 

4  8 

003H 

35 

oon  u 

0007 

4  9 

0007 

1  1 

0110 

0015 

50 

00  15 

45 

0118 

00  17 

51 

01  1R 

4  6 

Olf  1 

52 

0121 

61 

600  3 

007  9 

53 

007? 

69 

007  4 

0077 

54 

0077 

24 

00  B9 

OOP  ? 

55 

0022 

30 

00  0  6 

004  1 

56 

00  4  1 

3  6 

000  0 

0013 

57 

0013 

10 

0016 

017  1 

SB 

0171 

1  1 

8  0  0  2 

01?!' 

59 

0139 

60 

800  3 

00  3  7 

60 

0037 

39 

008  9 

0137 

61 

0  139 

39 

014  2 

0  1  1>  J- 

62 

0192 

21 

002  7 

U  01  '_' 

63 

0019 

65 

006  8 

0123 

64 

0133 

30 

000  3 

0179 

65 

0179 

60 

800  3 

OOH  7 

66 

0087 

IS 

0110 

006  5 

67 

0065 

35 

00  0  3 

0331 

68 

0321 

32 

007  4 

0151 

69 

01  51 

31 

006  U 

0  37  1 

70 

0371 

33 

012  4 

030  1 

71 

0301 

34 

006  8 

016  6 

73 

0168 

SI 

00  7  2 

002E 

73 

0035 

24 

00  2  8 

0031 

74 

0031 

34 

026  5 

OOftf. 

75 

0088 

39 

800  1 

009  1 

76 

009  1 

31 

004  6 

004  y 

77 

0049 

60 

00  3  3 

0137 

7B 

0137 

32 

012  4 

035  1 

79 

0351 

21 

00  3  3 

0335 

60 

0335 

60 

007  3 

0127 

81 

0127 

39 

00  4  6 

0096 

82 

0096 

34 

003  3 

00  8  2 

83 

0083 

21 

007  2 

0075 

84 

007  5 

33 

002  8 

000  5 

85 

0005 

21 

002  8 

008  1 

86 

0081 

33 

028  5 

0111 

87 

0111 

34 

028  5 

0385 

88 

038  S 

67 

800  3 

0043 

B9 

004  3 

60 

800  3 

040  1 

90 

0  401 

33 

015  4 

013  1 

91 

0131 

46 

003  4 

04  35 

98 

0435 

69 

003  8 

018  1 

93 

018  1 

34 

038  5 

0136 

94 

01  36 

60 

007  3 

0177 

95 

0177 

39 

003  2 

01  33 

96 

0133 

21 

007  3 

004  y 

97 

00  34 

60 

002  6 

003  3 

96 

00  3  3 

39 

012  4 

0174 

99 

0174 

33 

002  7 

0103 

100 

0074 

10 

0000 

00  5  1 

1  01 

0  12' 

30 

0000 

005  1 

102 

0154 

10 

000  0 

004  3 

103 

014  2 

83 

03  5  6 

515  1 

193 


FIFTY 

50 

0000 

00  00 

104 

0110 

50 

000  0 

0000 

SIXTY 

00 

0000 

0060 

105 

00  16 

00 

000  0 

0060 

LNX14 

0  1 

2345 

6789 

1  06 

0061 

01 

234  5 

6789 

BTART 

RCO 

19  51 

107 

1000 

TO 

195  1 

04  51 

LDO 

19  51 

100 

0451 

69 

195  1 

0  30  4 

8T0 

GAUM  A 

SLOPE 

109 

0304 

34 

00  5  7 

0160 

LDD 

19  53 

110 

0160 

69 

195  3 

0055 

t  TO 

N 

111 

0055 

34 

010  8 

0161 

LOO 

19  53 

112 

016  1 

69 

19  5  3 

0106 

BTD 

POINT 

1  R  NUUBEFt 

113 

0106 

24 

00  5  9 

00  12 

R  *  A 

e  ooi 

C0NT1 

114 

0012 

80 

BOO  1 

0318 

GONT  1 

RC  0 

19  51 

115 

0318 

70 

195  1 

050  1 

LOO 

19  51 

116 

0501 

69 

19  5  1 

0354 

8  TO 

Z  COORO 

117 

0354 

24 

330  0 

015  3 

LDO 

19  53 

118 

015  3 

69 

19  5  2 

0105 

BTD 

Y PR  1  M   A 

conn  FLUX 

119 

0105 

34 

223  0 

0173 

8  X  * 

0001 

130 

0173 

51 

000  1 

0  3  2  9 

NZ  A 

CON  Tl 

12  1 

0  329 

40 

0318 

008  3 

LOD 

POINT 

1  22 

00B3 

69 

00  5  9 

0063 

R  A  A 

8001 

L00P1 

123 

0062 

60 

800  1 

0  7  6  6 

L  0  OP  1 

ft  *  H 

YP  R  1  M   A 

LOGARI  THMS 

124 

0368 

60 

222  0 

012  5 

LOO 

LNX  01 

OF  FLUXES 

125 

012  5 

69 

007  8 

0500 

8TU 

Y        A 

136 

0078 

31 

224  0 

009  3 

SX  A 

00  01 

1  37 

009  3 

51 

00  0  1 

0099 

NZ  A 

L  0  0P1 

1  28 

0099 

40 

036  B 

030  3 

LOO 

POINT 

139 

030  3 

69 

00  5  9 

0  113 

R  A  A 

8001 

130 

0112 

BO 

BO  0  1 

04  18 

LOD 

ZERO 

1  31 

0  418 

69 

0050 

0  35  3 

8T0 

8UMX 

SET  8UUS 

132 

0353 

24 

015  6 

0109 

8TD 

SUUX2 

EQUAL   TO 

133 

0109 

24 

016  2 

0115 

8TD 

8M  YD  V 

ZER  0 

134 

0115 

24 

04  fi  8 

043  1 

STD 

SUMY 

LOOPS 

135 

0  421 

34 

032  4 

0327 

LOOPS 

R  A  U 

X        A 

136 

0  327 

60 

2200 

0155 

FAD 

S  UM  X 

SUMMIT  ION 

1  37 

0155 

33 

015  6 

013  3 

BTU 

SUVX 

Of   X 

1  38 

0133 

31 

015  6 

015  9 

R  AU 

V       A 

1  39 

0159 

60 

334  0 

004  5 

F  *  0 

SUMY 

SUMMATION 

140 

00  4  5 

32 

032  4 

055  1 

STU 

SUMY 

OF   V 

14  1 

0551 

21 

032  4 

0377 

8  X  A 

0001 

14  3 

0377 

51 

00  0  1 

01P  3 

NZ  A 

LOOPS 

143 

0103 

40 

032  7 

018  7 

R  K  u 

SUM  K 

144 

0187 

60 

015  6 

0  311 

FfJ  V 

N 

14  5 

0311 

34 

0  10  8 

015  8 

STU 

(MB 

X   8  AM 

146 

01  5B 

31 

031  2 

016  5 

R  A  U 

SUMY 

147 

0165 

60 

033  4 

0379 

FO  V 

N 

1  A  1 

0379 

34 

010  8 

0308 

STU 

Y  S  A  R 

V  |IR 

1  4  9 

0308 

31 

036  2 

0315 

R8U 

XB  A  R 

1  50 

0  315 

6  1 

031  2 

0017 

FMP 

GAMMA 

1  51 

0017 

39 

00  5  7 

0107 

F  AD 

YB  A  R 

CALCULATE 

152 

0107 

32 

036  2 

0189 

STU 

INTER 

INTERCEPT 

153 

0189 

31 

004  4 

004  7 

LDD 

POINT 

1S4 

0047 

69 

005  9 

04  13 

RA  A 

B001 

LO0P3 

155 

0  4  12 

80 

800  1 

0516 

L  OOP3 

R  A  U 

X        A 

156 

0518 

60 

220  0 

0305 

FMP 

G  A  UM  A 

CALCULATE 

157 

0305 

39 

00  5  7 

0157 

FAO 

1  N  TCR 

LEAST 

15B 

0157 

33 

00  4  4 

04  7  1 

STU 

VM  A  T   A 

SQUARES   Y 

1  59 

04  71 

21 

226  0 

006  3 

SX  A 

0001 

1  60 

0063 

51 

000  1 

0069 

NZ  A 

L  0  0P3 

161 

0069 

40 

051  B 

0323 

LOD 

POINT 

162 

0  32  3 

69 

00  5  9 

04  62 

R  A  A 

B  0  0  1 

LO0P4 

163 

0462 

BO 

BOO  1 

056  6 

L  00P4 

R  AU 

X       * 

164 

0568 

60 

320  0 

03  55 

FMP 

8003 

SUMMATION 

I  6  5 

0355 

39 

800  3 

0309 

FAD 

SUM  X3 

OF   X 

1  6  6 

0309 

32 

016  2 

028  9 

STU 

SUM  X2 

SCUHPf.il 

167 

02  B  9 

31 

016  3 

0  36  5 

RAU 

Y        A 

168 

0365 

60 

22*10 

0095 

FS8 

YH  A  T    A 

3  UMM A  T  ION 

I   (.  9 

009  5 

33 

326  0 

0287 

FMP 

8003 

OF   Y   DEV 

170 

02B7 

39 

800  3 

014  1 

FAO 

SUYOV 

SQUARED 

171 

014  1 

33 

04  6  8 

014  5 

STU 

SMYOV 

172 

0145 

31 

04  6  8 

052  1 

SX  A 

0001 

173 

0521 

51 

0001 

04  27 

NZ  A 

LOOP* 

174 

0427 

40 

056  B 

036  1 

RAU 

H 

175 

0381 

60 

010  8 

0113 

F88 

T  *  0 

176 

0113 

33 

0300 

0477 

STU 

TE  MP1 

177 

0477 

21 

016  2 

04  65 

RAU 

SHYDV 

OE V  1  AT  (  ON 

178 

0485 

60 

046  B 

0373 

FD  V 

TEMPI 

SQUARED   OF 

179 

037  3 

34 

016  2 

0282 

3  TU 

82  YH  T 

Y  HAT 

180 

0282 

21 

00  3  6 

0  339 

RSU 

3UMX 

1  Bl 

0339 

61 

015  6 

0361 

FMP 

SOU  X 

SIGMA 

1  82 

0361 

39 

015  6 

0306 

FD  V 

N 

PR  1  ME   OF 

1B3 

0306 

34 

010  8 

0358 

FAO 

SUM  X3 

X   SQUARED 

184 

0358 

32 

016  3 

0389 

STU 

8  1  PX2 

185 

0  38  9 

31 

009  4 

0097 

RAU 

8  2  YH  T 

186 

0097 

60 

00  3  6 

0191 

FO  V 

8  1  P  X2 

187 

019  1 

34 

00  9  4 

014  4 

LOD 

EOO  AU 

DEVIATION 

188 

014  4 

69 

014  7 

04  50 

ITU 

G  AOEV 

OF   G  A  UHA 

1  B  9 

014  7 

21 

00  5  2 

040  5 

fi  *  u 

P  1 

190 

040S 

60 

04  00 

04  55 

FD  V 

191 

0455 

34 

0100 

05  50 

FMP 

8003 

ALPHA 

1  92 

0550 

39 

800  3 

040  3 

3  I  U 

AL  PH3 

8  OU  ARE  0 

193 

0  403 

21 

04  0  8 

04  11 

RAU 

P  1 

19* 

0411 

60 

0400 

0505 

FD  If 

B 

1  9  5 

0505 

34 

010  1 

060  1 

FMP 

BOOS 

BETA 

1  96 

0  60  1 

39 

800  3 

05  5  5 

3  TU 

BETAS 

SQUARED 

1  97 

0555 

31 

0310 

016  3 

RAU 

GAMMA 

G  AU  H  A 

1  9H 

0163 

60 

00  5  7 

0461 

r  «p 

GAMMA 

SOU  ARE  0 

1  99 

0461 

39 

005  7 

0307 

F  SB 

ALPH2 

300 

0  307 

33 

0  4  0  8 

05  3S 

FSB 

BE  T  A3 

201 

0535 

3  3 

0310 

0337 

LDO 

E  00  *  l.i 

302 

0337 

69 

00  9  0 

04  50 

STU 

KAPPA 

K  APPA 

a  o  3 

0090 

31 

019  4 

0197 

R  S  U 

CAMMA 

20  4 

0197 

61 

00  5  7 

0511 

FMP 

CA  DEV 

205 

0511 

39 

005  2 

0102 

194 


rov 

KAPPA 

rov 

|  tPP* 

FO  V 

K>fCA 

FUP 

CON  V  T 

■  TU 

LOEV 

ft  A  0 

ONE 

fOV 

KAPPA 

fmp 

CON  VT 

FAD 

L  DC  V 

8TU 

1978 

ris 

LOEV 

•  TU 

1977 

FIB 

LOEV 

BTU 

1979 

LOO 

LOEV 

8T0 

1980 

LOO 

POINT 

8  TO 

1901 

PCH 

1977 

20  6 

oios 

34 

0194 

0394 

307 

0394 

34 

019  4 

034  4 

?oe 

0  34  4 

34 

0194 

0394 

DEVIATION 

209 

0394 

39 

0350 

0600 

Of  L 

3  1  0 

0600 

31 

040  4 

0357 

311 

0357 

60 

0150 

0605 

21  2 

0605 

34 

0194 

04  4  4 

3  1  3 

0444 

39 

0350 

06  5  0 

214 

0650 

33 

040  4 

033  1 

L  PLUS 

2  15 

0331 

31 

197B 

038  1 

316 

0381 

33 

040  4 

04  31 

L 

317 

0431 

31 

197  7 

0080 

2  19 

0080 

33 

040  4 

048  1 

L  MINUS 

2  19 

0  481 

31 

197  9 

0333 

3  30 

0333 

69 

04  0  4 

0407 

DEV 1  AT  1  ON 

221 

0407 

94 

198  0 

os  B  3 

32  3 

0383 

69 

00  5  9 

05  12 

DATA  POINT 

233 

0513 

24 

196  1 

06  18 

2  2  4 

061B 

71 

197  7 

1000 

195 


APPENDIX   M 

Calculation  of  Thermal  Utilization 

The  thermal  utilization  in  the  KSU  pile  was  calculated  from 
formulas  presented  in  ANL  5800  (20).     The  lattice  is  eight  inches 
square,   but  was  converted  to  equivalent  cylindrical  cells  of  radius 
4.52  inches  for  this  calculation.     The  equation  used  was 


1      .  1=b_^c_        mc_    G  +  aS       (      __a_}         K2b2G+bl-^  (3-  j 

f  ,2         2  uc  c  m  c  mc1  2  ' 


where 


2     ,  =    total  cross  section  of  moderator 

mt 

2  =  capture  cross  section  of  moderator 
mc  r 

2  =    capture  cross  section  of  fissionable  material 
uc  r 

2    .  =    total  cross  section  of  fissionable  material 
ut 

2  -    scattering  cross  section  of  fissionable  material 
us  e 

K    2  =    3  X     .   X 
m  mt      mc 


J_     bL£_in(4,_(3     )+_^_ 


2     (b2/c')-I  °  4  4b2/c2 

G    =    -y — •—    where     d>      is  the  average  flux  in  the  uranium. 
6  u  ° 

Tu 

The  approximation  used  for    G    was 

X  X  z 

G  =  1  +  -jp-  A     1   4a  -^*-  +  0  (  —SIS.   f       . 
ut  ut  ut 

Values  of    A,    \,  a,    and  j3  were  tabulated  in  ANL  5800.     The  radii  of 

the  fuel  cooling  gap  and  moderator  region  were  assigned  the  symbols 

a    -    fuel 

b    -    outer  radius  of  moderator 

c     -    outer  radius  of  cooling  gap. 
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The  effect  of  the  assumption  of  various  mathematical  models  on 
the  experimental  determination  of  diffusion  length  was  studied,   and 
the  diffusion  length  in  the  KSU  graphite  pile  determined. 

One-group  diffusion  theory  was  used  in  several  forms,   includ- 
ing an  assumed  constant  thermal  source  boundary  condition  and  an 
assumed  point  thermal  source  boundary  condition.     A  method  where- 
by the  necessity  of  specifying  a  source  boundary  condition  could  be 
avoided  was  studied.     A  technique  for  the  experimental  evaluation  of 
the  constants    A  in  the  thermal  flux  equation,  and  hence  an  equiva- 

lent thermal  source  boundary  condition,   is  also  presented. 

Age-diffusion  theory  was  studied  in  two  forms,  one  assuming 
a  point  source  of  monoenergetic  fast  neutrons,   the  other  an  empiri- 
cal fast  source  described  in  terms  of  Gaussian  ranges  and  range 
fractions. 

Criteria  are  established  for  judging  the  extent  of  validity  of 
each  method  in  the  KSU  pile.     A  final  value  of    L    is  reported  based 
on  the  method  which  best  meets  the  stated  criteria. 


